Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 503,596 T→A 7.5% F208L (TTT→TTA) ‡ PP_0414 → polyamine ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947503,5960TA 7.5% 290.6 / 10.4 93F208L (TTT→TTA) ‡PP_0414polyamine ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (49/37);  new base A (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.39e-01

TCGGCGACGCTGCAGGGGTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGAC  >  NC_002947/503447‑503745
                                                                                                                                                     |                                                                                                                                                         
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                                                                                                                                             >  1:1444524/1‑116 (MQ=255)
                                               cAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAggtggtgg                                                                                                                                             <  2:1444524/116‑1 (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                                                                                                                >  1:1138788/1‑135 (MQ=255)
                                                         gcTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAg                                                                                                                <  2:1138788/135‑1 (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                                                                 <  2:2374331/143‑1 (MQ=255)
                                                                                             ccGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCt                                                                 >  1:2374331/1‑143 (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                                                                          >  1:1626597/1‑113 (MQ=255)
                                                                                                                  aTCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTcg                                                                          <  2:1626597/113‑1 (MQ=255)
                                                                                                                  aTCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGAcc                                     >  1:967996/1‑150 (MQ=255)
                                                                                                                  aTCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGAcc                                     >  2:1403970/1‑150 (MQ=255)
                                                                                                                     tcGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCgccgcggcaatggtc                                  >  1:922107/1‑143 (MQ=255)
                                                                                                                     tcGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCgccgcggcaatggtc                                  <  2:922107/150‑8 (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                                                                    >  2:970548/1‑108 (MQ=255)
                                                                                                                             ccTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGaa                                                                    <  1:970548/108‑1 (MQ=255)
                                                                                                                             ccTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGc                          <  2:568190/150‑1 (MQ=255)
                                                                                                                             ccTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGc                          >  2:2086356/1‑150 (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                                                                                   >  2:573411/1‑77 (MQ=255)
                                                                                                                                             aCCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAg                                                                                   <  1:573411/77‑1 (MQ=255)
                                                                                                                                              ccTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTtgc         >  1:1779153/1‑150 (MQ=255)
                                                                                                                                              ccTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTtgc         >  1:293305/1‑150 (MQ=255)
                                                                                                                                              ccTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTtgc         >  2:1475952/1‑150 (MQ=255)
                                                                                                                                               cTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTtgct        >  1:2047593/1‑150 (MQ=255)
                                                                                                                                               cTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTtgct        >  2:2028404/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)
                                                                                                                                                     tGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGAc  <  2:1852977/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                         
TCGGCGACGCTGCAGGGGTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACT‑‑GCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGAC  >  NC_002947/503447‑503745

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: