Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,345,473 A→G 6.0% H29H (CAT→CAC yjdF ← DUF2238 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,345,4730AG 6.0% 183.9 / 2.0 67H29H (CAT→CACyjdFDUF2238 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (36/27);  new base G (2/2);  total (38/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

TGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCAAATCAGGGAAAATACTGACTATTCACCAGCAAGCGTGAAATTACAATCG  >  NC_000913/4345325‑4345620
                                                                                                                                                    |                                                                                                                                                   
tgtATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCa                                                                                                                                                     <  1:306877/149‑1 (MQ=255)
tgtATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCa                                                                                                                                                     <  2:208187/149‑1 (MQ=255)
     ggTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGc                                                                                                                                                <  1:87789/148‑1 (MQ=255)
            cgACCATCAGGAGGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATg                                                                                                                                         <  1:237777/149‑1 (MQ=255)
              aCCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCGTGGGCCGAAATGcc                                                                                                                                       >  1:245423/1‑149 (MQ=255)
              aCCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCGTGGGCCGAAATGcc                                                                                                                                       >  1:247504/1‑149 (MQ=255)
              aCCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACCTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGcc                                                                                                                                       <  1:670382/149‑1 (MQ=255)
                cATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCgg                                                                                                                                     >  1:792126/1‑149 (MQ=255)
                 aTCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGt                                                                                                                                    <  2:733109/149‑1 (MQ=255)
                  tCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGta                                                                                                                                   >  1:94143/1‑149 (MQ=255)
                  tCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGATTGTCATGGGCCGAAATGCCGGta                                                                                                                                   >  2:238953/1‑149 (MQ=255)
                         gatgGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCa                                                                                                                            <  2:787679/149‑1 (MQ=255)
                           tgGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAgg                                                                                                                          <  1:125362/149‑1 (MQ=255)
                           cggCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAgg                                                                                                                          <  1:125414/148‑1 (MQ=255)
                               gTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCGGCCAGGGGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTa                                                                                                                      >  2:690067/1‑149 (MQ=255)
                                 gAAGGCAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAAGCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAAc                                                                                                                    <  2:361783/149‑1 (MQ=255)
                                    ggAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGtt                                                                                                                 >  1:478898/1‑149 (MQ=255)
                                     gAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTa                                                                                                                >  2:698926/1‑149 (MQ=255)
                                       aaaTGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCTGAAATGCCGGTATAAATCAGGATTAACGTTAGc                                                                                                              <  1:514177/149‑1 (MQ=255)
                                              cGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTg                                                                                                       <  2:295336/149‑1 (MQ=255)
                                                 ataGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTgcgc                                                                                                    >  1:558394/1‑149 (MQ=255)
                                                   agagCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGTACGACAATAATCACCGGTGTCACTTTCATCAGCTAGTTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTg                                                                                                  <  1:507575/149‑1 (MQ=255)
                                                     agcagcGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAGATCAGGATTAACGTTAGCGTCAGTGCGCTGgt                                                                                                >  1:467571/1‑149 (MQ=255)
                                                       cagcGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCGGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGgtgt                                                                                              >  2:544588/1‑149 (MQ=255)
                                                         gcggcgTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTa                                                                                            >  1:115529/1‑149 (MQ=255)
                                                         gcggcgTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATAAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGGGCTGGTGTTa                                                                                            >  2:266093/1‑149 (MQ=255)
                                                          cggcgTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTaa                                                                                           >  2:645321/1‑149 (MQ=255)
                                                            gcgTTAACGGATAGCGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTGCATCAGCCAGGTGTGTTTGTCATGGGCCGAAATGCCAGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGa                                                                                         <  2:271999/149‑1 (MQ=255)
                                                             cgTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGaa                                                                                        >  1:688012/1‑149 (MQ=255)
                                                               ttAACGGATAACGTCTGGCGGGGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAAtt                                                                                      <  2:405069/149‑1 (MQ=255)
                                                                 aaCGGAGAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCATCGGTGTCACTTCCATCAGCCAGGTGAGTTTGCCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTaa                                                                                    <  1:687646/149‑1 (MQ=255)
                                                                   cGGATAACGTCTGGCGGTGGCAAGCAGTAGCTTCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTTATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACg                                                                                  <  2:447955/149‑1 (MQ=255)
                                                                   cGGATAACGTCGGTCGGTGGCAGGCAGTAGCTGCACGGCAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACg                                                                                  <  2:156976/149‑1 (MQ=255)
                                                                     gATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGc                                                                                <  2:486223/149‑1 (MQ=255)
                                                                     gATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGc                                                                                >  2:745167/1‑149 (MQ=255)
                                                                       tAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATGATCACCCGCGTCATTTCCATCACCCAGGTGAGTTTGTCATTGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCtt                                                                              <  2:86986/149‑1 (MQ=255)
                                                                       tAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCtt                                                                              <  2:781817/149‑1 (MQ=255)
                                                                        aaCGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATAAACGTTAGCGTCAGTGCGCCGGTGTTAAGACTGAAGGACTTg                                                                             >  2:620379/1‑149 (MQ=255)
                                                                             cTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGGTGGTGTTAAGAATTAACGGCTTGAgggt                                                                        >  1:229271/1‑149 (MQ=255)
                                                                                 cGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCGTGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCgg                                                                    <  2:245423/149‑1 (MQ=255)
                                                                                 cGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCGTGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCgg                                                                    <  2:247504/149‑1 (MQ=255)
                                                                                     ggCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCa                                                                >  2:801988/1‑149 (MQ=255)
                                                                                        aaGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATgg                                                             <  1:288955/149‑1 (MQ=255)
                                                                                         aGCAGTGGCTGCACGACAATAATCACCGGTGTCACTTCCAGCAGGCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGt                                                            <  2:410791/149‑1 (MQ=255)
                                                                                         aGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGt                                                            >  2:395931/1‑149 (MQ=255)
                                                                                                   gCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCaa                                                  >  1:366928/1‑149 (MQ=255)
                                                                                                    cACGACAGTAATCACCGGTATCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCaaa                                                 >  2:544050/1‑149 (MQ=255)
                                                                                                    cACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCaaa                                                 >  2:615429/1‑149 (MQ=255)
                                                                                                         cAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCAAATCAgg                                            >  2:581580/1‑149 (MQ=255)
                                                                                                               tCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCAAATCAGGGAAAAt                                      >  2:460885/1‑149 (MQ=255)
                                                                                                               tCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAATGTGCGGGTCACGGTGAGATCGCACATCAAGGAAAAt                                      >  2:423970/1‑149 (MQ=255)
                                                                                                                    ggTGTCACTTCCAGCAGCCAGGTGAGTTTGTCATGGGCCGAAATGGCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCAAATCAGGGAAAATACTGa                                 <  2:574403/149‑1 (MQ=255)
                                                                                                                           cTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCGGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCAACACAGGGAAAAGACTGACTATTCa                          >  2:649013/1‑149 (MQ=255)
                                                                                                                              ccATCAGCCAGGTGAGTTTGTCATGGGCCGAACTGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTAGTGTTAAGAATTAACGGCTTCAGTGTCCGGGTCATGGGGCGTTCGCAAATCACGGCAAATACTGACGATTccccc                       >  2:23229/1‑148 (MQ=255)
                                                                                                                              ccATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAATGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCCAATCAGGGAAAATACTAAGTATTCACCa                       >  2:350575/1‑149 (MQ=255)
                                                                                                                               cATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCAAATCAGGGAAAATACTGACTATTCACCAg                      >  1:78091/1‑149 (MQ=255)
                                                                                                                                 tCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTGATGGTGAGTTCGCAAATCAGGGAAAATACTGACTATTCACCAGCa                    >  2:58833/1‑149 (MQ=255)
                                                                                                                                     ccAGGTGAGATTGTCATGGGCCGAAATGTCGGTATAAATCAGGATTAACGTTAGCGTCAGAGCACTGGTGTGAAGAATTAGCGGCTTGAATGTGCGGGCCATGGTTAGGAGCTAAATAAGGGGACATACTCACTATACACCAGCAAGCg                >  2:380411/1‑149 (MQ=255)
                                                                                                                                        ggTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCGGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCAAATCAGGGCAAATACTGACTATTCCCCAGCAAGCGCGa             >  2:670237/1‑149 (MQ=255)
                                                                                                                                           gAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCAAATCAGGGAAAATACTGACTATTAACCAGCAAGCGTGAAAt          >  2:225596/1‑149 (MQ=255)
                                                                                                                                             gTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATTGTGAGTTCGCAAATCAGGGAAAATACTGACTCTTCACCAGCAAGCGTGAAATTa        >  1:31559/1‑149 (MQ=255)
                                                                                                                                             gTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCAAATCAGGGAAAATACTGACTATTCACCAGCAAGCGTGAAATTa        >  2:145178/1‑149 (MQ=255)
                                                                                                                                              tttGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCAAATCAGGGAAAATACTGACTATTCACCAGCAAGCGTGAAATTAc       >  1:298422/1‑149 (MQ=255)
                                                                                                                                              gttGTCAAGGGCCGAAAGGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCAAATCAGGGAAAATACTGACTATTCACCAGCAAGCGTGAAATTAc       <  1:396006/148‑1 (MQ=255)
                                                                                                                                              attGTCATGGACCCAAAGGCCGTAATAAATCAGGATTAACGTGGGCGTCAGGGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGGTCGCAAATCAGGGAAAATACTGACTATTCACCAGCAAGCGTGAAATTAc       >  2:670609/2‑149 (MQ=255)
                                                                                                                                                tGTCATGGGCGGAAATGCCGGTATAAATCAGGATTAACGTTAGAGTCGGTGCGCTGGTGTTAAGAATTACCGGCGTGAGCGTGCGGGTCATGGTGAGTTCGCAAATCAGGGAAAATACTGACTATTCACCAGCAAGCGTGAAATTACaa     <  2:142513/149‑1 (MQ=255)
                                                                                                                                                  tCATGGGCCGAAATGCCGGTATAAATCAGGATCAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCAAATCAGGGAAAATACTGACTATTCACCAGCAAGCGTGAAATTACAAt  <  1:652356/149‑2 (MQ=255)
                                                                                                                                                   cATGTGCCGAAATGCCGGTATACATCGGCATTAACGCTACCGTCTGGTCGCTCGTGTTAAGAATTATCGGCTTGAGTGTGCCGGTCATGGTGACTTCGCAAATCAGGGAAAATACTGACTATCCACCAGCAAGCGTGAAATTACAATCg  <  2:726279/149‑1 (MQ=255)
                                                                                                                                                   cATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCAAATCAGGGAAAATACTGACTATTCACCAGCAAGCGTGAAATTACAATCg  <  1:605560/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCAAATCAGGGAAAATACTGACTATTCACCAGCAAGCGTGAAATTACAATCG  >  NC_000913/4345325‑4345620

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: