Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,146,625 T→C 5.4% noncoding (37/169 nt) psrD → novel Novel sRNA, function unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,146,6250TC 5.4% 375.0 / 8.0 113noncoding (37/169 nt)psrDnovel Novel sRNA, function unknown
Reads supporting (aligned to +/- strand):  ref base T (58/48);  new base C (3/3);  total (61/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

Only 100 of 122 total aligned reads displayed.

TGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACA  >  NC_000913/1146476‑1146774
                                                                                                                                                     |                                                                                                                                                     
tGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATtgat                                                                                                                                                       <  1:3116895/150‑1 (MQ=255)
  aTTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATtgatga                                                                                                                                                     <  2:2423414/150‑1 (MQ=255)
  aTTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATtgatga                                                                                                                                                     >  1:2400473/1‑150 (MQ=255)
  aTTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATtgatga                                                                                                                                                     <  2:517517/150‑1 (MQ=255)
   ttCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTgatgat                                                                                                                                                    >  2:1099395/1‑150 (MQ=255)
   ttCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTgatgat                                                                                                                                                    >  2:2186175/1‑150 (MQ=255)
    tcgtcgCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATa                                                                                                                                                   >  1:2254979/1‑150 (MQ=255)
    tcgtcgCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGATAAAATAATTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATa                                                                                                                                                   <  1:2220380/150‑1 (MQ=255)
      gtcgCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGc                                                                                                                                                 <  1:2639775/150‑1 (MQ=255)
       tcgCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGcc                                                                                                                                                >  1:2938282/1‑150 (MQ=255)
         gCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCAt                                                                                                                                              >  2:1102272/1‑150 (MQ=255)
               tGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTaaa                                                                                                                                        <  2:862064/150‑1 (MQ=255)
                      tcgtcgACCTCGGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATgcg                                                                                                                                 <  2:1482747/150‑1 (MQ=255)
                        gtcgACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCt                                                                                                                               >  2:3251417/1‑150 (MQ=255)
                         tcgACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCtt                                                                                                                              >  1:1359878/1‑150 (MQ=255)
                         tcgACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCtt                                                                                                                              >  1:2290426/1‑150 (MQ=255)
                          cgACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTc                                                                                                                             <  2:57278/150‑1 (MQ=255)
                                  ggTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCa                                                                                                                     <  2:2264850/150‑1 (MQ=255)
                                   gTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCGGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAg                                                                                                                    >  2:2321770/1‑150 (MQ=255)
                                   gTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAg                                                                                                                    >  1:2110674/1‑150 (MQ=255)
                                   gTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAg                                                                                                                    >  2:274382/1‑150 (MQ=255)
                                     gctgcACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGt                                                                                                                  <  1:1936081/150‑1 (MQ=255)
                                       tgcACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTg                                                                                                                >  1:2094269/1‑150 (MQ=255)
                                        gcACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGa                                                                                                               >  1:893858/1‑150 (MQ=255)
                                         cacaTTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGAt                                                                                                              >  2:3116300/1‑150 (MQ=255)
                                            aTTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGAtgtt                                                                                                           >  1:709866/1‑150 (MQ=255)
                                             ttCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGAtgttt                                                                                                          >  1:639475/1‑150 (MQ=255)
                                                aaaaGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGAtgtttgtt                                                                                                       <  2:847697/150‑1 (MQ=255)
                                                   agagATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGACGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATc                                                                                                    >  2:2602207/1‑150 (MQ=255)
                                                   agagATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGACGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATc                                                                                                    >  2:2602290/1‑150 (MQ=255)
                                                   agagATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGACGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATc                                                                                                    >  2:2005414/1‑150 (MQ=255)
                                                       aTTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCa                                                                                                <  2:2400473/150‑1 (MQ=255)
                                                            cAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGTCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTTTTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGaa                                                                                           <  1:2016117/150‑1 (MQ=255)
                                                            cAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGaa                                                                                           >  2:866249/1‑150 (MQ=255)
                                                            cAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGaa                                                                                           <  2:1910515/150‑1 (MQ=255)
                                                            cAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGaa                                                                                           <  2:1722196/150‑1 (MQ=255)
                                                            cAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGaa                                                                                           >  2:1656607/1‑150 (MQ=255)
                                                            cAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTAAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGaa                                                                                           >  2:866089/1‑150 (MQ=255)
                                                              gTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACg                                                                                         >  2:2144314/1‑150 (MQ=255)
                                                                ttttCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGaa                                                                                       <  2:1453204/150‑1 (MQ=255)
                                                                   tCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAAt                                                                                    <  2:823581/150‑1 (MQ=255)
                                                                   tCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAAt                                                                                    <  1:2888460/150‑1 (MQ=255)
                                                                    ccAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATa                                                                                   >  1:1920012/1‑150 (MQ=255)
                                                                     cAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATaa                                                                                  >  2:643257/1‑150 (MQ=255)
                                                                     cAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATaa                                                                                  >  2:1223683/1‑150 (MQ=255)
                                                                      agaagaGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATaaa                                                                                 >  1:2347878/1‑150 (MQ=255)
                                                                      agaagaGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATaaa                                                                                 >  1:68092/1‑150 (MQ=255)
                                                                       gaagaGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAg                                                                                <  1:2716479/150‑1 (MQ=255)
                                                                       gaagaGCGCGACGACATGGCGAGGTGGAGGCTAAAATAAGTTTAGTCATTTTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAg                                                                                <  1:2431313/150‑1 (MQ=255)
                                                                         agagCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCa                                                                              <  2:1180857/150‑1 (MQ=255)
                                                                         agagCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCa                                                                              <  1:385184/150‑1 (MQ=255)
                                                                                gACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAAccc                                                                       >  1:140454/1‑150 (MQ=255)
                                                                                gACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAAccc                                                                       >  1:2125442/1‑150 (MQ=255)
                                                                                    ccATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGTACGAAAATAAAGCAGTAACCCGCaa                                                                   >  1:489000/1‑150 (MQ=255)
                                                                                    ccATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCaa                                                                   >  1:1390680/1‑150 (MQ=255)
                                                                                         gCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAAtgtgt                                                              <  1:1307127/150‑1 (MQ=255)
                                                                                         gCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAAtgtgt                                                              <  2:2296816/150‑1 (MQ=255)
                                                                                         gCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAAtgtgt                                                              <  1:3196782/150‑1 (MQ=255)
                                                                                         gCGAGGTGGAGGCTAAAATAAGTTTAGCCATATTTTTTTCCATCAGATATAGCGTATTGAAGATAGCCATTTTAAACTATGCGCTTCGATTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAAtgtgt                                                              <  2:184505/150‑1 (MQ=255)
                                                                                         gCGAGGTGGAGCCTAAAATAAGTTTACCCAAATTATTTTCCATCAGATATAGCCTAATGATGATAGCCATTTTAAACTATGCGCTACGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAAtgtgt                                                              <  2:184858/150‑1 (MQ=255)
                                                                                           gAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGc                                                            <  2:361618/150‑1 (MQ=255)
                                                                                           gAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGACGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGc                                                            <  1:2005414/150‑1 (MQ=255)
                                                                                           gAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGACGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGc                                                            <  1:2602207/150‑1 (MQ=255)
                                                                                           gAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGACGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGc                                                            <  1:2602290/150‑1 (MQ=255)
                                                                                                 gAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGTACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAatta                                                      >  2:719240/1‑150 (MQ=255)
                                                                                                 gAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGActta                                                      >  2:103533/1‑150 (MQ=255)
                                                                                                 gAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAatta                                                      >  2:2471914/1‑150 (MQ=255)
                                                                                                 gAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAatta                                                      >  2:103541/1‑150 (MQ=255)
                                                                                                  aGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAattat                                                     <  1:2947299/150‑1 (MQ=255)
                                                                                                  aGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAattat                                                     <  1:2946706/150‑1 (MQ=255)
                                                                                                   ggCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAattatt                                                    <  2:763666/150‑1 (MQ=255)
                                                                                                    gCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTg                                                   >  2:921867/1‑150 (MQ=255)
                                                                                                    gCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTg                                                   <  2:639475/150‑1 (MQ=255)
                                                                                                     cTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCTGGTTGATGTTTGTTATCAGAACTGAACGAAAATAAAGCAGTAACCCGCAATGAGTGCGAATTATTgg                                                  >  1:3047632/1‑150 (MQ=255)
                                                                                                          aTAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCaaaa                                             <  2:3231871/150‑1 (MQ=255)
                                                                                                             aGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGc                                          >  2:1278679/1‑150 (MQ=255)
                                                                                                                  aGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACca                                     >  1:3103190/1‑150 (MQ=255)
                                                                                                                    gCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACcaca                                   >  2:1819159/1‑150 (MQ=255)
                                                                                                                     cATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAg                                  <  1:2055992/150‑1 (MQ=255)
                                                                                                                       tatTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGTTGTTTGTTATCAGCACTGAACGTAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGc                                >  2:259690/1‑150 (MQ=255)
                                                                                                                       tatTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGc                                >  2:2045556/1‑150 (MQ=255)
                                                                                                                       tatTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGc                                >  1:2231353/1‑150 (MQ=255)
                                                                                                                       tatTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGc                                >  2:1994534/1‑150 (MQ=255)
                                                                                                                       tatTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGc                                >  1:2231295/1‑150 (MQ=255)
                                                                                                                          ttttttCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGc                             >  2:757102/1‑150 (MQ=255)
                                                                                                                              ttCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAGCCTGTGGTTGCCttt                         <  1:3116300/150‑1 (MQ=255)
                                                                                                                               tCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAGCCTGTGGTTGCCtttt                        <  2:2254979/150‑1 (MQ=255)
                                                                                                                               tCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAGCCTGTGGTTGCCtttt                        <  1:253413/150‑1 (MQ=255)
                                                                                                                                ccATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCttttt                       <  2:2938282/150‑1 (MQ=255)
                                                                                                                                  aTCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTAAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCt                     <  1:2210785/150‑1 (MQ=255)
                                                                                                                                   tCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCtt                    >  2:2993002/1‑150 (MQ=255)
                                                                                                                                   tCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCtt                    >  1:1706384/1‑150 (MQ=255)
                                                                                                                                   tCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCtt                    >  2:2620911/1‑150 (MQ=255)
                                                                                                                                   tCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCtt                    >  2:2168877/1‑150 (MQ=255)
                                                                                                                                    cAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCttt                   <  2:2878244/150‑1 (MQ=255)
                                                                                                                                     aGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTg                  >  1:364082/1‑150 (MQ=255)
                                                                                                                                     aGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTg                  >  1:364134/1‑150 (MQ=255)
                                                                                                                                      gatcgaGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTTTCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGa                 <  2:3122317/145‑1 (MQ=255)
                                                                                                                                                     tgatAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACa  >  1:2146407/1‑150 (MQ=255)
                                                                                                                                                     tgatAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACa  >  1:635292/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACA  >  NC_000913/1146476‑1146774

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