Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 843,338 C→G 6.1% A715P (GCG→CCG)  ybiO ← mechanosensitive channel protein, intermediate conductance

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913843,3380CG 6.1% 390.7 / 14.0 114A715P (GCG→CCG) ybiOmechanosensitive channel protein, intermediate conductance
Reads supporting (aligned to +/- strand):  ref base C (58/49);  new base G (3/4);  total (61/53)
Fisher's exact test for biased strand distribution p-value = 7.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

Only 100 of 121 total aligned reads displayed.

TCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTGCTTAAGCCGACAATCCCGGCAAAATTCGGTTC  >  NC_000913/843192‑843485
                                                                                                                                                  |                                                                                                                                                   
tCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGgcgcgc                                                                                                                                                   >  3:729153/1‑149 (MQ=255)
tCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGgcgcgc                                                                                                                                                   >  3:161708/1‑149 (MQ=255)
tCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGgcgc                                                                                                                                                     <  3:725526/147‑1 (MQ=255)
tCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGgcgc                                                                                                                                                     >  4:725526/1‑147 (MQ=255)
   cTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGaa                                                                                                                                                <  3:274077/149‑1 (MQ=255)
    tttATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAAc                                                                                                                                               <  4:638177/149‑1 (MQ=255)
       aTTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGGTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAAcgcc                                                                                                                                            >  1:641895/1‑149 (MQ=255)
        ttGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAAcgccc                                                                                                                                           >  3:649639/1‑149 (MQ=255)
            gtGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCGCCCGGCGCTAACGGTTCAGCCGGGGTCGGGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCa                                                                                                                                       >  2:372585/2‑149 (MQ=255)
            cTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCCCCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCa                                                                                                                                       >  2:284052/1‑149 (MQ=255)
            cTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAAcgccc                                                                                                                                           >  2:772286/1‑145 (MQ=255)
            cTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAAcgccc                                                                                                                                           <  1:772286/145‑1 (MQ=255)
             tGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAg                                                                                                                                      >  3:455209/1‑149 (MQ=255)
              ggACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGGGCGTTCGCGCGCCAgt                                                                                                                                     <  2:655166/149‑2 (MQ=255)
                  cTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGTGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGa                                                                                                                                 >  2:222968/1‑149 (MQ=255)
                  cTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGa                                                                                                                                 >  4:883790/1‑149 (MQ=255)
                  cTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGa                                                                                                                                 >  1:702679/1‑149 (MQ=255)
                   tttATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGaa                                                                                                                                >  4:725853/1‑149 (MQ=255)
                    ttATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGaaa                                                                                                                               <  4:811169/149‑1 (MQ=255)
                    ttATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGaaa                                                                                                                               <  3:883790/149‑1 (MQ=255)
                       tGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGt                                                                                                                            >  2:202948/1‑149 (MQ=255)
                        gAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGGGCGTTCGCGCGCCaa                                                                                                                                      <  1:508310/138‑2 (MQ=255)
                        gAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGGGCGTTCGCGCGCCaa                                                                                                                                      >  2:508310/1‑137 (MQ=255)
                        gAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAAc                                                                                                                                               <  3:868887/129‑1 (MQ=255)
                        gAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAAc                                                                                                                                               >  4:868887/1‑129 (MQ=255)
                          acgacgAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGGGCGTTCGCGCGCCAgt                                                                                                                                     >  3:886474/1‑136 (MQ=255)
                           cgacgaGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGttttt                                                                                                                        <  4:604258/149‑1 (MQ=255)
                            gacgaGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGtttttt                                                                                                                       >  2:175442/1‑149 (MQ=255)
                            gacgaGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGtttttt                                                                                                                       <  2:867337/149‑1 (MQ=255)
                            gacgaGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGtttttt                                                                                                                       <  1:916408/149‑1 (MQ=255)
                            gacgaGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGtttttt                                                                                                                       >  2:435045/1‑149 (MQ=255)
                              cgaGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCTGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGATTTTTCa                                                                                                                     <  4:680037/149‑1 (MQ=255)
                              cgaGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCa                                                                                                                     <  4:58232/149‑1 (MQ=255)
                              cgaGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCa                                                                                                                     >  4:622088/1‑149 (MQ=255)
                               gaGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCAc                                                                                                                    <  2:12057/149‑1 (MQ=255)
                                 gcagcGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCt                                                                                                                  >  2:1094197/1‑149 (MQ=255)
                                   agcGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTgg                                                                                                                >  4:949777/1‑149 (MQ=255)
                                     cgccgcCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCt                                                                                                              <  4:625086/149‑1 (MQ=255)
                                       ccgccgTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGt                                                                                                            >  3:125051/1‑149 (MQ=255)
                                          ccgTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGa                                                                                                         <  3:989818/149‑1 (MQ=255)
                                           cgatttGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAg                                                                                                        >  2:630973/4‑149 (MQ=255)
                                           cgTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAg                                                                                                        <  4:318371/149‑1 (MQ=255)
                                           cgTTTGGTGAATCGCCAGCGTTAAAGCGTTGGTGCCCCCGGCGGTAACGGTGCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCGGATAAGACTGCACGGGCGCGCGAACGCCCGCCAGGGCGAAATGTTTGTACACCTGGCCGTCGAg                                                                                                        <  2:888606/149‑1 (MQ=255)
                                            gTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAgg                                                                                                       <  2:657707/149‑1 (MQ=255)
                                               ttGTCAATCGCCAGCGTTAAAGCGGTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTATCGAGGGCa                                                                                                    >  3:913169/1‑149 (MQ=255)
                                                   cAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGc                                                                                                >  1:492836/1‑149 (MQ=255)
                                                    aaTCGCCAGCGTTAAAGCTTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCTGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGTTCGAAATGTTTTTTCACCTGGCCGTCGAGGGCAAAGCg                                                                                               >  4:337223/1‑149 (MQ=255)
                                                    aaTCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCg                                                                                               >  1:563482/1‑149 (MQ=255)
                                                    aaTCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCg                                                                                               >  3:1059936/1‑149 (MQ=255)
                                                        gCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTAcc                                                                                           >  1:1063908/1‑149 (MQ=255)
                                                         ccAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCg                                                                                          >  2:593048/1‑149 (MQ=255)
                                                          cAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGt                                                                                         <  4:1070366/149‑1 (MQ=255)
                                                          cAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGt                                                                                         >  3:897011/1‑149 (MQ=255)
                                                           aGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTc                                                                                        <  4:913169/149‑1 (MQ=255)
                                                                 aaaGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGt                                                                                  >  1:312033/1‑149 (MQ=255)
                                                                  gaGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGtt                                                                                 <  2:203243/148‑1 (MQ=255)
                                                                  aaGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGtt                                                                                 <  3:350455/149‑1 (MQ=255)
                                                                   aGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACGTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGttt                                                                                <  2:641895/149‑1 (MQ=255)
                                                                   aGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAAcgcc                                                                                                                                            >  2:156735/1‑89 (MQ=255)
                                                                    gCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTg                                                                               >  2:1011566/1‑149 (MQ=255)
                                                                     cGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGa                                                                              <  4:872961/149‑1 (MQ=255)
                                                                     cGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGa                                                                              <  1:1011566/149‑1 (MQ=255)
                                                                     cGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGa                                                                              <  3:562753/149‑1 (MQ=255)
                                                                      gTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTAAg                                                                             <  1:164789/149‑1 (MQ=255)
                                                                       ttGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCTCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGt                                                                            >  3:608148/1‑149 (MQ=255)
                                                                          gTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAACGCGTACCGTCCACTGTTTGAGTGGc                                                                         >  4:333406/1‑149 (MQ=255)
                                                                          gTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGc                                                                         >  4:91808/1‑149 (MQ=255)
                                                                              ccccGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCg                                                                     >  1:1009264/1‑149 (MQ=255)
                                                                                ccGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCgtg                                                                   >  2:705843/1‑149 (MQ=255)
                                                                                ccGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCgtg                                                                   >  1:550978/1‑149 (MQ=255)
                                                                                     ggTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGaa                                                              <  2:364253/149‑1 (MQ=255)
                                                                                        aaCGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGa                                                           >  1:377042/1‑149 (MQ=255)
                                                                                          cGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGaaa                                                         <  1:454712/149‑1 (MQ=255)
                                                                                          cGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGaaa                                                         <  4:729153/149‑1 (MQ=255)
                                                                                          cGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGaaa                                                         <  1:162835/149‑1 (MQ=255)
                                                                                            gTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAaca                                                       <  4:971654/149‑1 (MQ=255)
                                                                                             ttCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAacac                                                      >  2:403901/1‑149 (MQ=255)
                                                                                             ttCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAacac                                                      >  1:1032549/1‑149 (MQ=255)
                                                                                                aGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCa                                                   <  1:435045/149‑1 (MQ=255)
                                                                                                  ccGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAgt                                                 <  2:492836/149‑1 (MQ=255)
                                                                                                  ccGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAgt                                                 <  3:725853/149‑1 (MQ=255)
                                                                                                                ggAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTAtt                                   >  1:316844/1‑149 (MQ=255)
                                                                                                                ggAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTAtt                                   <  2:550978/149‑1 (MQ=255)
                                                                                                                       gcagcaCCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTGCTTAAg                            <  1:331604/149‑1 (MQ=255)
                                                                                                                        cagcaCCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTGCTTAAGc                           <  2:377042/149‑1 (MQ=255)
                                                                                                                          gcaCCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTGCTTAAGCCg                         >  1:117600/1‑149 (MQ=255)
                                                                                                                              cTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAAc                                                             <  3:664237/109‑1 (MQ=255)
                                                                                                                              cTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTGCTTAAGCCGACaa                     >  3:820920/1‑149 (MQ=255)
                                                                                                                              cTGATAAGTCTGCACTGGCGCGCGAACGACCGCCAGGTCGAAATGTATTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTTTTTGAGTGGCAGCGTGGTGAAc                                                             >  4:664237/1‑109 (MQ=255)
                                                                                                                                gATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTGCTTAAGCCGACAATc                   >  1:896416/1‑149 (MQ=255)
                                                                                                                                 aTAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTGCTTAAGCCGACAATcc                  >  2:192678/1‑149 (MQ=255)
                                                                                                                                 aTAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTGCTTAAACCGACAATcc                  <  1:669086/149‑1 (MQ=255)
                                                                                                                                    aGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTGCTTAAGCCGACAATCCCgg               <  4:369492/149‑1 (MQ=255)
                                                                                                                                     gTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTGCTTAAGCCGACAATCCCGGc              >  4:792236/1‑149 (MQ=255)
                                                                                                                                        tGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTGCTTAAGCCGACAATCCCGGCaaa           >  2:708386/1‑149 (MQ=255)
                                                                                                                                        tGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTGCTTAAGCCGACAATCCCGGCaaa           >  3:404703/1‑149 (MQ=255)
                                                                                                                                         gaccTGGCGGGCGTTCGCGCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTGCTTAAGCCGACAATCCCGGCaaaa          >  3:921236/4‑149 (MQ=255)
                                                                                                                                         gCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTGCTTAAGCCGACAATCCCGGCaaaa          <  2:358362/149‑1 (MQ=255)
                                                                                                                                               gcgcgcGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTGCTTAAGCCGACAATCCCGGCAAAATTCGGt    <  2:485881/149‑1 (MQ=255)
                                                                                                                                                cgcgcgAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTGCTTAAGCCGACAATCCCGGCAAAATTCGGtt   <  4:143455/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTGCTTAAGCCGACAATCCCGGCAAAATTCGGTTC  >  NC_000913/843192‑843485

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