Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,601,825 G→C 5.7% P88A (CCC→GCC) ‡ ftsX ← putative ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,601,8250GC 5.7% 365.0 / 11.3 106P88A (CCC→GCC) ‡ftsXputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (47/53);  new base C (4/2);  total (51/55)
Fisher's exact test for biased strand distribution p-value = 4.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.84e-01

Only 100 of 111 total aligned reads displayed.

ACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGATTTACCGGTTTTGGCGAGGATTTC  >  NC_000913/3601677‑3601973
                                                                                                                                                    |                                                                                                                                                    
aCGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTggg                                                                                                                                                      >  1:602525/1‑149 (MQ=255)
aCGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTggg                                                                                                                                                      >  3:922833/1‑149 (MQ=255)
    cTTGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGc                                                                                                                                                  >  1:318983/1‑149 (MQ=255)
     ttGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACCGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCcgcg                                                                                                                                                 >  2:143623/1‑149 (MQ=255)
        cTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCa                                                                                                                                              <  4:522129/149‑1 (MQ=255)
          cGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCaga                                                                                                                                            <  1:1066808/149‑1 (MQ=255)
          cGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCaga                                                                                                                                            <  3:354802/149‑1 (MQ=255)
          cGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCaga                                                                                                                                            >  3:468209/1‑149 (MQ=255)
              cTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGCCCAGCGTctgttata                                                                                                                                        <  2:91978/149‑8 (MQ=255)
               tGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAAt                                                                                                                                       <  2:122200/149‑1 (MQ=255)
                 cAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATgg                                                                                                                                     <  4:12866/149‑1 (MQ=255)
                   aCTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCg                                                                                                                                   >  2:555027/1‑149 (MQ=255)
                    cTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGa                                                                                                                                  <  4:231104/149‑1 (MQ=255)
                    cTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGa                                                                                                                                  >  3:145229/1‑149 (MQ=255)
                     tgtgCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGAt                                                                                                                                 >  2:732517/1‑149 (MQ=255)
                     tgtgCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGAt                                                                                                                                 <  3:157633/149‑1 (MQ=255)
                       tgCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATaa                                                                                                                               <  1:156150/149‑1 (MQ=255)
                         ccaccaCGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGc                                                                                                                                                  >  3:843362/1‑128 (MQ=255)
                         ccaccaCGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGc                                                                                                                                                  <  4:843362/128‑1 (MQ=255)
                         ccaccaCGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAAcc                                                                                                                             >  2:885160/1‑149 (MQ=255)
                         ccaccaCGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAAcc                                                                                                                             >  1:1029099/1‑149 (MQ=255)
                          caccacGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCa                                                                                                                            <  3:70767/149‑1 (MQ=255)
                           accacGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCAt                                                                                                                           >  1:474261/1‑149 (MQ=255)
                            ccacGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGa                                                                                                                                  <  1:731894/141‑1 (MQ=255)
                            ccacGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGa                                                                                                                                  >  2:731894/1‑141 (MQ=255)
                              acGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCAc                                                                                                                        <  4:922833/149‑1 (MQ=255)
                                        cagcGTCATCGTCCAGCATTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCTTCGCCGCCTGGTTAACGTTTTTGAACACCAGATAACAGACGCTGGGTAGCGTCAGAGAAATGGCGATAACCTTCACCGTTAAAAAc                                                                                                              <  4:729158/149‑1 (MQ=255)
                                          gcGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGt                                                                                                            <  1:1036134/149‑1 (MQ=255)
                                          gcGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGt                                                                                                            <  1:116875/149‑1 (MQ=255)
                                           cGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTg                                                                                                           <  1:180611/149‑1 (MQ=255)
                                                tCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCACCTGGTTAACGTTTTTGTACGCCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGaa                                                                                                      <  3:316164/149‑1 (MQ=255)
                                                      aGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATGATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGttt                                                                                                <  4:912285/149‑1 (MQ=255)
                                                      aGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCaga                                                                                                                                            >  4:621511/1‑105 (MQ=255)
                                                      aGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCaga                                                                                                                                            <  3:621511/105‑1 (MQ=255)
                                                      aGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGttt                                                                                                <  4:1034336/149‑1 (MQ=255)
                                                       gCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTg                                                                                               >  1:551097/1‑149 (MQ=255)
                                                            ttttGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCtttt                                                                                          <  3:642833/149‑1 (MQ=255)
                                                            ttttGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCtttt                                                                                          <  4:577005/149‑1 (MQ=255)
                                                             tttGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTc                                                                                         <  4:110310/149‑1 (MQ=255)
                                                                  aGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATc                                                                                    >  2:354005/1‑149 (MQ=255)
                                                                  aGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATc                                                                                    >  4:32240/1‑149 (MQ=255)
                                                                    aTAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCt                                                                                  <  3:357057/149‑1 (MQ=255)
                                                                    aTAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCt                                                                                  >  3:992457/1‑149 (MQ=255)
                                                                     tAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTg                                                                                 <  1:986245/149‑1 (MQ=255)
                                                                      aaaCAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGc                                                                                >  1:562132/1‑149 (MQ=255)
                                                                        aCAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCaa                                                                              >  2:533595/1‑149 (MQ=255)
                                                                          aGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTAGACCATATAACAGACGCTGTGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGCCGAAGGGTTTGCTTTTCAGATCCTGCAATg                                                                            >  4:526641/1‑149 (MQ=255)
                                                                           gTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATgc                                                                           <  4:468209/149‑1 (MQ=255)
                                                                           gTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATgc                                                                           <  4:662235/149‑1 (MQ=255)
                                                                                       gACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGc                                                               <  1:555027/149‑1 (MQ=255)
                                                                                              aataCTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAgcgc                                                        >  2:789421/1‑149 (MQ=255)
                                                                                              aataCTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAgcgc                                                        >  1:283138/1‑149 (MQ=255)
                                                                                                taCTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCAc                                                      <  4:369847/149‑1 (MQ=255)
                                                                                                   tGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTg                                                   >  3:489587/1‑149 (MQ=255)
                                                                                                      gTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTc                                                >  3:944557/1‑149 (MQ=255)
                                                                                                       tCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACAGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGGACCTGTTCg                                               <  2:775396/149‑1 (MQ=255)
                                                                                                        cgccgcCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGt                                              <  1:533595/149‑1 (MQ=255)
                                                                                                         gccgccTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGtt                                             <  3:837559/149‑1 (MQ=255)
                                                                                                         gccgccTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGtt                                             >  1:837427/1‑149 (MQ=255)
                                                                                                            gccTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATTACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGaa                                          >  4:129511/1‑149 (MQ=255)
                                                                                                            gccTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGaa                                          >  2:810212/1‑149 (MQ=255)
                                                                                                             ccTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGaaa                                         <  2:559467/149‑1 (MQ=255)
                                                                                                              cTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTAAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGaaaa                                        >  4:318439/1‑149 (MQ=255)
                                                                                                                ggTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACg                                      >  2:20165/1‑149 (MQ=255)
                                                                                                                 gTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCAGCGTTAAAAACGTGGCGAACAGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGGGCACCTGTTCGTTGAAAACGt                                     >  3:647264/1‑149 (MQ=255)
                                                                                                                 gTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGt                                     <  2:1000786/149‑1 (MQ=255)
                                                                                                                 gTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGt                                     >  4:287718/1‑149 (MQ=255)
                                                                                                                 gTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGt                                     >  1:69963/1‑149 (MQ=255)
                                                                                                                 gTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTCGAAGGCATAGCGCACCTGTTCGTTGATAACGt                                     >  1:968826/1‑149 (MQ=255)
                                                                                                                  ttAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGGATAGCGCACCTGTTCGTTGAAAACGtt                                    >  1:912598/1‑149 (MQ=255)
                                                                                                                  ttAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGtt                                    >  1:346523/1‑149 (MQ=255)
                                                                                                                       gTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGttt                               >  4:606840/1‑149 (MQ=255)
                                                                                                                          tttGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGc                            <  1:895253/149‑1 (MQ=255)
                                                                                                                              tACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCtttt                                                                                          <  1:334863/83‑1 (MQ=255)
                                                                                                                              tACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCtttt                                                                                          >  2:334863/1‑83 (MQ=255)
                                                                                                                              tACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGAtt                        <  1:354005/149‑1 (MQ=255)
                                                                                                                                cacCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTc                                                >  4:766544/1‑123 (MQ=255)
                                                                                                                                cacCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTc                                                <  3:766544/123‑1 (MQ=255)
                                                                                                                                cacCAGATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGATTTa                      <  2:967294/149‑1 (MQ=255)
                                                                                                                                 acCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGATTTAc                     >  3:423630/1‑149 (MQ=255)
                                                                                                                                  ccATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGATTTAcc                    <  4:145229/149‑1 (MQ=255)
                                                                                                                                     tataACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGc                            >  2:896791/1‑138 (MQ=255)
                                                                                                                                     tataACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGc                            <  1:896791/138‑1 (MQ=255)
                                                                                                                                     tataACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGATTTACCGGt                 >  3:259737/1‑149 (MQ=255)
                                                                                                                                     tataACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGATTTACCGGt                 >  2:97442/1‑149 (MQ=255)
                                                                                                                                       taACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAATGTTGGTTTTGCGATTTACCGGttt               >  2:695181/1‑149 (MQ=255)
                                                                                                                                        ctctgACGCTGCCCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTCTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGATTTACCGGtttt              >  1:380360/5‑149 (MQ=255)
                                                                                                                                        aaCAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGATTTACCGGtttt              <  3:627702/149‑1 (MQ=255)
                                                                                                                                        aaCAGACGCTGGGCACCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGACGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGATTTACCGGtttt              <  2:311636/149‑1 (MQ=255)
                                                                                                                                         tctgACGCTGCCCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGATTTACCGGTTTTg             >  1:272586/4‑149 (MQ=255)
                                                                                                                                         aCAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTg                                   <  3:705178/127‑1 (MQ=255)
                                                                                                                                         aCAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGATTTACCGGTTTTg             <  1:885160/149‑1 (MQ=255)
                                                                                                                                         aCAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGATTTACCGGTTTTg             >  4:827925/1‑149 (MQ=255)
                                                                                                                                         aCAGACGCTGGGCAGCGTAAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTg                                   >  4:705178/1‑127 (MQ=255)
                                                                                                                                           tgACGCTGCCCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGATTTACCGGTTTTGGc           >  2:679108/2‑149 (MQ=255)
                                                                                                                                           aGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGATTTACCGGTTTTGGc           >  3:473050/1‑149 (MQ=255)
                                                                                                                                             aCGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGATTTACCGGTTTTGGCGa         <  4:259737/149‑1 (MQ=255)
                                                                                                                                              cGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGATTTACCGGTTTTGGCGAg        <  1:695181/149‑1 (MQ=255)
                                                                                                                                                    gCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGTGATTTACCGGTTTTGGCGAGGATTTc  <  4:896284/149‑1 (MQ=255)
                                                                                                                                                    gCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGATTTACCGGTTTTGGCGAGGATTTc  <  1:893466/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGATTTACCGGTTTTGGCGAGGATTTC  >  NC_000913/3601677‑3601973

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: