Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,286,965 A→G 10.5% intergenic (+439/+101) narI → / ← rttR nitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,286,9650AG10.5% 325.9 / 16.2 133intergenic (+439/+101)narI/rttRnitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript
Reads supporting (aligned to +/- strand):  ref base A (47/62);  new base G (5/7);  total (56/77)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 134 total aligned reads displayed.

ATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGG  >  NC_000913/1286819‑1287112
                                                                                                                                                  |                                                                                                                                                   
 tCCTTCCCCGCTTGGGCAGAATCTTTGATTGCGGATTCGTTTGAGAATTACGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCt                                                                                                                                                 <  2:4064356/150‑1 (MQ=255)
 tCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCt                                                                                                                                                 <  1:2252400/150‑1 (MQ=255)
    gtcaccGCTTGGGCAGAATATAAGATTGCGGATTCGTATGAGAATTCCTGGCCTTTAGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGa                                                                                                                                              <  2:1545447/146‑1 (MQ=14)
      ccccGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGc                                                                                                                                            <  2:2511866/150‑1 (MQ=255)
        ccGCTTGGGCAGAATATTTTATTGCGGATTCGTTTGAGAATTCCGGGTCTTTTGTAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCg                                                                                                                                          <  2:1821455/150‑1 (MQ=255)
        ccGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCg                                                                                                                                          >  2:852489/1‑150 (MQ=255)
        ccGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCg                                                                                                                                          <  2:3618309/150‑1 (MQ=255)
        ccGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACCCTCTCTCGCGGACGCACGAGCCg                                                                                                                                          >  1:482450/1‑150 (MQ=255)
        ccGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTAGTTGCCGGCGGCAACGCTCTCTCGCTTTCGCTCGAGCCg                                                                                                                                          >  1:3079619/1‑150 (MQ=255)
           cTTGGGCAGAATTTTTGATTGCGGATTCGTTTGAGCATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAAc                                                                                                                                       <  2:3267070/150‑1 (MQ=255)
           cTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAAc                                                                                                                                       <  1:1689625/150‑1 (MQ=255)
                  aGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTGTCGCTGACGCTCGAGCCGAACCTTAGTc                                                                                                                                <  2:2148252/150‑1 (MQ=255)
                  aGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTc                                                                                                                                <  2:2148359/150‑1 (MQ=255)
                  aGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCATAGTc                                                                                                                                >  1:2271912/1‑150 (MQ=255)
                  aGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCACCGCTCCGCGCATCGACCTTCGGGTCGTTGCCTGCGGCAACGCTCTATCTCTGGCGCTCGGGCCGAACATTAGTc                                                                                                                                >  2:2369907/1‑150 (MQ=14)
                    aaTATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGa                                                                                                                              <  1:3704170/150‑1 (MQ=255)
                        tttGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCTCTGCGCGCATCGCCCTTCGGGTCGTTTACTGCGTCAACGCTCTCTCGCACACGCTACAGCCGAACCATACTCAAAGCt                                                                                                                          >  2:1317391/1‑150 (MQ=14)
                         ttGATTGCGGACTCGTTTGAGAATTCCTGTGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCtt                                                                                                                         <  2:2905960/150‑1 (MQ=255)
                                  gATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCtt                                                                                                                <  1:3815497/150‑1 (MQ=255)
                                  gATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCtt                                                                                                                <  1:3815311/150‑1 (MQ=255)
                                     tCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAAAGCGCTCTCCCAGTCGCTCGAGCCGAAACGAAGAAGAAGCATCAAATCCTTcac                                                                                                             >  1:3311318/1‑148 (MQ=14)
                                      cGTTTGAGAATTCCGTGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTcccc                                                                                                            <  2:1602906/150‑1 (MQ=255)
                                      cGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTcccc                                                                                                            <  1:783661/150‑1 (MQ=255)
                                        tttGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCAGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGc                                                                                                          >  1:1650660/1‑150 (MQ=255)
                                         ttGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCa                                                                                                         <  1:3507823/150‑1 (MQ=255)
                                           gagaATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCACTCGCAGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCAACCTTCCCCGCaag                                                                                                       >  2:3547937/1‑148 (MQ=255)
                                             gaATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCCACCCTCTCTCGCTGACGCACGAGCCGAACCATAGTCGAAGAGTCACAACCAACCCCGCAgggg                                                                                                     >  1:32148/1‑150 (MQ=14)
                                               aTTCCGGGGCTTTTGAAAGTGTTGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCTAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCa                                                                                                   <  1:1458930/150‑1 (MQ=255)
                                                 tCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGa                                                                                                 <  1:77681/150‑1 (MQ=255)
                                                  tcgtgccTTTAGCAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGaa                                                                                                <  2:105124/144‑1 (MQ=11)
                                                  ccGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGaa                                                                                                <  1:2517375/150‑1 (MQ=40)
                                                  ccGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGaa                                                                                                <  1:2431868/150‑1 (MQ=40)
                                                    ggggCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCCCTTCGCCCATCGGGTCGTTTCCTGCCGCAACGCTCTCTCGCTGACCCAAGAGGCGCACCTTCGTCGAAGCTTCTCATCCTTCCCCCCAGGGGCAGAata                                                                                              >  2:2306434/1‑150 (MQ=11)
                                                    ggggCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGATCTGCGCGCTTCTCCCATCGGGTCGATTCCTGCCGCAAAGCTCTCTCGCTGACCCACCAGCCGAACATTAGTAGAAGCATATCATACTACCCCGCATGGGCAGAata                                                                                              >  2:196804/1‑150 (MQ=14)
                                                    ggggCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCACATCTCCGACCATCGCCCAACTGGTACATTCCTGACCCAAAGCACTCTCGCTGACCCACCAGCCGCACCTTAGAATAATCATATCATACTTCCCCCCAGGGGCataata                                                                                              >  2:196827/1‑150 (MQ=9)
                                                     gggCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTAAGAGCTGAGCTCTTAGCCCTTCGGGTCGTTGCCGGCGGCAACGCTATCTCGCTGACGCGCGAGGCGAACCATAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAggagat                                                                                             >  2:961167/1‑145 (MQ=11)
                                                       gCTTTTGTAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAttt                                                                                           <  2:3588260/150‑1 (MQ=38)
                                                         ttttGATTGTGAGGGTGGTGGGGGAAGGATTTCTCAGCGCTGCGCGCTTCGCCGTTGGTGTCGTTTCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGa                                                                                         <  2:3472205/150‑1 (MQ=11)
                                                                   gATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTc                                                                               <  2:3375239/150‑1 (MQ=39)
                                                                     tggtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGc                                                                             <  1:2063755/150‑1 (MQ=40)
                                                                       gtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCtt                                                                           <  1:1183009/150‑1 (MQ=255)
                                                                          gtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTgag                                                                        <  1:3898082/150‑1 (MQ=255)
                                                                                 aaGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCACCGCTCTCTCGCTGACGCTCGAGCCGAACCTTCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAg                                                                 <  2:3734622/150‑1 (MQ=255)
                                                                                 aaGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAg                                                                 <  2:2616362/150‑1 (MQ=255)
                                                                                 aaGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAg                                                                 <  2:3152722/150‑1 (MQ=255)
                                                                                    gATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCGTCCCCGCGTGGGCCGAGTATTTGATTGCGGAGTCGCTTGAGAGGTCAGGGa                                                              >  2:1704814/1‑150 (MQ=14)
                                                                                       tACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTACCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACtt                                                           <  1:2663687/150‑1 (MQ=11)
                                                                                         cTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACtttt                                                         <  1:3560636/150‑1 (MQ=255)
                                                                                         cTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTGGGCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGGTTGGGAGTTCAGGGACtttt                                                         >  2:936771/1‑150 (MQ=255)
                                                                                          tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTg                                                        <  1:3315439/150‑1 (MQ=255)
                                                                                          tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAAGATTTGATTGCGGATTCGGTGGTGAGTTCAGGGACTTTTg                                                        >  2:2205175/1‑150 (MQ=255)
                                                                                             gcgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGaaa                                                     <  1:3761773/150‑1 (MQ=255)
                                                                                             gcgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGaaa                                                     <  1:3761745/150‑1 (MQ=255)
                                                                                             gcgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGaaa                                                     <  2:215753/150‑1 (MQ=255)
                                                                                               gcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGt                                                   <  1:2310000/150‑1 (MQ=255)
                                                                                               gcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGt                                                   <  2:814477/150‑1 (MQ=255)
                                                                                                 tGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCGTAGTCGAAGCTTCTCATCCATCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGGGAGTTCAGGGACTTTTGAAAGTGA                                                 >  2:3588926‑M1/1‑150 (MQ=255)
                                                                                                     cgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGT                                             <  2:2610594‑M1/150‑1 (MQ=255)
                                                                                                        ttCGCCCTTCGGGTCGTTGCCGGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTTATTGCGGATTCGCTGGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGg                                          >  2:2849753‑M1/1‑149 (MQ=255)
                                                                                                         tCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTG                                         <  1:3547937‑M1/150‑1 (MQ=255)
                                                                                                          cGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCAAGTCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGG                                        <  1:2542416‑M1/150‑1 (MQ=255)
                                                                                                           gCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGG                                       <  2:1381797‑M1/150‑1 (MQ=255)
                                                                                                            cccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGG                                      <  2:1499071‑M1/150‑1 (MQ=255)
                                                                                                            cccTTCGGGTCGTTGCCGGCGGCAACGCTCTCTCGCTGGCGCGCGAGGCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCCGAAGGTTTGTTTGCGGGTTGGGTTGAGGGTTGGGGGACTGTTGAAATGGAGGGTGgggggg                                      >  2:2472891/1‑150 (MQ=14)
                                                                                                              cTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCCAGCCGAACCATAGTCGAATCATCACATCCATCCCCCCATGGGCACAATATTTGATTTCGGAATCGCTTGAGAGTTCAGGGACTTTTTAAAGTGATGGTGGTGGGGGt                                    >  2:717320‑M1/1‑149 (MQ=255)
                                                                                                               ttCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAA                                   <  2:3548920‑M1/150‑1 (MQ=255)
                                                                                                               ttCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAA                                   <  2:3548923‑M1/150‑1 (MQ=255)
                                                                                                                tCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAG                                  <  1:2563544‑M1/150‑1 (MQ=255)
                                                                                                                tCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAG                                  >  2:3080190‑M1/1‑150 (MQ=255)
                                                                                                                tCGGGTCGTTGCCGGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCGTAGTCGAAGCTTATCATCCATACCCCCATGGGCAGAATATTTTATTTCGGATTAGGTTGTTAGTTCAGGGACTTTTTAAAGTGATGGTGGTGGGGGGAG                                  >  2:702919‑M1/1‑150 (MQ=255)
                                                                                                                  gggTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGA                                <  1:1014912‑M1/150‑1 (MQ=255)
                                                                                                                  gggTCGTTGCCTGCGGCAACCCTCTCACGCTGACGCACGAGCCCAACCGTATTCGCAGCGTCTAATACATACCCCCATGGGCAGAATATTTAATTGCGGATTAGGTGGGGAGTTCAGGGACTTTTGAAAGTGATGGTGGGGGGGGAAggt                                >  2:1431515/1‑149 (MQ=2)
                                                                                                                   ggTCGTTGCCTGCGGCAACGCACTCTCGCTGACGCTCGAGCCGAAACTTAGTCGAAGCTTCTCAAACTTACCCGCAGGGGCAGAATATTTTATTTCGGATTCGCGTGAGAGGTCAGGGACTTTTGTAAGTGAGGGGGGGGGGGGGAGGAt                               >  2:1145879/1‑150 (MQ=14)
                                                                                                                   ggTCGGTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGCACCTTCGTCGAAGCTTCTCCTCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGATTTCCGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAT                               <  1:3954846‑M1/150‑1 (MQ=255)
                                                                                                                    gTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCTAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATT                              <  2:3655247‑M1/150‑1 (MQ=255)
                                                                                                                    gTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGATCGAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGATTGAGATTTCATTGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATT                              <  1:2680623‑M1/150‑1 (MQ=255)
                                                                                                                    gTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATT                              <  1:3023507‑M1/150‑1 (MQ=255)
                                                                                                                    gTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATT                              <  1:725287‑M1/150‑1 (MQ=255)
                                                                                                                      cGTTGCGTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCGTCCCCGGAGGGGGAGAATATTTGATTGCGGATTCGCTTGAGTGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTtc                            <  1:1317391‑M1/150‑3 (MQ=255)
                                                                                                                       gTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGGTGGGGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTcga                           >  2:2105808‑M1/1‑147 (MQ=255)
                                                                                                                        ttGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGTATATTTGATTGCGGATTCGCTTGTGAGTTCAGGGACTTTTGAATGTGATGGTGGTGGGGGAAGGATTACTc                          >  2:927851/1‑150 (MQ=255)
                                                                                                                        ttGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTgaga                                                                       <  1:2232105/105‑1 (MQ=255)
                                                                                                                        ttGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTgaga                                                                       >  2:2232105/1‑105 (MQ=255)
                                                                                                                        ttGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTagta                          >  2:3750017/1‑147 (MQ=255)
                                                                                                                           ccTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAGTCGAAGCTTAACATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGCAGGATTACTCAgc                       >  2:860646/1‑150 (MQ=255)
                                                                                                                           ccTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGGTGGGGAGTTCAGGGACTTTTGAAAGTGATGGTGGGGGGGGAAGGATTcgacacg                       >  2:2223604‑M1/1‑143 (MQ=255)
                                                                                                                            cTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCCGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGGGATGGTGGTGGGGGAAGGATTACTCAgcg                      >  1:3113245/1‑150 (MQ=255)
                                                                                                                                  cAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCCGACGAGTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAACGTGATGGTGGGGGGGGAAGGATTACTCAGCGCTgcgc                >  2:1569382/1‑150 (MQ=17)
                                                                                                                                   aaCGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGaa                                                                                                >  2:1813996/1‑69 (MQ=255)
                                                                                                                                   aaCGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTgcgcg               >  2:1813980/1‑150 (MQ=255)
                                                                                                                                   aaCGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCACCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTgcgcg               >  2:1813925/1‑150 (MQ=255)
                                                                                                                                    aCGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTATCAACCTTCCCCCCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTTATGGTGGTGGGGGAAGGATTACTCAGCGCTgcgcgc              >  1:89105/1‑150 (MQ=255)
                                                                                                                                       ctctctCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgcgcttcgcgcttc           >  2:413047/1‑150 (MQ=255)
                                                                                                                                       ctctctCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTc           >  2:2104759/1‑150 (MQ=255)
                                                                                                                                         ctctcGCTGTCGCTCGAGCGGCACCATAGTCGAATACACTCATGCCTACCCTCTTTGGCATAATCTCTGTTTTCGGATTTGATTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGc         <  2:3425170/150‑1 (MQ=9)
                                                                                                                                             cGCTGGCGCTCCAGTCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATTTTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCtt     <  2:2769312/150‑1 (MQ=11)
                                                                                                                                              gCTGGCGCTCTAGCCGAACCTTCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTc    <  1:1524612/150‑1 (MQ=255)
                                                                                                                                               cTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGAGCTGCGCGCTTCGCCCTTCg   >  1:4048496/1‑150 (MQ=255)
                                                                                                                                                tGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCgg  <  2:4155006/150‑1 (MQ=255)
                                                                                                                                                tGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCgg  <  2:2107134/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
ATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGG  >  NC_000913/1286819‑1287112

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