Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 503,595 T→C 7.5% F208S (TTT→TCT) ‡ PP_0414 → polyamine ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947503,5950TC 7.5% 272.0 / 8.3 93F208S (TTT→TCT) ‡PP_0414polyamine ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (49/37);  new base C (4/3);  total (53/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

TCGGCGACGCTGCAGGGGTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTG  >  NC_002947/503447‑503743
                                                                                                                                                    |                                                                                                                                                        
tCGGCGACGCTGCAGGGGTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGttt                                                                                                                                                         <  2:569141/150‑1 (MQ=255)
   gCGACGCTGCAGGGGTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGAc                                                                                                                                                     <  1:1246881/150‑1 (MQ=255)
   gCGACGCTGCAGGGGTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCGTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGAc                                                                                                                                                     >  2:1565860/1‑150 (MQ=255)
       cGCTGCAGGGGTTCAACTACGATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAgg                                                                                                                                                 <  1:418657/150‑1 (MQ=255)
          tGCAGGGGTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAggtgg                                                                                                                                              >  1:348006/1‑150 (MQ=255)
                ggTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCCCCTCGTTTGTCGAGGTGGTGGTGa                                                                                                                                        <  2:979015/150‑1 (MQ=255)
                     aaCTACAGTCTGGTGCGTGCGGCGGCCAGCCTGGGCGGTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTg                                                                                                                                   <  2:841020/150‑1 (MQ=255)
                     aaCTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGttt                                                                                                                                                         >  2:168229/1‑129 (MQ=255)
                     aaCTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGttt                                                                                                                                                         <  1:168229/129‑1 (MQ=255)
                     aaCTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTg                                                                                                                                   >  2:431657/1‑150 (MQ=255)
                      aCTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTCAAACGAggtggcaaa                                                                                                                                          <  2:2347391/142‑5 (MQ=255)
                      aCTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTCAAACGAggtggcaaa                                                                                                                                          >  1:2347391/1‑138 (MQ=255)
                      aCTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCATCCTTGTTTGACGAGGTGGTGGTGACGCTGt                                                                                                                                  <  1:1643599/150‑1 (MQ=255)
                        tACCATATGGTGCGCGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTTGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTc                                                                                                                                <  1:290048/150‑1 (MQ=255)
                         aCAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTcc                                                                                                                               >  1:1505992/1‑150 (MQ=255)
                           aaTCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGAc                                                                                                                                       >  2:2236097/1‑140 (MQ=255)
                           aaTCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGAc                                                                                                                                       <  1:2236097/140‑1 (MQ=255)
                             tCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGt                                                                                                                                  <  1:1556949/143‑1 (MQ=255)
                             tCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGt                                                                                                                                  >  2:1556949/1‑143 (MQ=255)
                             tCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTc                                                                                                                                <  2:1755796/145‑1 (MQ=255)
                             tCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTc                                                                                                                                >  1:1755796/1‑145 (MQ=255)
                              cTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCgcc                                                                                                                          >  1:1035739/1‑150 (MQ=255)
                                ggTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCgccgg                                                                                                                        <  1:2272389/150‑1 (MQ=255)
                                  tgcgtgcgGCGGCCAGCCTGGGCGCTTCACCTCTGCTAGCTTTCCGCCGGGTGACCCTGCCGCTGATCGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGGGGTGGTGACGCTGTTCCTCgccgggc                                                                                                                      >  1:1946948/1‑150 (MQ=255)
                                     gtgcgGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCgccgggccgg                                                                                                                   >  2:2369721/1‑150 (MQ=255)
                                          gcggcCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAg                                                                                                              <  2:348006/150‑1 (MQ=255)
                                          gcggcCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAg                                                                                                              <  1:60148/150‑1 (MQ=255)
                                              ccAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTc                                                                                                                             >  1:1094026/1‑131 (MQ=255)
                                              ccAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTc                                                                                                                             <  2:1094026/131‑1 (MQ=255)
                                               cAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAggtggtgg                                                                                                                                           <  2:1444524/116‑1 (MQ=255)
                                               cAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAggtggtgg                                                                                                                                           >  1:1444524/1‑116 (MQ=255)
                                                     gggCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGGTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑Gcc                                                                                                 >  1:1255334/1‑150 (MQ=255)
                                                       gcgcTTCACCTCTGCTGGCTTTCCGCAGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCcgc                                                                                               >  2:884353/1‑150 (MQ=255)
                                                         gttttaCCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCCCCTCGCTTGACGGGGTGGTCGTCAGCCTGTTCCTCGCCGGGCCGGAGCAGGCCACACT‑‑GCCGCgg                                                                                            <  2:84552/145‑1 (MQ=255)
                                                         gcTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAg                                                                                                              <  2:1138788/135‑1 (MQ=255)
                                                         gcTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAg                                                                                                              >  1:1138788/1‑135 (MQ=255)
                                                           ttCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCa                                                                                          >  1:205219/1‑150 (MQ=255)
                                                              aCCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGAt                                                                                       >  2:2350820/1‑150 (MQ=255)
                                                                ctctGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGt                                                                                     >  1:568190/1‑150 (MQ=255)
                                                                 tctGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGtt                                                                                    >  1:407327/1‑150 (MQ=255)
                                                                    gctgGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAg                                                                                 <  2:1035739/150‑1 (MQ=255)
                                                                     ctgGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGc                                                                                <  1:1407589/150‑1 (MQ=255)
                                                                        gCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCGTCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGTGCGAGGGGTGGTGACGCTGTTCGTCGCCGGGCCGGAGCCGGGAATACT‑‑GCTGCTGCAGATGTTCAGCGGt                                                                             >  1:1852977/1‑150 (MQ=255)
                                                                        gCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTCAAACGAggtgg                                                                                                                                              >  1:226901/1‑88 (MQ=255)
                                                                        gCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTCAAACGAggtgg                                                                                                                                              <  2:226901/88‑1 (MQ=255)
                                                                          tttCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTAt                                                                           >  1:867638/1‑150 (MQ=255)
                                                                                       aCCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTg                                                              >  2:896983/1‑150 (MQ=255)
                                                                                        cccTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGa                                                             <  1:2369721/150‑1 (MQ=255)
                                                                                             ccGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCt                                                               >  1:2374331/1‑143 (MQ=255)
                                                                                             ccGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCt                                                               <  2:2374331/143‑1 (MQ=255)
                                                                                              cGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGa                                                       >  2:1665816/1‑150 (MQ=255)
                                                                                              cGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGa                                                       >  1:468600/1‑150 (MQ=255)
                                                                                               gCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGAc                                                      <  2:1946948/150‑1 (MQ=255)
                                                                                                        gcCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCcgc                                             <  2:1255334/150‑1 (MQ=255)
                                                                                                             gtgtCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATATTCAGCGGTATTCGCGAAAACCTGAGCCCGACCATTGCCgcggcgg                                        >  2:1125078/1‑150 (MQ=255)
                                                                                                                  aTCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTcg                                                                        <  2:1626597/113‑1 (MQ=255)
                                                                                                                  aTCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTcg                                                                        >  1:1626597/1‑113 (MQ=255)
                                                                                                                  aTCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGAcc                                   >  2:1403970/1‑150 (MQ=255)
                                                                                                                  aTCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGAcc                                   >  1:967996/1‑150 (MQ=255)
                                                                                                                     tcGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCgccgcggcaatggtc                                <  2:922107/150‑8 (MQ=255)
                                                                                                                     tcGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCgccgcggcaatggtc                                >  1:922107/1‑143 (MQ=255)
                                                                                                                     tcGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATATTGAGCGGGACTCGCGAGAAGCTGAGCCCGACCATTGGAGCGGCGGCGACCTTg                                >  1:752175/1‑150 (MQ=255)
                                                                                                                      cGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGAGAACAGC‑‑ACCGCGTCAGATGTTCAGCGGTATTCGCCAGACCCTGAGCCCGACCATTGCCGCTGCGGCGACCTTGc                               >  1:1022393/1‑150 (MQ=255)
                                                                                                                       gggAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCt                              <  2:407327/150‑1 (MQ=255)
                                                                                                                        ggAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCcgc                                                                                               <  2:1319346/85‑1 (MQ=255)
                                                                                                                        ggAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCcgc                                                                                               >  1:1319346/1‑85 (MQ=255)
                                                                                                                             ccTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGATGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGc                        >  2:2125658/1‑150 (MQ=255)
                                                                                                                             ccTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGaa                                                                  <  1:970548/108‑1 (MQ=255)
                                                                                                                             ccTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGaa                                                                  >  2:970548/1‑108 (MQ=255)
                                                                                                                             ccTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGc                        >  2:2086356/1‑150 (MQ=255)
                                                                                                                             ccTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGc                        <  2:568190/150‑1 (MQ=255)
                                                                                                                               tGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGt                      >  1:2379446/1‑150 (MQ=255)
                                                                                                                                    gcctttgccACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCg                 <  1:1565860/150‑1 (MQ=255)
                                                                                                                                     tcaccaccACCTCGTCAAACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCgg                >  1:1967832/7‑150 (MQ=255)
                                                                                                                                       tttgccACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACt                                                                                                      <  2:630617/65‑1 (MQ=255)
                                                                                                                                       tttgccACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACt                                                                                                      >  1:630617/1‑65 (MQ=255)
                                                                                                                                          gccACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGtgttg           <  1:2359578/150‑1 (MQ=255)
                                                                                                                                           ccACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTTTGCCGGCGGCGACCTTGCTGATTGCGTGCTCGGtgttgt          >  1:633266/1‑150 (MQ=255)
                                                                                                                                           ccACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGc                               >  2:1929071/1‑129 (MQ=255)
                                                                                                                                           ccACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGc                               <  1:1929071/129‑1 (MQ=255)
                                                                                                                                           ccACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGtgttgt          >  2:1145623/1‑150 (MQ=255)
                                                                                                                                           ccACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGtgttgt          <  1:2058968/150‑1 (MQ=255)
                                                                                                                                           ccACCTCGTCAAACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGtgttgt          <  1:1794889/150‑1 (MQ=255)
                                                                                                                                            cACCTCGTCAAACGAGGTGGTGGTAACGGTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCGGCTGTTAAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGtgttgtg         >  2:548442/1‑149 (MQ=255)
                                                                                                                                             aCCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAg                                                                                 <  1:573411/77‑1 (MQ=255)
                                                                                                                                             aCCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAg                                                                                 >  2:573411/1‑77 (MQ=255)
                                                                                                                                              ccTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTtgc       >  1:293305/1‑150 (MQ=255)
                                                                                                                                              ccTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTtgc       >  1:1779153/1‑150 (MQ=255)
                                                                                                                                              ccTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTtgc       >  2:1475952/1‑150 (MQ=255)
                                                                                                                                               cTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTtgct      >  2:2028404/1‑150 (MQ=255)
                                                                                                                                               cTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTtgct      >  1:2047593/1‑150 (MQ=255)
                                                                                                                                                 cGTTTGACGCGGTGGTTGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACCCTAGGGTGCTGGCGAAG‑‑CTGGAGTAGTAGCGAGAACCGGAGAACGCCCATTGCCGCAGCGGCGGCATTGCTGATTGCGTTCTCGGTGTTGTtgctgc    <  1:2350820/150‑1 (MQ=255)
                                                                                                                                                 cGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTtgctgc    <  1:2086356/150‑1 (MQ=255)
                                                                                                                                                 cGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTtgctgc    <  2:967996/150‑1 (MQ=255)
                                                                                                                                                  gTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTgctgct   >  2:424496/1‑150 (MQ=255)
                                                                                                                                                  gTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCAGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTgctgct   >  2:66260/1‑150 (MQ=255)
                                                                                                                                                   tttGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGctgctg  >  1:2165125/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                        
TCGGCGACGCTGCAGGGGTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTG  >  NC_002947/503447‑503743

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