Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,854,053 A→G 6.5% intergenic (+23/‑248) rlmI → / → ilvD ribosomal RNA large subunit methyltransferase I/dihydroxy‑acid dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,854,0530AG 6.5% 1262.4 / 47.6 402intergenic (+23/‑248)rlmI/ilvDribosomal RNA large subunit methyltransferase I/dihydroxy‑acid dehydratase
Reads supporting (aligned to +/- strand):  ref base A (224/150);  new base G (13/13);  total (237/165)
Fisher's exact test for biased strand distribution p-value = 4.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.21e-01

Only 100 of 405 total aligned reads displayed.

CCACCTGGACCGCAACATGCAGCTGCTCGAGCGAGGCGGTCAGGGCCCGGATCACCCGGTACACCCGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTATGAGCTCTGGTCAAGCACGCGGCGATCCCGCGCGGTCT  >  NC_002947/5853904‑5854202
                                                                                                                                                     |                                                                                                                                                     
  ggatGGACCGCAACATGCAGCTGCTCGAGCGAGGCGGTCAGGGCCCGGATCACCCGGTACACCCGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGc                                                                                                                                                     <  1:4590039/147‑1 (MQ=255)
   ccTGGACCGCAACATGCAGCTGCTCGAGCGAGGCGGTCAGGGCCCGGATCACCCGGTACACCCGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCg                                                                                                                                                    <  2:2161983/150‑1 (MQ=255)
     tGGACCGAAACATGCAGCTGCTCGAGCGAGGCGGTCAGGGCCCGGATCACCCGGTACACCCGGCCATCGACGAAACCCGCTACATCAAAAGCATAACCTGCCGGTAACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGc                                                                                                                                                  >  1:6094288/1‑150 (MQ=255)
      ggACCGCAACATGCAGCTGCTCGAGCGAGGCGGTCAGGGCCCGGATCACCCGGTACACCCGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGcc                                                                                                                                                 <  1:3224288/150‑1 (MQ=255)
        aCCGCAACATGCAGCTGCTCGAGCGAGGCGGTCAGGGCCCGGATCACCCGGTACACCCGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCTGGGCGGCCCt                                                                                                                                               <  1:355695/150‑1 (MQ=255)
        aCCGCAACATGCAGCTGCTCGAGCGAGGCGGTCAGGGCCCGGATCACCCGGTACACCCGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCt                                                                                                                                               >  1:4241806/1‑150 (MQ=255)
          cGCAACATGCAGCTGCTCGAGCGAGGCGGTCAGGGCCCGGATCACCCGGTACACCCGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCttt                                                                                                                                             >  2:4832124/1‑150 (MQ=255)
           gCAACATGCAGCTGCTCGAGCGAGGCGGTCAGGGCCCGGATCACCCGGTACACCCGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCTGGGCGGCCCtttt                                                                                                                                            <  1:3713995/150‑2 (MQ=255)
            cAACATGCAGCTGCTCGAGCGAGGCGGTCAGGGCCCGGATCACCCGGTACACCCGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCTGGGCGGCCCttttt                                                                                                                                           <  2:2978802/150‑3 (MQ=255)
              aCATGCAGCTGCTCGAGCGAGGCGGTCAGGGCCCGGATCACCCGGTACACCCGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCTGGGCGGCCCtttttgg                                                                                                                                         <  2:3446437/150‑5 (MQ=255)
                aTGCAGCTGCTCGAGCGAGGCGGTCAGGGCCCGGATCACCCGGTACACCCGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGcc                                                                                                                                       <  2:3755262/150‑1 (MQ=255)
                aTGCAGCTGCTCGAGCGAGGCGGTCAGGGCCCGGATCACCCGGTACACCCGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGcc                                                                                                                                       <  2:2345766/150‑1 (MQ=255)
                aTGCAGCTGCTCGAGCGAGGCGGTCAGGGCCCGGATCACCCGGTACACCCGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGcc                                                                                                                                       <  2:709987/150‑1 (MQ=255)
                    aGCTGCTCGAGCGAGGCGGTCAGGGCCCGGATCACCCGGTACACCCGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCaa                                                                                                                                   >  2:4392303/1‑150 (MQ=255)
                    aGCTGCTCGAGCGAGGCGGTCAGGGCCCGGATCACCCGGTACACCCGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCaa                                                                                                                                   <  2:5387025/150‑1 (MQ=255)
                         ctCGAGCGAGGCGGTCAGGGCCCGGATCACCCGGTACACCCGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATcccc                                                                                                                              >  1:5786538/1‑150 (MQ=255)
                         ctCGAGCGAGGCGGTCAGGGCCCGGATCACCCGGTACACCCGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATcccc                                                                                                                              >  1:1017166/1‑150 (MQ=255)
                         ctCGAGCGAGGCGGTCAGGGCCCGGATCACCCGGTACACCCGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATcccc                                                                                                                              >  1:4366288/1‑150 (MQ=255)
                          tCGAGCGAGGCGGTCAGGGCCCGGATCACCCGGTACACCCGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATccccc                                                                                                                             >  2:423802/1‑150 (MQ=255)
                          tCGAGCGAGGCGGTCAGGGCCCGGATCACCCGGTACACCCGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATccccc                                                                                                                             >  1:6142693/1‑150 (MQ=255)
                            gagcgagGCGGTCAGGGCCCGGATCACCCGGTACACCCGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCg                                                                                                                           >  1:4668209/1‑150 (MQ=255)
                                gagGCGGTCAGGGCCCGGATCACCCGGTACACCCGGCCATCGCCGAAACCCGCTACATCAACAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCaa                                                                                                                       <  2:6845048/150‑1 (MQ=255)
                                  ggCGGTCAGGGCCCGGATCACCCGGTACACCCGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAAcc                                                                                                                     <  1:6785453/150‑1 (MQ=255)
                                     ggTCAGGGCCCGGATCACCCGGTACACCCGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATc                                                                                                                  <  1:5056112/150‑1 (MQ=255)
                                           ggCCCGGATCACCCGGTACACCCGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGTAATCCCCCCGCCAACCATCCACATc                                                                                                            >  1:1363099/1‑150 (MQ=255)
                                           ggCCCGGATCACCCGGTACACCCGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCAACCACATc                                                                                                            >  2:178064/1‑150 (MQ=255)
                                            gCCCGGATCACCCGGGACACCCGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCt                                                                                                           <  2:3110681/150‑1 (MQ=255)
                                                   tCACCCGGTACACCCGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTcc                                                                                                    >  1:2139637/1‑150 (MQ=255)
                                                        cGGTACACCCGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCAc                                                                                               <  2:3898590/150‑1 (MQ=255)
                                                         ggTACACCCGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCAcc                                                                                              >  1:2576448/1‑150 (MQ=255)
                                                           tACACCCGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACcgc                                                                                            >  1:5267209/1‑150 (MQ=255)
                                                              acCCGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGcc                                                                                         >  1:1450444/1‑150 (MQ=255)
                                                               cccGGCCATCGCCGAAACCCGCTACAGCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCa                                                                                        >  2:1889116/1‑150 (MQ=255)
                                                                ccGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCGTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAg                                                                                       >  2:5160470/1‑150 (MQ=255)
                                                                   gCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCgg                                                                                    <  1:1567926/150‑1 (MQ=255)
                                                                    ccATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGgt                                                                                   >  2:6558114/1‑150 (MQ=255)
                                                                    ccATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGgt                                                                                   >  2:778067/1‑150 (MQ=255)
                                                                    ccATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGgt                                                                                   <  1:4416847/150‑1 (MQ=255)
                                                                     cATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGgtg                                                                                  <  2:669981/150‑1 (MQ=255)
                                                                      aTCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGgtgt                                                                                 >  1:2036444/1‑150 (MQ=255)
                                                                       tCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTa                                                                                <  2:356860/150‑1 (MQ=255)
                                                                          ccGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGaa                                                                             >  1:1745853/1‑150 (MQ=255)
                                                                            gAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATc                                                                           <  2:571172/150‑1 (MQ=255)
                                                                                cccGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGcc                                                                       >  1:1820816/1‑150 (MQ=255)
                                                                                 ccGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCt                                                                      >  2:4520847/1‑150 (MQ=255)
                                                                                 ccGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCt                                                                      >  2:2765928/1‑150 (MQ=255)
                                                                                 ccGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCt                                                                      >  2:433179/1‑150 (MQ=255)
                                                                                 ccGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCt                                                                      >  2:6233237/1‑150 (MQ=255)
                                                                                   gCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTAt                                                                    >  2:4637773/1‑150 (MQ=255)
                                                                                    cTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTAtt                                                                   >  2:6014899/1‑150 (MQ=255)
                                                                                        aTCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATc                                                               >  2:5922571/1‑150 (MQ=255)
                                                                                        aTCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATc                                                               >  2:2711036/1‑150 (MQ=255)
                                                                                         tCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGTGCCAGCGGTGTAGAATCGGCCTATTCATCg                                                              >  1:287563/1‑150 (MQ=255)
                                                                                         tCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCg                                                              >  2:6355292/1‑150 (MQ=255)
                                                                                         tCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCg                                                              >  1:410229/1‑150 (MQ=255)
                                                                                         tCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCg                                                              >  1:3848628/1‑150 (MQ=255)
                                                                                            aaaGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCa                                                           >  2:6452520/1‑150 (MQ=255)
                                                                                            aaaGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCa                                                           >  2:1220032/1‑150 (MQ=255)
                                                                                             aaGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAg                                                          >  1:739675/1‑150 (MQ=255)
                                                                                                 aTCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCAt                                                      >  2:2952583/1‑150 (MQ=255)
                                                                                                    aCCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATccc                                                   <  2:1724918/150‑1 (MQ=255)
                                                                                                         ccGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCcggcg                                              >  1:4805466/1‑150 (MQ=255)
                                                                                                          cGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCcggcgg                                             <  1:5369805/150‑1 (MQ=255)
                                                                                                               aCTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTa                                        >  1:6194756/1‑150 (MQ=255)
                                                                                                                cTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACCTCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTAt                                       <  1:3161742/150‑1 (MQ=255)
                                                                                                                cTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTAt                                       <  1:628957/150‑1 (MQ=255)
                                                                                                                 tGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTATg                                      >  2:5448700/1‑150 (MQ=255)
                                                                                                                 tGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTATg                                      >  2:2820874/1‑150 (MQ=255)
                                                                                                                 tGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTATg                                      >  2:119997/1‑150 (MQ=255)
                                                                                                                  gCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTATGa                                     <  1:6431544/150‑1 (MQ=255)
                                                                                                                  gCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTATGa                                     <  1:6247309/150‑1 (MQ=255)
                                                                                                                     aaaCAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTATGAGct                                  >  1:3068271/1‑150 (MQ=255)
                                                                                                                       aCAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTCTGAGctct                                >  1:2283132/1‑150 (MQ=255)
                                                                                                                       aCAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTATGAGctct                                >  2:5651797/1‑150 (MQ=255)
                                                                                                                        cAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCCCATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTATGAGCTCTg                               <  1:5077954/150‑1 (MQ=255)
                                                                                                                        cAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTATGAGCTCTg                               <  2:1517389/150‑1 (MQ=255)
                                                                                                                        cAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTATGAGCTCTg                               <  2:6094157/150‑1 (MQ=255)
                                                                                                                            tGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTATGAGCTCTGGTCa                           >  1:2726956/1‑150 (MQ=255)
                                                                                                                             gATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCgcgc                                                                                          >  2:495449/1‑86 (MQ=255)
                                                                                                                             gATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTATGAGCTCTGGTCaa                          <  1:2668651/150‑1 (MQ=255)
                                                                                                                             gATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTATGAGCTCTGGTCaa                          <  2:337534/150‑1 (MQ=255)
                                                                                                                                 tCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTATGAGCTCTGGTCAAGCAc                      >  2:6668619/1‑150 (MQ=255)
                                                                                                                                   tGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTATGAGCTCTGGTCAAGCAcgc                    >  1:5427396/1‑150 (MQ=255)
                                                                                                                                   ggcaggaaaGGGCCGCTGGGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTATGAGCTCTGGTCAAGCAcgc                    >  2:5276307/7‑150 (MQ=255)
                                                                                                                                    gcaggaaaGGGCCGCTGGGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTATGAGCTCTGGTCAAGCAcgcg                   >  1:1312623/6‑150 (MQ=255)
                                                                                                                                    gCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTATGAGCTCTGGTCAAGCAcgcg                   >  1:1069085/1‑150 (MQ=255)
                                                                                                                                    gCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTATGAGCTCTGGTCAAGCAcgcg                   >  1:4511553/1‑150 (MQ=255)
                                                                                                                                      cAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTATGAGCTCTGGTCAAGCACgcggc                 >  1:6643323/1‑150 (MQ=255)
                                                                                                                                       aaaaaGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTATGAGCTCTGGTCAAGCACgcggcg                >  2:1005731/1‑150 (MQ=255)
                                                                                                                                       aaaaaGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTATGAGCTCTGGTCAAGCACgcggcg                >  1:696272/1‑150 (MQ=255)
                                                                                                                                       aaaaaGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTATGAGCTCTGGTCAAGCACgcggcg                >  1:4937774/1‑150 (MQ=255)
                                                                                                                                       aaaaaGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGGTTATGAGCTCTGGTCAAGCACgcggcg                >  2:2049455/1‑150 (MQ=255)
                                                                                                                                        gaaaGGGCCGCTGGGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATc                                                                           <  1:5675379/89‑1 (MQ=39)
                                                                                                                                        gaaaGGGCCGCTGGGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTCTGAGCTCTGGTCAAGCACGCGGCGa               >  1:6047795/2‑150 (MQ=255)
                                                                                                                                        aaaaGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTATGAGCTCTgg                              >  1:3427256/1‑135 (MQ=255)
                                                                                                                                        aaaaGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTATGAGCTCTGGTCAAGCACGCGGCGa               <  1:5970628/150‑1 (MQ=255)
                                                                                                                                           aGGGCCGCTGGGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCAACGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTATGAGCTCTGGTCAAGCACGCGGCGATcc            >  1:2011997/1‑150 (MQ=255)
                                                                                                                                             ggCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTCTGAGCTCTGGTCAAGCACGCGGCGATCCcg          >  2:5002536/1‑150 (MQ=255)
                                                                                                                                               ccgccCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTATGAGCTCTGGTCAAGCACGCGGCGATCCcgcg        <  2:242686/150‑1 (MQ=255)
                                                                                                                                                cgccCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTATGAGCTCTGGTCAAGCACGCGGCGATCCcgcgc       <  1:2478036/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCACCTGGACCGCAACATGCAGCTGCTCGAGCGAGGCGGTCAGGGCCCGGATCACCCGGTACACCCGGCCATCGCCGAAACCCGCTACATCAAAAGCATCACCTGCCGGTTACTGCCAAACAGCTGATATCTGCCAAAAAGGGCCGCCCAGCGGCCCTTTCCTGCCGCAATCCCCCCGCCAACCATCCACATCTGCATTCCCTCACCGCGCCAGCGGTGTAGAATCGGCCTATTCATCGCCAGTCATCCCCGGCGGGTTTATGAGCTCTGGTCAAGCACGCGGCGATCCCGCGCGGTCT  >  NC_002947/5853904‑5854202

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