Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 550,552 C→G 7.3% intergenic (+30/‑112) tufB → / → rpsJ elongation factor Tu‑B/30S ribosomal protein S10

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947550,5520CG 7.3% 248.4 / 10.9 82intergenic (+30/‑112)tufB/rpsJelongation factor Tu‑B/30S ribosomal protein S10
Reads supporting (aligned to +/- strand):  ref base C (37/39);  new base G (3/3);  total (40/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.28e-01

AACATTCAGATGACTGTTACCCTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTATCAGGTTGAAGGCTTT  >  NC_002947/550403‑550701
                                                                                                                                                     |                                                                                                                                                       
aacaTTCAGATGACTGTTACCCTGATCAAGACCAACGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTc                                                                                                                                                         >  1:1148326/1‑150 (MQ=255)
  caTTCAGATGACTGTTACCCTGATCAAGACCATCGCAATGGAAGACGGTCTGCGATTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAg                                                                                                                                                       >  2:572117/1‑150 (MQ=255)
   aTTCAGATGACTGTTACCCTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGc                                                                                                                                                      >  1:377645/1‑150 (MQ=255)
     tCAGATGACTGTTACCCTGCAAACGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCgg                                                                                                                                                    <  2:540448/150‑1 (MQ=255)
      cAGATGACTGTTACCCTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTGAGCggg                                                                                                                                                   <  1:963826/150‑1 (MQ=255)
       aGATGACTGTTACCCTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTGAGCgggg                                                                                                                                                  <  1:1393732/150‑1 (MQ=255)
         aTGACTGTTACCCTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGt                                                                                                                                               >  2:115324/1‑150 (MQ=255)
         aTGACTGTTACCCTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGt                                                                                                                                               <  2:444370/150‑1 (MQ=255)
           gACTGTTACCCTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGttt                                                                                                                                             <  1:1900934/150‑1 (MQ=255)
               gTTACCCTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGttttttt                                                                                                                                         <  1:867753/150‑1 (MQ=255)
                ttACCCTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGtttttttt                                                                                                                                        >  2:1847982/1‑150 (MQ=255)
                  aCCCTGATCAAGACCATCGCAATGGAAGCCGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGATTTTTTTAt                                                                                                                                      >  2:748847/1‑150 (MQ=255)
                   cccTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTAtt                                                                                                                                     >  2:1097744/1‑150 (MQ=255)
                    ccTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTg                                                                                                                                    >  1:1361098/1‑150 (MQ=255)
                      tGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTCATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTggg                                                                                                                                  <  2:140400/150‑1 (MQ=255)
                       gATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGt                                                                                                                                 <  2:2091224/150‑1 (MQ=255)
                           aaGACCCTCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTTATTGGGTTGa                                                                                                                              <  2:59433/150‑1 (MQ=255)
                           aaGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGac                                                                                                                             >  1:772376/1‑150 (MQ=255)
                             gACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGacac                                                                                                                           >  2:394829/1‑150 (MQ=255)
                                  tCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAAtt                                                                                                                      >  1:265439/1‑150 (MQ=255)
                                   cGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTg                                                                                                                     <  1:1835633/150‑1 (MQ=255)
                                      aaTGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCgggg                                                                                                                                                  >  1:369234/1‑119 (MQ=255)
                                      aaTGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCgggg                                                                                                                                                  <  2:369234/119‑1 (MQ=255)
                                      aaTGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCggg                                                                                                                                                   >  1:1305913/1‑118 (MQ=255)
                                      aaTGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCggg                                                                                                                                                   <  2:1305913/118‑1 (MQ=255)
                                      aaTGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTgggg                                                                                                                  <  1:352603/150‑1 (MQ=255)
                                        tGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGc                                                                                                                <  1:1982855/150‑1 (MQ=255)
                                          gAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGt                                                                                                              >  1:936997/1‑150 (MQ=255)
                                           aaGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTc                                                                                                             <  1:882491/150‑1 (MQ=255)
                                             gACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCta                                                                                                           <  1:394829/150‑1 (MQ=255)
                                              aCGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCtat                                                                                                          <  2:1325783/150‑1 (MQ=255)
                                              aCGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCtat                                                                                                          <  1:621580/150‑1 (MQ=255)
                                              aCGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCtat                                                                                                          <  1:804298/150‑1 (MQ=255)
                                                       gcTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACg                                                                                                 >  2:316455/1‑150 (MQ=255)
                                                        cTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGc                                                                                                <  2:1148326/150‑1 (MQ=255)
                                                         ttCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGcc                                                                                               >  1:1706624/1‑150 (MQ=255)
                                                         ttCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGcc                                                                                               <  1:2255829/150‑1 (MQ=255)
                                                               aTCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCttt                                                                                         >  1:974897/1‑150 (MQ=255)
                                                                    tGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTggg                                                                                                                   >  1:50365/1‑119 (MQ=255)
                                                                    tGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTggg                                                                                                                   <  2:50365/119‑1 (MQ=255)
                                                                     gAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCgg                                                                                   >  2:396518/1‑150 (MQ=255)
                                                                      aaGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCggg                                                                                  <  1:1097744/150‑1 (MQ=255)
                                                                           ggTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAg                                                                             <  1:106648/150‑1 (MQ=255)
                                                                              cGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGc                                                                         <  2:974897/150‑1 (MQ=255)
                                                                               gTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGaa                                                                                                      >  1:1165036/1‑121 (MQ=255)
                                                                               gTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGaa                                                                                                      <  2:1165036/121‑1 (MQ=255)
                                                                                tACCGTCGGCGCCGGAGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGATTTTTTTATTGGGTTGACACTAATTGGGGGCGTATATAGAATCACGCCTCCTTTTAGCGGGCGTGATGCGt                                                                       >  2:394137/1‑150 (MQ=255)
                                                                                    gTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGt                                                                   <  1:572117/150‑1 (MQ=255)
                                                                                          gcCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGTGGGCGTAGTGCGTCCGTTGGGaa                                                             <  1:316455/150‑1 (MQ=255)
                                                                                            cGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTTTAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTg                                                                 >  2:1119823/1‑144 (MQ=255)
                                                                                            cGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTg                                                                 <  1:1119823/144‑1 (MQ=255)
                                                                                                   gtgGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGa                                                    <  2:772376/150‑1 (MQ=255)
                                                                                                   gtgACCAAAACCATCGAGTCTTCACTCAAGGGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGa                                                    <  2:1667866/150‑1 (MQ=255)
                                                                                                    tgGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAg                                                   >  1:948158/1‑150 (MQ=255)
                                                                                                      gCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTc                                                 >  2:1759252/1‑150 (MQ=255)
                                                                                                       ccAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCt                                                <  1:115324/150‑1 (MQ=255)
                                                                                                        cAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTg                                               >  1:2326637/1‑150 (MQ=255)
                                                                                                              catcGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATcc                                         >  2:212718/1‑150 (MQ=255)
                                                                                                                tcGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTg                                                      >  1:14647/1‑135 (MQ=255)
                                                                                                                tcGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTg                                                      <  2:14647/135‑1 (MQ=255)
                                                                                                                tcGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCaa                                       >  1:879745/1‑150 (MQ=255)
                                                                                                                tcGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCaa                                       >  2:2115007/1‑150 (MQ=255)
                                                                                                                tcGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGAGTCTATAGACTCACGCCTCCTTTTAGCGTGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCaa                                       <  2:265439/150‑1 (MQ=255)
                                                                                                                     tAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCaa                                  <  2:1361098/150‑1 (MQ=255)
                                                                                                                       aTCACTCGATCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCaaaa                                >  2:405036/1‑150 (MQ=255)
                                                                                                                       aTCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAAt                                      <  1:1879685/144‑1 (MQ=255)
                                                                                                                       aTCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAAt                                      >  2:1879685/1‑144 (MQ=255)
                                                                                                                       aTCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCaaaa                                >  2:2292078/1‑150 (MQ=255)
                                                                                                                          aCTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCt                                                                                                            <  1:1395004/72‑1 (MQ=255)
                                                                                                                          aCTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCt                                                                                                            >  2:1395004/1‑72 (MQ=255)
                                                                                                                                  gATTGAAAAACGCCACGGCCCGCGGCGGATTTCTTCACTGGCTGGACACTGAGTGGGGGCGTCTATATAATCACGCTTCCTTTTAGCGGGAGTAGTGCGTCCGCTGCGAACAGCCTGGAGTCTGAAATCCACTGCAAAATCAGGAAATcc                     >  2:1780051/1‑150 (MQ=255)
                                                                                                                                  gATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATcc                     >  2:1087488/1‑150 (MQ=255)
                                                                                                                                  gATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATcc                     <  1:2115007/150‑1 (MQ=255)
                                                                                                                                  aataaaaaaaaCCCCCGCTGAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATcc                     >  2:1942905/6‑150 (MQ=255)
                                                                                                                                         aaaaCCCCCGCTGAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTACAGATCGGAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTATCAg              >  1:145060/1‑150 (MQ=255)
                                                                                                                                         aaaaCCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAGGGGAAGGCAAAATCAGCAAATCCGTATCAg              >  2:2046043/1‑150 (MQ=255)
                                                                                                                                         aaaaCCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTATCAg              >  2:2029790/1‑150 (MQ=255)
                                                                                                                                          aaaCCCCCGCTCGGCGGGGGTCTTATTAATTGGGTGGACACTAAGTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTAGCAgg             <  2:879745/150‑1 (MQ=255)
                                                                                                                                            aCCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTATCAGGtt           <  2:1706624/150‑1 (MQ=255)
                                                                                                                                              ccccGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAATCTGAAATCCAATGCAAAATCAGCAAATCCGTATCAGGTTGa         >  2:1644193/1‑150 (MQ=255)
                                                                                                                                                ccGCTGAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTc                                                 >  1:755013/1‑108 (MQ=255)
                                                                                                                                                ccGCTGAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTc                                                 <  2:755013/108‑1 (MQ=255)
                                                                                                                                                ccGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTATCAGGTTGAAg       <  1:2029790/150‑1 (MQ=255)
                                                                                                                                                     cAGCGGGGGTTTTTTTTATTGGGTTGCCCCTAATTGGGGGCGTCTATAGAATCACGCCTCCTTAAGCCCCAGGAGTTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTATCAGGTTGAAGGCttt  <  1:405036/150‑1 (MQ=255)
                                                                                                                                                     cAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTATCAGGTTGAAGGCttt  <  2:948158/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
AACATTCAGATGACTGTTACCCTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTATCAGGTTGAAGGCTTT  >  NC_002947/550403‑550701

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