Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 503,596 T→A 11.3% 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,5960TA11.3% 176.3 / 12.9 62F208L (TTT→TTA) ‡PP_0414polyamine ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (33/22);  new base A (4/3);  total (37/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

GCGACGCTGCAGGGGTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTG  >  NC_002947/503450‑503740
                                                                                                                                                  |                                                                                                                                                  
gCGACGCTGCAGGGGTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGAc                                                                                                                                                 >  2:1183822/1‑150 (MQ=255)
       tGCAGGGGTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAggtgg                                                                                                                                          >  2:417785/1‑150 (MQ=255)
         cAGGGGTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAggtggtg                                                                                                                                        <  1:611497/150‑1 (MQ=255)
           ggggTTCAACTTCAATCTGGTGCGGGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTACCACCTCGTCAAACGAggtggcaaa                                                                                                                                      <  2:368450/150‑5 (MQ=255)
            gggTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGgtggtggtg                                                                                                                                     >  1:899173/1‑150 (MQ=255)
             ggTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGa                                                                                                                                    <  1:573297/150‑1 (MQ=255)
                  aaCTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCCCCTCTGCTGGCGTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTg                                                                                                                               <  1:1163764/150‑1 (MQ=255)
                   aCTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGt                                                                                                                              >  1:1110374/1‑150 (MQ=255)
                   aCTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTACCACCTCGTCAAACGAggtggc                                                                                                                                         >  2:196443/1‑138 (MQ=255)
                   aCTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTACCACCTCGTCAAACGAggtggc                                                                                                                                         <  1:196443/139‑2 (MQ=255)
                     tACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAggtggtaa                                                                                                                                       >  2:1035094/1‑137 (MQ=255)
                     tACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAggtggtaa                                                                                                                                       <  1:1035094/139‑3 (MQ=255)
                          tCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCgc                                                                                                                       >  2:192678/1‑150 (MQ=255)
                               tgcgtgcgGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCgccgggc                                                                                                                  <  1:417785/150‑1 (MQ=255)
                                          gcCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCa                                                                                                       >  2:327923/1‑150 (MQ=255)
                                           ccAGCGTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCaa                                                                                                      <  1:638004/150‑1 (MQ=255)
                                           ccAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCaa                                                                                                      <  2:1110374/150‑1 (MQ=255)
                                               ccTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACt                                                                                                  >  2:442855/1‑150 (MQ=255)
                                                 tGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGc                                                                                                <  2:183631/150‑1 (MQ=255)
                                                      gcTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCgg                                                                                           >  1:950913/1‑150 (MQ=255)
                                                           aCCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGAt                                                                                      >  1:214787/1‑150 (MQ=255)
                                                                tgctgGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAgg                                                                                                         <  1:320434/126‑1 (MQ=255)
                                                                tgctgGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAgg                                                                                                         >  2:320434/1‑126 (MQ=255)
                                                                tgctgGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCa                                                                                 >  2:1103434/1‑150 (MQ=255)
                                                                 gctgGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAg                                                                                <  2:1127020/150‑1 (MQ=255)
                                                                  ctgGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGc                                                                               <  2:899173/150‑1 (MQ=255)
                                                                   tgGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCg                                                                              >  2:1009014/1‑150 (MQ=255)
                                                                   tgGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCg                                                                              >  1:1030018/1‑150 (MQ=255)
                                                                            gccgGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTcgcg                                                                     >  1:1419538/1‑150 (MQ=255)
                                                                                ggTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGaa                                                               <  1:192678/150‑1 (MQ=255)
                                                                                     cccTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCa                                                                                 <  1:419365/129‑1 (MQ=255)
                                                                                     cccTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCa                                                                                 >  2:419365/1‑129 (MQ=255)
                                                                                      ccTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAg                                                                                                                                              >  1:561349/1‑67 (MQ=255)
                                                                                      ccTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAg                                                                                                                                              <  2:561349/67‑1 (MQ=255)
                                                                                       cTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGc                                                        >  2:1295374/1‑150 (MQ=255)
                                                                                          ccGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTc                                                                                                                         >  1:568140/1‑84 (MQ=255)
                                                                                          ccGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTc                                                                                                                         <  2:568140/84‑1 (MQ=255)
                                                                                                  tGCGCCCGGTCTCATCTCGGGAGCTCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGc                                             <  1:280038/150‑1 (MQ=255)
                                                                                                     gcCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCcgc                                          >  2:418768/1‑150 (MQ=255)
                                                                                                       ccGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCa                                                                                                           >  2:1273971/1‑85 (MQ=255)
                                                                                                       ccGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCa                                                                                                           <  1:1273971/85‑1 (MQ=255)
                                                                                                             tCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCga                                                                    <  2:1191175/118‑1 (MQ=255)
                                                                                                             tCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCga                                                                    >  1:1191175/1‑118 (MQ=255)
                                                                                                               aTCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGAcc                                >  2:694094/1‑150 (MQ=255)
                                                                                                                  tcGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTg                             >  1:327771/1‑150 (MQ=255)
                                                                                                                  tcGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACAAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTg                             >  1:762824/1‑150 (MQ=255)
                                                                                                                    gggAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCt                           <  2:1382448/150‑1 (MQ=255)
                                                                                                                     ggAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTg                          >  1:454716/1‑150 (MQ=255)
                                                                                                                        gCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGAtt                       <  2:214787/150‑1 (MQ=255)
                                                                                                                            tGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCTGACTATTGCTGCGGCGGGGACCTTGCTGATTGCGt                   >  2:74816/1‑150 (MQ=255)
                                                                                                                            tGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGt                   <  2:1030018/150‑1 (MQ=255)
                                                                                                                              tttgcctttgccACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTtc                 >  2:1394303/1‑150 (MQ=255)
                                                                                                                              tttgcctttgccACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGAGGTATCTGGAGTTTGCGCGGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTtc                 >  2:1106822/1‑150 (MQ=255)
                                                                                                                                 tcctttgccACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCg              <  1:1207699/149‑1 (MQ=255)
                                                                                                                                     ccaccACCTCGTCAAACGAGGTGGCGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGgt            <  1:419415/145‑1 (MQ=255)
                                                                                                                                     ccaccACCTCGTCAAACGAGGTGGCGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGgt            >  2:419415/4‑148 (MQ=255)
                                                                                                                                       gccACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGGACTTCCTGATTGCCTTCTTTGtgttt        >  2:227208/1‑149 (MQ=255)
                                                                                                                                       accACCTCGTCAAACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGtgttg        >  1:353658/2‑150 (MQ=255)
                                                                                                                                        ccACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGtgttgt       >  1:1011138/1‑150 (MQ=255)
                                                                                                                                        ccACCTCGTCAAACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGGAAGCGCAGCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGtgttgt       >  2:37983/1‑150 (MQ=255)
                                                                                                                                          aCCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGttgttg     >  2:458624/1‑150 (MQ=255)
                                                                                                                                             tCGTTTGACGAGGTGGTGGTGACGCTGTTATTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTtgctg  >  1:749556/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GCGACGCTGCAGGGGTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTG  >  NC_002947/503450‑503740

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