Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,800,517 T→A 5.7% intergenic (+42/‑206) 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,5170TA 5.7% 890.7 / 31.4 259intergenic (+42/‑206)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 (117/127);  new base A (8/7);  total (125/134)
Fisher's exact test for biased strand distribution p-value = 7.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

Only 100 of 264 total aligned reads displayed.

GCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACAGAG  >  NC_000913/2800372‑2800665
                                                                                                                                                 |                                                                                                                                                    
     gcgGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCgg                                                                                                                                              <  3:147650/149‑1 (MQ=255)
       ggATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTAACCCCGAAGACAGGCTTCgggg                                                                                                                                            <  2:791315/149‑1 (MQ=255)
        gATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGt                                                                                                                                           >  2:1703308/1‑149 (MQ=255)
           tCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTc                                                                                                                                        >  3:355376/1‑149 (MQ=255)
             ttttGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCtt                                                                                                                                      >  4:1004066/1‑149 (MQ=255)
               ttGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCtttt                                                                                                                                    >  2:1161235/1‑149 (MQ=255)
                tGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCttttt                                                                                                                                   <  2:364111/149‑1 (MQ=255)
                tGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCttttt                                                                                                                                   <  1:1277738/149‑1 (MQ=255)
                tGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCttttt                                                                                                                                   <  3:649637/149‑1 (MQ=255)
                  ttCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGc                                                                                                                                 >  4:641523/1‑149 (MQ=255)
                  ttCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGc                                                                                                                                 >  4:305993/1‑149 (MQ=255)
                    cGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCt                                                                                                                               <  4:1376751/149‑1 (MQ=255)
                      tGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTgg                                                                                                                             <  3:2112052/149‑1 (MQ=255)
                      tGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTgg                                                                                                                             >  4:1351222/1‑149 (MQ=255)
                       gAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGt                                                                                                                            >  4:1401374/1‑149 (MQ=255)
                        aaaGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTg                                                                                                                           >  3:699305/1‑149 (MQ=255)
                           gACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCttttt                                                                                                                                   <  3:1788945/138‑1 (MQ=255)
                           gACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGCGAAt                                                                                                                        >  4:274838/1‑149 (MQ=255)
                              ccTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATATGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATcac                                                                                                                     >  1:437114/1‑149 (MQ=255)
                              ccTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATcac                                                                                                                     >  4:793843/1‑149 (MQ=255)
                              ccTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATcac                                                                                                                     >  2:2168729/1‑149 (MQ=255)
                              ccTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATcac                                                                                                                     >  1:959338/1‑149 (MQ=255)
                              ccTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATcac                                                                                                                     >  1:2460655/1‑149 (MQ=255)
                                 ggTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACaaa                                                                                                                  <  1:748063/149‑1 (MQ=255)
                                 ggTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACaaa                                                                                                                  <  2:437114/149‑1 (MQ=255)
                                      aGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATcc                                                                                                             <  3:1067827/149‑1 (MQ=255)
                                       gCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATccc                                                                                                            <  1:179865/149‑1 (MQ=255)
                                       gCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATccc                                                                                                            >  3:1717516/1‑149 (MQ=255)
                                        cGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCCTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATcccc                                                                                                           <  4:417611/149‑1 (MQ=255)
                                           gggTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTa                                                                                                        >  4:234729/1‑149 (MQ=255)
                                             gTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTAc                                                                                                       <  1:894006/148‑1 (MQ=255)
                                             gTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTAc                                                                                                       >  2:894006/1‑148 (MQ=255)
                                              tACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCggg                                                                                                                                             <  1:2086213/109‑1 (MQ=255)
                                                cagcagCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCg                                                                                                   <  3:1356956/149‑1 (MQ=255)
                                                  gcagcTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTc                                                                                                 >  3:1558309/1‑149 (MQ=255)
                                                   cagcTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCa                                                                                                >  3:1318066/1‑149 (MQ=255)
                                                    agcTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCggg                                                                                                                                             <  4:408024/103‑1 (MQ=255)
                                                      cTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGc                                                                                             <  3:1908060/149‑1 (MQ=255)
                                                        gCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTc                                                                                           <  3:274838/149‑1 (MQ=255)
                                                        gCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTc                                                                                           <  2:1346950/149‑1 (MQ=255)
                                                          ggtaggtaTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCat                                                                                         >  2:854631/1‑149 (MQ=255)
                                                             aggtaTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATc                                                                                      <  3:414598/149‑1 (MQ=255)
                                                              ggtaTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATcc                                                                                     >  4:1594888/1‑149 (MQ=255)
                                                              ggtaTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATcc                                                                                     <  4:699305/149‑1 (MQ=255)
                                                                taTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAg                                                                                   <  1:1687355/149‑1 (MQ=255)
                                                                taTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAg                                                                                   <  1:2152218/149‑1 (MQ=255)
                                                                  tttttATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAggg                                                                                 >  2:2110832/1‑149 (MQ=255)
                                                                  tgtttATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAggg                                                                                 >  1:584984/3‑149 (MQ=255)
                                                                         cGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCGTTGACCCCGAATCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTc                                                                          <  2:1491701/149‑1 (MQ=255)
                                                                            cgcACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGAc                                                                       <  2:718457/149‑1 (MQ=255)
                                                                                cTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACt                                                                   <  1:299986/149‑1 (MQ=255)
                                                                                 tGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATccc                                                                                                            <  2:1836992/107‑1 (MQ=255)
                                                                                 tGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTATTGCGTGGTGAATCACAATAATCCCCCTACCCCGTCACGCTCATATGCAGGGTAATTTCTACCACTa                                                                  >  1:1944141/1‑149 (MQ=255)
                                                                                   ttAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTAtt                                                                <  4:1269986/149‑1 (MQ=255)
                                                                                         aTGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAGAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCtat                                                          >  3:1322593/1‑149 (MQ=255)
                                                                                          tGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCtata                                                         >  3:886759/1‑149 (MQ=255)
                                                                                             gcAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTatatat                                                      >  3:530939/1‑149 (MQ=255)
                                                                                             gcAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTatatat                                                      >  4:2271897/1‑149 (MQ=255)
                                                                                               aaaTCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATtg                                                    <  1:1324680/149‑1 (MQ=255)
                                                                                                     tAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGt                                              <  4:1818505/149‑1 (MQ=255)
                                                                                                      aGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGtt                                             <  1:938134/149‑1 (MQ=255)
                                                                                                      aGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGtt                                             <  3:181350/149‑1 (MQ=255)
                                                                                                        tGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGa                                           <  2:10366/149‑1 (MQ=255)
                                                                                                        tGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGa                                           >  2:741433/1‑149 (MQ=255)
                                                                                                          cAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAAt                                         <  4:892312/149‑1 (MQ=255)
                                                                                                           aaaaaTGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATc                                        >  3:797608/1‑149 (MQ=255)
                                                                                                              aaTGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCttt                                     <  2:1250998/149‑1 (MQ=255)
                                                                                                               aTGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATccccc                                                                                                          <  1:1173637/79‑1 (MQ=255)
                                                                                                               aTGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCtttt                                    <  2:2231123/149‑1 (MQ=255)
                                                                                                                tGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCttttt                                   <  1:2332438/149‑1 (MQ=255)
                                                                                                                 gATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGtt                                             <  1:1022861/138‑1 (MQ=255)
                                                                                                                 gATAATAAATACGCGTCTTTAACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGtt                                             >  2:1022861/1‑138 (MQ=255)
                                                                                                                    aataaataCGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCaa                               <  3:234729/149‑1 (MQ=255)
                                                                                                                    aataaataCGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCaa                               <  1:2279117/149‑1 (MQ=255)
                                                                                                                     ataaataCGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAAc                              <  2:262916/149‑1 (MQ=255)
                                                                                                                      taaataCGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTGCTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACt                             >  1:1957015/1‑149 (MQ=255)
                                                                                                                       aaataCGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTa                            <  3:2474587/149‑1 (MQ=255)
                                                                                                                       aaataCGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTa                            >  2:1395105/1‑149 (MQ=255)
                                                                                                                       aaataCGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTa                            <  2:2280193/149‑1 (MQ=255)
                                                                                                                       aaataCGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTa                            <  3:2206414/149‑1 (MQ=255)
                                                                                                                       aaataCGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTa                            >  2:938512/1‑149 (MQ=255)
                                                                                                                         ataCGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACa                          <  4:1484511/149‑1 (MQ=255)
                                                                                                                              cgTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAg                     <  1:1657991/149‑1 (MQ=255)
                                                                                                                              cgTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAg                     <  4:358824/149‑1 (MQ=255)
                                                                                                                                tCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCttttt                                   <  4:1115394/133‑1 (MQ=255)
                                                                                                                                tCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCttttt                                   >  3:1115394/1‑133 (MQ=255)
                                                                                                                                aaagaaaCCCCGAAGACAGGCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCtt                                      >  4:1321300/7‑130 (MQ=255)
                                                                                                                                  tttGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTc                                                                          <  3:1888507/92‑1 (MQ=255)
                                                                                                                                  tttGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTc                                                                          >  4:1888507/1‑92 (MQ=255)
                                                                                                                                  tttGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTAtc                 <  2:1957015/149‑1 (MQ=255)
                                                                                                                                       ccccGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGt            >  1:2201385/1‑149 (MQ=255)
                                                                                                                                       ccccGAAGACAGGCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGt                                              >  4:1921180/1‑115 (MQ=255)
                                                                                                                                        cccGAAGACAGGCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTAACGCCCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTa           >  1:27410/1‑149 (MQ=255)
                                                                                                                                         ccGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTc                                                                          <  2:997587/85‑1 (MQ=255)
                                                                                                                                          cGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACa                          <  3:1377862/132‑1 (MQ=255)
                                                                                                                                                cTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACaga   <  4:2063409/149‑1 (MQ=255)
                                                                                                                                                cTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACaga   >  2:1101234/1‑149 (MQ=255)
                                                                                                                                                 tGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCTGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACagag  >  4:1721609/1‑149 (MQ=255)
                                                                                                                                                 tGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACagag  <  1:2416094/149‑1 (MQ=255)
                                                                                                                                                 tGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACagag  >  1:1557734/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACAGAG  >  NC_000913/2800372‑2800665

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