Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 550,552 C→G 11.0% 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,5520CG11.0% 265.0 / 22.1 91intergenic (+30/‑112)tufB/rpsJelongation factor Tu‑B/30S ribosomal protein S10
Reads supporting (aligned to +/- strand):  ref base C (38/43);  new base G (6/4);  total (44/47)
Fisher's exact test for biased strand distribution p-value = 5.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.37e-01

AACATTCAGATGACTGTTACCCTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTATCAGGTTGAAGGC  >  NC_002947/550403‑550698
                                                                                                                                                     |                                                                                                                                                  
aacaTTCAGATGACTGTTACCCTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTc                                                                                                                                                    >  1:1606157/1‑150 (MQ=255)
   aTTCAGATGACTGTTACCCTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGc                                                                                                                                                 >  1:188818/1‑150 (MQ=255)
    ttCAGATGACTGTTACCCTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCg                                                                                                                                                >  1:601847/1‑150 (MQ=255)
      cAGATGACTGTTACCCTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCggg                                                                                                                                              <  1:1638048/150‑1 (MQ=255)
       aGATGACTGTTACCCTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCgggg                                                                                                                                             <  2:946260/150‑1 (MQ=255)
        gATGACTGTTACCCTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCATAGCTATCCGTGAAGGCTGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCCATTGAAAAAACCCCCGCTCAGCggggg                                                                                                                                            >  2:1010588/1‑150 (MQ=255)
         aTGACTGTTACCCTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGt                                                                                                                                           >  1:886312/1‑150 (MQ=255)
           gACTGTTACCCTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGttt                                                                                                                                         <  2:1740863/150‑1 (MQ=255)
              tGTTACCCTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGtttttt                                                                                                                                      <  1:56290/150‑1 (MQ=255)
              tGTTACCCTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGtttttt                                                                                                                                      >  2:1529157/1‑150 (MQ=255)
                ttACCCTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGtttttttt                                                                                                                                    >  2:23572/1‑150 (MQ=255)
                ttACCCTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGtttttttt                                                                                                                                    <  2:373202/150‑1 (MQ=255)
                    ccTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTg                                                                                                                                <  2:1298962/150‑1 (MQ=255)
                     cTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTgg                                                                                                                               >  1:599329/1‑150 (MQ=255)
                     cTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTgg                                                                                                                               <  2:599329/150‑1 (MQ=255)
                      tGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTggg                                                                                                                              >  2:1440525/1‑150 (MQ=255)
                       gATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGt                                                                                                                             <  1:1440525/150‑1 (MQ=255)
                        aTCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGtt                                                                                                                            <  2:886312/150‑1 (MQ=255)
                           aaGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTGAGCgg                                                                                                                                               >  2:776493/1‑128 (MQ=255)
                           aaGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTGAGCgg                                                                                                                                               <  1:776493/128‑1 (MQ=255)
                             gACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGacac                                                                                                                       >  1:262872/1‑150 (MQ=255)
                                cATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTaa                                                                                                                    >  2:1362862/1‑150 (MQ=255)
                                cATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTaa                                                                                                                    <  1:387180/150‑1 (MQ=255)
                                cATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTaa                                                                                                                    >  1:893281/1‑150 (MQ=255)
                                    gCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTgg                                                                                                                <  2:1010109/150‑1 (MQ=255)
                                    gCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTATTGGGTTGACACTAATTggg                                                                                                               >  1:1378584/1‑150 (MQ=255)
                                          gAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGt                                                                                                          <  2:1606157/150‑1 (MQ=255)
                                          gAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGt                                                                                                          >  2:143272/1‑150 (MQ=255)
                                             gACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGaca                                                                                                                        <  1:997277/133‑1 (MQ=255)
                                             gACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGaca                                                                                                                        >  2:997277/1‑133 (MQ=255)
                                               cGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTAGCGTCGGCGCCGGCATTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCtata                                                                                                     >  1:626780/1‑150 (MQ=255)
                                                   cTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATc                                                                                                <  2:1378584/150‑1 (MQ=255)
                                                            gCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGcctcc                                                                                        <  2:601847/150‑1 (MQ=255)
                                                            gCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGcctcc                                                                                        <  1:1653264/150‑1 (MQ=255)
                                                                       aGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGc                                                                             <  1:143272/150‑1 (MQ=255)
                                                                          cGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTa                                                                          >  1:776206/1‑150 (MQ=255)
                                                                          cGGTCGTACCGTCGGCGCCGGCGTTGCGGCCAAAATCATCGAGTAATCACTCTACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTa                                                                          <  2:892318/150‑1 (MQ=255)
                                                                            gTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGt                                                                        <  1:23572/150‑1 (MQ=255)
                                                                                   cGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGt                                                                        <  2:730551/143‑1 (MQ=255)
                                                                                   cGTCGGCGCCGGCGTTGTGGCCAAAAGCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGt                                                                        >  1:730551/1‑143 (MQ=255)
                                                                                     tCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCggg                                                                                                                                              <  1:1369867/71‑1 (MQ=255)
                                                                                     tCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCggg                                                                                                                                              >  2:1369867/1‑71 (MQ=255)
                                                                                      cGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTg                                                              <  1:1362862/150‑1 (MQ=255)
                                                                                       ggCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTgg                                                             <  1:100881/150‑1 (MQ=255)
                                                                                             ggCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAg                                                       >  1:913004/1‑150 (MQ=255)
                                                                                                 ttGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTg                                                   <  2:188818/150‑1 (MQ=255)
                                                                                                     ggCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGt                                               <  1:1904197/150‑1 (MQ=255)
                                                                                                     ggCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCGCCTTTTAGCGGGCGTGGTGCGTCCGTTGGGAACAGCCTGGAGt                                               <  2:262872/150‑1 (MQ=255)
                                                                                                        cAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGATACTAATTGGGGGCCTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCTTCCGTTGGGAACAGCCTGGAGTCTg                                            >  1:820442/1‑150 (MQ=255)
                                                                                                            atcatcGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGGCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGACTGGAGTATGaaaa                                        >  2:372344/1‑149 (MQ=255)
                                                                                                               gtcGAGTGATTACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTCTATTGGGGTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACCGCCTGGAGTCTGAAATCCa                                     >  2:255804/2‑150 (MQ=255)
                                                                                                                tcGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCaa                                    >  1:1013696/1‑150 (MQ=255)
                                                                                                                tcGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCaa                                    <  2:104069/150‑1 (MQ=255)
                                                                                                                tcGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCaa                                    <  2:776206/150‑1 (MQ=255)
                                                                                                                tcGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCaa                                    >  2:343542/1‑150 (MQ=255)
                                                                                                                tcGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCaa                                    <  1:1426886/150‑1 (MQ=255)
                                                                                                                tcGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCaa                                    <  1:1426679/150‑1 (MQ=255)
                                                                                                                     tAATCACTCGACCGATTGAAAAAGCCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCaa                               <  2:141325/150‑1 (MQ=255)
                                                                                                                       aTCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTggg                                                                                                               <  1:1039580/68‑1 (MQ=255)
                                                                                                                       aTCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTggg                                                                                                               >  2:1039580/1‑68 (MQ=255)
                                                                                                                           cgcGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGATGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTGTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGACTGAAATCCAATGCAAAACCAg                         >  2:1763376/3‑150 (MQ=255)
                                                                                                                            tCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTa                                                                                   >  1:1050461/1‑91 (MQ=255)
                                                                                                                            tCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTa                                                                                   <  2:1050461/91‑1 (MQ=255)
                                                                                                                            tCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGa                                                 <  2:579115/125‑1 (MQ=255)
                                                                                                                            tCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGa                                                 >  1:579115/1‑125 (MQ=255)
                                                                                                                            tCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGaaa                                         >  1:1123273/1‑133 (MQ=255)
                                                                                                                            tCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGaaa                                         <  2:1123273/133‑1 (MQ=255)
                                                                                                                              gACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCaa                      >  1:1104083/1‑150 (MQ=255)
                                                                                                                               aCCGATTGAAAAAGCCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCaaa                     <  1:1010588/150‑1 (MQ=255)
                                                                                                                                  gATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATcc                  >  2:1593981/1‑150 (MQ=255)
                                                                                                                                       aaaaaaCCCCCGCTGAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTATc             >  2:632158/1‑150 (MQ=255)
                                                                                                                                       aaaaaaCCCCCGCTGAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTATc             >  2:621242/1‑150 (MQ=255)
                                                                                                                                       aaaaaaCCCCCGATCCGCCGGGGGTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTATc             >  2:643984/1‑150 (MQ=255)
                                                                                                                                        aaaaaCCCCCGCTCAGCGGGGGATTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTATCa            >  1:1825172/1‑150 (MQ=255)
                                                                                                                                         aaaaCCCCCGCTGAGCGGGGGTTTTTTTTATTGGGGTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTATCAg           >  1:1086585/1‑150 (MQ=255)
                                                                                                                                         aaaaCCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAGCAGCCTGGAGTCTGCAATCGGATGCAAAATCAGCTCATCCGCATccg           >  2:99488/1‑148 (MQ=255)
                                                                                                                                          aaaCCCCCGCTGAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGc                                                                                            >  1:697741/1‑68 (MQ=255)
                                                                                                                                          aaaCCCCCGCTGAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGc                                                                                            <  2:697741/68‑1 (MQ=255)
                                                                                                                                          aaaCCCCCGCTGAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGa                                                 <  2:589163/111‑1 (MQ=255)
                                                                                                                                          aaaCCCCCGCTGAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGa                                                 >  1:589163/1‑111 (MQ=255)
                                                                                                                                          aaaCCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGtt                                                               >  1:1919660/1‑97 (MQ=255)
                                                                                                                                          aaaCCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGtt                                                               >  1:1919897/1‑97 (MQ=255)
                                                                                                                                          aaaCCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGtt                                                               <  2:1919897/97‑1 (MQ=255)
                                                                                                                                          aaaCCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGtt                                                               <  2:1919660/97‑1 (MQ=255)
                                                                                                                                          aaaCCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTATCAgg          <  2:1862694/150‑1 (MQ=255)
                                                                                                                                            aCCCCCGCTGAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTATCAGGtt        <  1:632158/150‑1 (MQ=255)
                                                                                                                                               cccGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTg                                                   >  2:436046/1‑104 (MQ=255)
                                                                                                                                               cccGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTg                                                   <  1:436046/104‑1 (MQ=255)
                                                                                                                                                ccGCTCAGCGGGGGTTTTTTTTATTTGGTTTATACTATTTTTTGGTTGCTTTCGAATTACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTATCAGGTTGAAg    >  2:1809643/1‑150 (MQ=255)
                                                                                                                                                ccGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTATCAGGTTGAAg    <  1:1529157/150‑1 (MQ=255)
                                                                                                                                                 cGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTATCAGGTTGAAgg   <  1:1593981/150‑1 (MQ=255)
                                                                                                                                                  gCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCTATGCAAAATCAGCAAATCCGTATCAGGTTGTAggt  <  2:1013696/150‑2 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
AACATTCAGATGACTGTTACCCTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTATCAGGTTGAAGGC  >  NC_002947/550403‑550698

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: