Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,800,519 T→G 5.5% intergenic (+44/‑204) ygaM → / → nrdH putative membrane‑anchored DUF883 family ribosome‑binding protein/hydrogen donor for NrdEF electron transport system; glutaredoxin‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,800,5190TG 5.5% 375.5 / 6.7 128intergenic (+44/‑204)ygaM/nrdHputative membrane‑anchored DUF883 family ribosome‑binding protein/hydrogen donor for NrdEF electron transport system; glutaredoxin‑like protein
Reads supporting (aligned to +/- strand):  ref base T (64/56);  new base G (6/2);  total (70/58)
Fisher's exact test for biased strand distribution p-value = 2.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 128 total aligned reads displayed.

GCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACAGAGA  >  NC_000913/2800377‑2800666
                                                                                                                                              |                                                                                                                                                   
gcgGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCggg                                                                                                                                              <  1:724208/150‑1 (MQ=255)
gcgGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCggg                                                                                                                                              >  1:1582263/1‑150 (MQ=255)
 cgGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCgggg                                                                                                                                             <  1:1099169/150‑1 (MQ=255)
   gATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGtt                                                                                                                                           >  2:581165/1‑150 (MQ=255)
       cTTTTGTTCGTGAAAGACCCTGGTGTAGCGTCGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATTATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCTGGGTTTCtt                                                                                                                                       >  1:1248752/1‑150 (MQ=255)
           tGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACTGCTGATGATGTAGGTATTTTTATCGGCGCTCTGTTAAGCATTCGCAATTCGTTGTGCAAAAATGATAATTAATACGCGTCTTTGACCCCGAAGCCTGTCTTGGGGGTTTCTTTTTg                                                                                                                                   <  2:14233/150‑1 (MQ=255)
            gTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGc                                                                                                                                  <  2:816609/150‑1 (MQ=255)
             ttggTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGcc                                                                                                                                 >  2:1966257/4‑150 (MQ=255)
              tCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCt                                                                                                                                >  1:1557687/1‑150 (MQ=255)
              tCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCt                                                                                                                                <  1:194718/150‑1 (MQ=255)
               cGTGAAAGACCCTGGGGCAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTTACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTg                                                                                                                               <  2:241104/150‑1 (MQ=255)
                gTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTgg                                                                                                                              >  2:1683136/1‑150 (MQ=255)
                 tGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGt                                                                                                                             <  2:744490/150‑1 (MQ=255)
                 tGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGt                                                                                                                             >  1:418189/1‑150 (MQ=255)
                   aaaGACCCTGGTGTAGAGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGAAAATAAATACGCGGCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGa                                                                                                                           <  2:1198959/150‑1 (MQ=255)
                    aaGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCGTGGTGaa                                                                                                                          <  1:427187/150‑1 (MQ=255)
                      gACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATc                                                                                                                        <  1:1477097/150‑1 (MQ=255)
                        cccTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTAAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGCTTTCTTTTTTCCTGGTTAATaac                                                                                                                      >  1:273158/1‑150 (MQ=255)
                         ccTGTTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGAATGCTAAACTCGTAGTGCAAAAATGATAATAAATACGCGACTTTGACCCCTAAGCATTTCTTCGGGTTTTCTTTTTGCCTGGTGATTcaca                                                                                                                     >  2:684588/1‑150 (MQ=255)
                             gtgtAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACaaaaa                                                                                                                 <  1:1986164/150‑1 (MQ=255)
                                tAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATcc                                                                                                              <  2:671828/150‑1 (MQ=255)
                                 cgCGGGGGTACAGCAGCGGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTACTGCAAAAATGATAATAAATCCGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATccc                                                                                                             <  2:308533/149‑1 (MQ=255)
                                 aGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATccc                                                                                                             >  2:107124/1‑150 (MQ=255)
                                  gCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTTGTGAATCACAAAAATcccc                                                                                                            >  2:1813275/1‑150 (MQ=255)
                                      gggTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCACGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTAc                                                                                                        >  1:992422/1‑150 (MQ=255)
                                      gggTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAAAAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGCTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTAc                                                                                                        >  2:111265/1‑150 (MQ=255)
                                       ggTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTAcc                                                                                                       >  1:257337/1‑150 (MQ=255)
                                           cagcagCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGt                                                                                                   >  2:1309613/1‑150 (MQ=255)
                                               agcTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACg                                                                                               <  2:36167/150‑1 (MQ=255)
                                                 cTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAAGCACAAAAATCCCCCTACCCCGTCACGCt                                                                                             <  2:1037109/150‑1 (MQ=255)
                                                   gCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCa                                                                                           >  1:849107/1‑150 (MQ=255)
                                                   gCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCa                                                                                           <  2:385608/150‑1 (MQ=255)
                                                             tttttATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGt                                                                                 >  2:884671/1‑150 (MQ=255)
                                                              ttttATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCCCGCTCATATCCAGGGTa                                                                                >  2:304065/1‑150 (MQ=255)
                                                              ttttATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTa                                                                                <  1:640692/150‑1 (MQ=255)
                                                                 tATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACACCGTCACGCTCATATCCAGGGTAAtt                                                                             >  1:141647/1‑150 (MQ=255)
                                                                     ggCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCTGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCCCGCTCATATCCAGGGTAATTTCGa                                                                         >  2:82840/1‑150 (MQ=255)
                                                                      gcgcACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGAc                                                                        <  2:38434/150‑1 (MQ=255)
                                                                       cgcACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGAcc                                                                       <  2:475910/150‑1 (MQ=255)
                                                                         cACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCAc                                                                     <  1:1487172/150‑1 (MQ=255)
                                                                             gTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTAtt                                                                 >  1:523644/1‑150 (MQ=255)
                                                                             gTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTAtt                                                                 >  1:394236/1‑150 (MQ=255)
                                                                             gTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTAtt                                                                 >  1:368926/1‑150 (MQ=255)
                                                                             gTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATAACCAGGGTAATTTCGACCACTAtt                                                                 >  1:6107/1‑150 (MQ=255)
                                                                                aaGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGGCGTTGACCCCGAAGCCTGGCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGc                                                              <  1:166746/150‑1 (MQ=255)
                                                                                 aGCATGCGCAAATCGCCGTGCAAAAATGAAAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGCCCACTATTTGCt                                                             <  2:324280/150‑1 (MQ=255)
                                                                                   cATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCtat                                                           <  2:1697960/150‑1 (MQ=255)
                                                                                   cATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCtat                                                           <  1:1878805/150‑1 (MQ=255)
                                                                                   cATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGAGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCtat                                                           <  1:30680/150‑1 (MQ=255)
                                                                                     tGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCtatat                                                         >  2:1989419/1‑150 (MQ=255)
                                                                                     tGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCtatat                                                         >  1:57066/1‑150 (MQ=255)
                                                                                     tGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCtatat                                                         <  2:596907/150‑1 (MQ=255)
                                                                                      gcgcAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCtatata                                                        <  1:935725/150‑1 (MQ=255)
                                                                                       cgcAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTatatat                                                       >  2:478441/1‑150 (MQ=255)
                                                                                       cgcAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTatatat                                                       >  1:1524129/1‑150 (MQ=255)
                                                                                       cgcAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTatatat                                                       >  2:1292005/1‑150 (MQ=255)
                                                                                       cgcAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGACACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTatatat                                                       >  2:216509/1‑150 (MQ=255)
                                                                                          aaaTCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATtgt                                                    >  1:1603733/1‑150 (MQ=255)
                                                                                              cGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTgg                                                <  1:427067/150‑1 (MQ=255)
                                                                                               gTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGt                                               >  2:383593/1‑150 (MQ=255)
                                                                                                 aGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAATCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTg                                             <  1:605970/150‑1 (MQ=255)
                                                                                                 aGTGCAAAAATGATAATAAATACGCGTCTTTGACACCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGTATCAAAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGAGCACTATTTGCTATATATTGTGTGGTTg                                             >  1:252785/1‑150 (MQ=255)
                                                                                                   tGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGat                                           >  1:1320199/1‑149 (MQ=255)
                                                                                                   tGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGaa                                           >  1:1313628/1‑150 (MQ=255)
                                                                                                   tGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGaa                                           >  1:838160/1‑150 (MQ=255)
                                                                                                   tGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGaa                                           <  2:641446/150‑1 (MQ=255)
                                                                                                    gCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAAt                                          >  1:825305/1‑150 (MQ=255)
                                                                                                    gCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAAt                                          >  1:1273513/1‑150 (MQ=255)
                                                                                                    gCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCGTTTTGCCTGGTGAATCACAAAAATCCCCCTCCCCCGGCACGCTCATATCCAGGGTAATGTCGGCCACTATTTGCTAGATATTGTGTGGATGAAt                                          >  1:408304/1‑150 (MQ=255)
                                                                                                     cAAAAATGATTATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGATCATATCCAGGGTAATTTCGACAACTATTTGCTATATATTGAGTGGTTGAATc                                         >  1:27753/1‑150 (MQ=255)
                                                                                                      aaaaaTGATAATAAATACGCGTCTTTGACCACGAAGCCTGTCTTCGCGGTTTCTTTTTGACTGGTGAATCACAAAAATCCCACTACCCCGTCACGCTCAAATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGGGGAATcg                                        >  1:1508044/1‑149 (MQ=255)
                                                                                                         aaTGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGCCCAATATTTGCTATATATTGTGTGGTTGAATCtttt                                     >  2:267678/1‑150 (MQ=255)
                                                                                                           tGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCtttttt                                   <  1:1547397/150‑1 (MQ=255)
                                                                                                               aataaataCGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCaat                               <  2:825305/150‑2 (MQ=255)
                                                                                                                ataaataCGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAAca                              >  2:866675/1‑149 (MQ=255)
                                                                                                                      aCGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCt                        <  1:1382039/150‑1 (MQ=255)
                                                                                                                      aCGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCt                        <  1:1038285/150‑1 (MQ=255)
                                                                                                                       cgcgTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTa                       >  2:1889040/1‑150 (MQ=255)
                                                                                                                        gcgTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTTGTGAATCACAAAAATCCCCCTACGCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAg                      <  1:1075626/150‑1 (MQ=255)
                                                                                                                        gcgTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAACATCCCCCTACCCCGTCACGCTCATATCCAGCGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAg                      >  1:1753211/1‑150 (MQ=255)
                                                                                                                         cgTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGt                     <  2:578250/150‑1 (MQ=255)
                                                                                                                          aaaagaaaCCCCGAAGACAGGCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTa                    >  1:2008202/8‑150 (MQ=255)
                                                                                                                            cTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTTCATCTAGTAtc                  >  2:1865845/1‑150 (MQ=255)
                                                                                                                             tttGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAATTCCCCCTACCCCGTCACGCGCATATCCAGGGGATTTTCGACCACTATTGGCTATATATTGTGTGGTTGAAGCTTTTTTCAACTACATCTAGTTtct                 >  1:1379531/1‑150 (MQ=255)
                                                                                                                             tttGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTAtct                 >  2:555481/1‑150 (MQ=255)
                                                                                                                               tGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATctct               <  1:148960/150‑1 (MQ=255)
                                                                                                                                gACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTg              >  1:890408/1‑150 (MQ=255)
                                                                                                                                aaCCCCGAAGACAGGCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCACTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTAGCTATATATTGAGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTg              >  1:1916359/2‑150 (MQ=255)
                                                                                                                                 aCCCCGAAGACAGGCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGt             >  1:1852109/1‑150 (MQ=255)
                                                                                                                                  ccccGAAGACAGGCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTGCGACCCCTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGtt            >  2:1291562/1‑149 (MQ=255)
                                                                                                                                  ccccGAAGACAGGCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAAGCCCCCTCCCCCGTCCCGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGGGTGGTTGAATCTTTTTTCAACTACCTCTAGTATCTCTgga            >  1:775368/1‑148 (MQ=255)
                                                                                                                                   cccGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCCTATCCAGGGTAATTTCGACCACTAATTGCTATATATAGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTAt           >  1:297106/1‑150 (MQ=255)
                                                                                                                                    ccGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATc          >  1:1078563/1‑150 (MQ=255)
                                                                                                                                       aaGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAAc       <  2:380879/150‑1 (MQ=255)
                                                                                                                                       aaGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGGTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAAc       <  1:1966257/150‑1 (MQ=255)
                                                                                                                                        aGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACa      <  2:894059/150‑1 (MQ=255)
                                                                                                                                        aGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACa      <  1:404926/150‑1 (MQ=255)
                                                                                                                                          ccTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCATCTACATCTAGTATCTCTGTATCAACaga    >  1:986878/1‑150 (MQ=255)
                                                                                                                                          ccTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACaga    >  1:622823/1‑150 (MQ=255)
                                                                                                                                           cTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACagag   <  2:57066/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
GCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACAGAGA  >  NC_000913/2800377‑2800666

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