Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,389,957 C→G 6.5% P25A (CCG→GCG)  PP_1212 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,389,9570CG 6.5% 285.6 / 9.9 92P25A (CCG→GCG) PP_1212hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (42/44);  new base G (3/3);  total (45/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.78e-01

TGCTTGCTAACCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAGGCGGCGACAAGGCCGACTGA  >  NC_002947/1389811‑1390106
                                                                                                                                                  |                                                                                                                                                     
tgcttgctAACCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGATCAGCTGGGAAGTGCTGCAGAAGATCAGCgccgcgccgc                                                                                                                                                    >  1:1307229/1‑150 (MQ=255)
 gcttgctAACCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCt                                                                                                                                                   >  1:1174561/1‑150 (MQ=255)
 gcttgctAACCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCt                                                                                                                                                   <  2:2365052/150‑1 (MQ=255)
  cttgctAACCATATATAATCCCTCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTTCTGCAGAAGATCAGCGCCGCGCCGCTg                                                                                                                                                  >  2:1798924/1‑150 (MQ=255)
    tgctAACCATATATAATCCCGCTTTTTGCTGCCCCCCGACAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGCAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGAt                                                                                                                                                >  2:1265686/1‑150 (MQ=255)
     gctAACCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGGCGCGGCGCTGATc                                                                                                                                               <  1:1287978/150‑1 (MQ=255)
         aCCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAgcgccgcgc                                                                                                                                                       >  2:314933/1‑138 (MQ=255)
         aCCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAgcgccgcgc                                                                                                                                                       <  1:314933/138‑1 (MQ=255)
          ccATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTg                                                                                                                                          <  1:1063907/150‑1 (MQ=255)
           cATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGGCGCGGCGCTGATcttc                                                                                                                                            >  2:1559711/1‑144 (MQ=255)
           cATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGGCGCGGCGCTGATcttc                                                                                                                                            <  1:1559711/147‑4 (MQ=255)
            atagatAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTCGACGACACACACCCATGCTCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGcc                                                                                                                                        >  2:1064220/1‑150 (MQ=255)
              atataATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCg                                                                                                                                      <  2:1739751/150‑1 (MQ=255)
                ataATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGGCGCGGCGCTGATc                                                                                                                                               <  2:895128/139‑1 (MQ=255)
                ataATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGGCGCGGCGCTGATc                                                                                                                                               >  1:895128/1‑139 (MQ=255)
                    tCCCGCTTTTTGCTGCTCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGt                                                                                                                                >  2:1574546/1‑150 (MQ=255)
                      ccGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGcc                                                                                                                                   >  1:2310635/1‑145 (MQ=255)
                      ccGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGcc                                                                                                                                   <  2:2310635/145‑1 (MQ=255)
                      ccGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCa                                                                                                                              >  1:811401/1‑150 (MQ=255)
                      ccGCTTTTTGCTGACACAAGGAAGACAGAGTTATGGAAGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCa                                                                                                                              <  2:1307229/150‑1 (MQ=255)
                       cGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAg                                                                                                                             <  1:2164538/150‑1 (MQ=255)
                       cGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAg                                                                                                                             >  1:2098276/1‑150 (MQ=255)
                         cTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGgc                                                                                                                           <  2:811401/150‑1 (MQ=255)
                         cTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGgc                                                                                                                           >  1:1918572/1‑150 (MQ=255)
                         cTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGgc                                                                                                                           >  1:1918423/1‑150 (MQ=255)
                         cTTTTTCCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGCGCACCCCATGGAAGTGCTGCAGAAGATCACCCCCGCGCCGCTGCTCGATTGCCCGGCCTGTCAGgc                                                                                                                           >  2:1151762/1‑150 (MQ=255)
                                 tgcCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCAc                                                                                                                   >  2:2425509/1‑150 (MQ=255)
                                  gcCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACt                                                                                                                  >  1:2259028/1‑150 (MQ=255)
                                  gcCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACt                                                                                                                  <  2:2292932/150‑1 (MQ=255)
                                     ccccGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGaa                                                                                                               >  2:1767734/1‑150 (MQ=255)
                                     ccccGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGaa                                                                                                               >  1:337994/1‑150 (MQ=255)
                                              ccagcGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAgctgct                                                                                                      >  2:2212939/1‑150 (MQ=255)
                                              ccagcGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAgctgct                                                                                                      <  2:1419000/150‑1 (MQ=255)
                                               cagcGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGctgctg                                                                                                     >  2:465283/1‑150 (MQ=255)
                                               cagcGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGACTGGCCTGTCAGCCACCGGCACCGAAACACctgctg                                                                                                     >  2:1580415/1‑150 (MQ=255)
                                               cagcGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATAAATGTGCAACCTTCTAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGctgctg                                                                                                     >  2:585133/1‑150 (MQ=255)
                                                agcGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGt                                                                                                    >  2:2295074/1‑150 (MQ=255)
                                                    ttATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTTCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGt                                                                                                <  1:1798924/150‑1 (MQ=255)
                                                    ttATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGccc                                                                                                                                       >  1:1201672/1‑111 (MQ=255)
                                                    ttATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGccc                                                                                                                                       <  2:1201672/111‑1 (MQ=255)
                                                    ttATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGACCCCGCCGTTCCGCATGGAATTGGGGCAGGGGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGt                                                                                                >  1:1486853/1‑150 (MQ=255)
                                                        ggTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCa                                                                                            <  1:2212939/150‑1 (MQ=255)
                                                         gTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAg                                                                                           >  1:2399994/1‑150 (MQ=255)
                                                            gACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCt                                                                                        <  2:2399994/150‑1 (MQ=255)
                                                             aCGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCtt                                                                                       >  1:1611051/1‑150 (MQ=255)
                                                                   gaCACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAgg                                                                                          <  2:120902/141‑1 (MQ=255)
                                                                   gaCACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAgg                                                                                          >  1:120902/1‑141 (MQ=255)
                                                                   gaCACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCt                                                                                 >  2:287985/1‑150 (MQ=255)
                                                                   gaCACCCATGCCCCGTTATGACTATCAATGTGCGTCCTGCGACCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCt                                                                                 <  2:829051/150‑1 (MQ=255)
                                                                    acacCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCGGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCCGGCTTTCGCCTg                                                                                >  2:1238537/1‑150 (MQ=255)
                                                                      acCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCGGGGTGTTGCTAGCCCCTTCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAg                                                                              >  2:1241849/1‑150 (MQ=255)
                                                                      acCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAg                                                                              <  2:1918572/150‑1 (MQ=255)
                                                                      acCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAg                                                                              <  2:1918423/150‑1 (MQ=255)
                                                                        ccATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCg                                                                            >  2:2237550/1‑150 (MQ=255)
                                                                        ccATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCg                                                                            <  2:2225460/150‑1 (MQ=255)
                                                                             ccccTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAAc                                                                       >  1:586432/1‑150 (MQ=255)
                                                                               ccTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTACTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACgg                                                                     >  1:1300965/1‑150 (MQ=255)
                                                                                cTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGt                                                                    <  2:535931/150‑1 (MQ=255)
                                                                                  ttATGACTATCAATGTGCGTCCTGGGAGCACCGCATGGGAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCCCTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTg                                                                  <  2:77886/150‑1 (MQ=255)
                                                                                  ttATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTg                                                                  <  1:2425509/150‑1 (MQ=255)
                                                                                      gACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTAc                                                              >  2:2076314/1‑150 (MQ=255)
                                                                                      gACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTAc                                                              >  1:601565/1‑150 (MQ=255)
                                                                                         tATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCAAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGaa                                                           >  1:762582/1‑150 (MQ=255)
                                                                                          aTCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGaaa                                                          <  2:412440/150‑1 (MQ=255)
                                                                                           tCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGCCCGGCGCCGGCCAATAGAGAACGTCTAGCTATGGTAGTCTGTCGCCTGAGCGGTAACGGTTGGTACGAAAc                                                         <  1:646164/150‑1 (MQ=255)
                                                                                             aaTGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCg                                                       >  2:1311827/1‑150 (MQ=255)
                                                                                               tgtgCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGAc                                                     <  2:762582/150‑1 (MQ=255)
                                                                                               tgtgCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGAc                                                     <  2:601565/150‑1 (MQ=255)
                                                                                                  gCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGATGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTc                                                  <  1:1239781/150‑1 (MQ=255)
                                                                                                    gTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCaa                                                <  1:1151762/150‑1 (MQ=255)
                                                                                                       cTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGAc                                             >  2:415891/1‑150 (MQ=255)
                                                                                                        tGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGAcc                                            >  1:673711/1‑150 (MQ=255)
                                                                                                          cGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCgg                                          <  2:78635/150‑1 (MQ=255)
                                                                                                           gAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTACCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCggg                                         <  2:1800267/150‑1 (MQ=255)
                                                                                                            aGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCgggg                                        <  2:2098276/150‑1 (MQ=255)
                                                                                                               aCCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGa                                     <  1:1767734/150‑1 (MQ=255)
                                                                                                                ccGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCTTGTGCGGTAACGGTTgg                                                                 <  2:1920444/121‑1 (MQ=255)
                                                                                                                ccGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTgg                                                                 >  1:1920444/1‑121 (MQ=255)
                                                                                                                ccGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGaa                                    <  1:585133/150‑1 (MQ=255)
                                                                                                                  gCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGaaaa                                  >  2:2319744/1‑150 (MQ=255)
                                                                                                                   cATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTACTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGaaaaa                                 <  2:1300965/150‑1 (MQ=255)
                                                                                                                       gAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAAt                             <  1:287985/150‑1 (MQ=255)
                                                                                                                         aGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCt                           <  2:1611051/150‑1 (MQ=255)
                                                                                                                         aGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCt                           <  1:465283/150‑1 (MQ=255)
                                                                                                                          gTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTg                          <  2:673711/150‑1 (MQ=255)
                                                                                                                           tgctgcAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTgg                         <  1:2319744/150‑1 (MQ=255)
                                                                                                                                cAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAgctg                                                                                                        >  2:81574/1‑66 (MQ=255)
                                                                                                                                cAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAgctg                                                                                                        <  1:81574/66‑1 (MQ=255)
                                                                                                                                 aatcgaTCAGCGGCGCGGCGCTGATCTATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAggcg                   >  1:1359283/5‑150 (MQ=255)
                                                                                                                                   aagaTCAGCGCCGCGCCGCTGATCGATTGCTCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAggcggc                 <  1:1265686/150‑1 (MQ=255)
                                                                                                                                    agaTCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTTGGAGCTAGGGGTGCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAGgcggcg                <  2:191376/150‑1 (MQ=255)
                                                                                                                                     gaTCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAGGCGGCGa               >  2:1422018/1‑150 (MQ=255)
                                                                                                                                      aTCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAGGCGGCGAc              <  1:216127/150‑1 (MQ=255)
                                                                                                                                       tCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAggcggc                 <  2:611037/146‑1 (MQ=255)
                                                                                                                                       tCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAggcggc                 >  1:611037/1‑146 (MQ=255)
                                                                                                                                          gcgccgcgccGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTGGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAGGCGGCGACAAgg          >  2:2038023/1‑150 (MQ=255)
                                                                                                                                            gccgcgccgcTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAGGCGGCGACAAGGcc        <  2:337994/150‑1 (MQ=255)
                                                                                                                                                 gccgcTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTgg                         <  2:1056899/128‑1 (MQ=255)
                                                                                                                                                 gccgcTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTgg                         >  1:1056899/1‑128 (MQ=255)
                                                                                                                                                  ccgcTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAGGCGGCGACAAGGCCGACTGa  >  1:2021373/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TGCTTGCTAACCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAGGCGGCGACAAGGCCGACTGA  >  NC_002947/1389811‑1390106

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