Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,286,965 A→T 6.8% 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,9650AT 6.8% 529.8 / 4.7 206intergenic (+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 (71/107);  new base T (3/14);  total (80/126)
Fisher's exact test for biased strand distribution p-value = 1.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 210 total aligned reads displayed.

CTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGT  >  NC_000913/1286816‑1287114
                                                                                                                                                     |                                                                                                                                                     
   aTCCTTCCCCGCTTGGGCATATTATTTGATTTCGTATTCGTTTGAGAATTCCGGTGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGc                                                                                                                                                    <  2:5242031/150‑1 (MQ=14)
   aTCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGGTGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACCCACTAACGCAGACGc                                                                                                                                                    >  1:5178587/1‑150 (MQ=14)
       ttCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGa                                                                                                                                                >  2:6078444/1‑150 (MQ=255)
       ttCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGa                                                                                                                                                <  1:60849/150‑1 (MQ=255)
       ttCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTAGCAAGCGGCAACGCTCACACGCTGACGCTcca                                                                                                                                                >  1:4945414/1‑148 (MQ=14)
       ttCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCGGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCTCTCGCGCGGACCCTcca                                                                                                                                                >  2:2022050/1‑148 (MQ=14)
        tCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTTAGAAATCCGGGGCTTTTGAAAGTGATGGTGGAGGGGGAAGGATTACTCAGCGCAGCGCGCATCGCCCTACGGGTCGTAGACAGAGTCAACGCACTCAAGCTGACGcaagag                                                                                                                                               >  1:1655357/1‑145 (MQ=14)
        tCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCTCTATCGCTGACGCTCgcg                                                                                                                                               >  2:2817770/1‑148 (MQ=255)
         ccccGCTTGGGCAGACTATTTGATTGCGGATTCGTTTGAGTATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGc                                                                                                                                              <  1:4823010/150‑1 (MQ=255)
         ccccGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGc                                                                                                                                              <  1:2782397/150‑1 (MQ=255)
          cccGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACGCTCTCACGCTGACGCTCGAGcc                                                                                                                                             >  1:5538985/1‑150 (MQ=255)
          cccGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCGGTTGAAAGTGATGGGGGTGGGGGACGGATTACTCCCCGCGGCGCCCATCGCCCGACGGGTCGAGCCCTCCGGCCACCCGCTCTCGCGCACGCGCGAGcc                                                                                                                                             >  2:6432632/1‑150 (MQ=9)
             gCTTGGGCAGAATCTTTGCTGGCGGATTCGGGTGAGAAGTCCGGGGCTTTTGCCCGGGCTGGTGGTGGGGGACGGATTACTCCGCGCGGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGaa                                                                                                                                          <  2:4971272/150‑1 (MQ=14)
                  ggCATAATATTTGATTGCGGATTAGTTTGAGAATTCCGGTGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTa                                                                                                                                     <  2:363785/150‑1 (MQ=255)
                   gCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAg                                                                                                                                    <  2:3880914/150‑1 (MQ=255)
                   gCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTAGCCCTTCGGGTCGTGGGCAGCGGCCACGCGCACTCCCGGGCGCGCGAGGCGTAccgtgg                                                                                                                                    >  1:5567920/1‑146 (MQ=11)
                   gCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGAGCTGAGAGCTTCGCACATCGGGTCGTTGCCTTCGGCAAAGCGCTATCGCTGACGCTAAAGACGCACCtttc                                                                                                                                    >  2:3958707/1‑148 (MQ=14)
                   gCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTAACAGCTGCGAGCGTCGCCCTACCGGAAGTTGACCGCTGCAAAGCACTCACGGGGTAGCGCGAGGCGTACCTTAg                                                                                                                                    >  1:129949/1‑150 (MQ=9)
                     aGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACGCTCTCTCGCAGACGCTCGAGCCGAACCATAGTc                                                                                                                                  >  2:2435210/1‑150 (MQ=255)
                     aGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCAGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACCCACACACGCAGACGCACCAGCCCAACCCTaata                                                                                                                                  >  1:1188697/1‑147 (MQ=14)
                     aGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGCGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTc                                                                                                                                  >  2:4085552/1‑150 (MQ=255)
                        atatTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGaa                                                                                                                               <  1:4402146/150‑1 (MQ=255)
                               agtGCGGATTCGTGTGAGCAGTCCGGGGCGTTTGAACGTGATGGTGGGGGGGGACGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTGGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTGAGTCGAAGCTtctc                                                                                                                        <  1:4277665/148‑1 (MQ=14)
                                   cGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATcc                                                                                                                    <  2:6451900/150‑1 (MQ=255)
                                   cGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGAGGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCGTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATcc                                                                                                                    <  2:1230602/150‑1 (MQ=255)
                                      agtCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTc                                                                                                                 <  2:4151270/148‑1 (MQ=255)
                                        tCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGACGAACCTTAGTCGAAGCTTCTCATCCATccc                                                                                                               >  2:1897809/1‑150 (MQ=255)
                                        tCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTccc                                                                                                               <  1:748312/150‑1 (MQ=255)
                                        gcGTTTGAGAATTCCGGGGCTTCTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTccc                                                                                                               <  2:4247071/149‑1 (MQ=255)
                                           tttGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACCCTCTCTCGCTGACGCACGAGCCGAACCATTGTCGAAGCTTCTAATCCTTCCCCGc                                                                                                            >  1:3960360/1‑150 (MQ=14)
                                              gagaATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACCCACTCACGCAGACGCACCAGCCGAACCGAAGTCGAAGCTTAACAAACATCCCCGCaag                                                                                                         >  1:1900104/1‑148 (MQ=14)
                                              gagaATTCCGGGGCAATAGAAACTGAAGCTGGTGGCGGAAGGATTACTCAGCGCTGCGAGCTTCGCCCTACGGCTATTTTCCTGCCCCAAATCACTCTCTCAGACCATCCAGCCGAACCTTACTCGAAGCTTCTCATCCTTCCCCGCATg                                                                                                         <  2:5012497/150‑1 (MQ=9)
                                               agtaTTGCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATgg                                                                                                        <  1:674185/147‑1 (MQ=255)
                                               agaATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCCACGCACACGCGCTGACGCTCGAGCCGAACCGTAGTCGAAGCTTCTCATACATCCCCCCCTgg                                                                                                        >  1:4677235/1‑150 (MQ=14)
                                               agaATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATgg                                                                                                        <  2:5347678/150‑1 (MQ=255)
                                                     ggttgtgttttGTAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGaa                                                                                                  <  1:4142683/143‑1 (MQ=255)
                                                         ggCTTTTGTAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAtt                                                                                              <  2:3411443/150‑1 (MQ=38)
                                                            ttttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGa                                                                                           <  1:4599327/150‑1 (MQ=38)
                                                               tGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCAGCGGCAACGCACTCTCGCTGACGCACGAGCCGAAACATAGTAGAAGCTTCTAATCCTTCCCCGCATGGGCAGCATAATTTATTg                                                                                        >  1:6464305/1‑150 (MQ=11)
                                                                 aaaGTGATGGTGGTGGGGGAAGGAATACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAAAGCTCACACGCTGACGCTCGGGCCGAACCTTAGTCGAAGCGTCTCATCCTTCCCCGCCGGGGCAGGATAtttgtttggg                                                                                      >  1:3966998/1‑148 (MQ=11)
                                                                     tGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACCCTCACTCGCTGACGCGCGAGCCGAACCGTAGTCCAAGCGTCACATCCGTCCCCCCAAGGGCAGAATATTTTATTGCGGAtt                                                                                  >  2:5543237/1‑150 (MQ=11)
                                                                            ggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTga                                                                           <  1:3987461/150‑1 (MQ=255)
                                                                              tggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTgaga                                                                         <  2:5335422/150‑1 (MQ=255)
                                                                              tggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTgaga                                                                         <  1:625988/150‑1 (MQ=255)
                                                                              tggGGGAAGGATGACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTgaga                                                                         <  1:155628/150‑1 (MQ=255)
                                                                               gggggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGagag                                                                        <  2:1123919/150‑1 (MQ=255)
                                                                                 gggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGtt                                                                      <  2:1655357/150‑1 (MQ=255)
                                                                                  ggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCATAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCACAAAATTTGATTGCGGATTCGGTTGGGAGTTc                                                                     >  2:3153541/1‑150 (MQ=14)
                                                                                    aaGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAg                                                                   <  2:3516053/150‑1 (MQ=255)
                                                                                    aaGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAg                                                                   <  1:2469983/150‑1 (MQ=255)
                                                                                    aaGGATTACTCAGCGCTGCGCGCTTCGCCCGTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAg                                                                   <  1:2355337/150‑1 (MQ=255)
                                                                                     aGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAgg                                                                  <  2:2915655/150‑1 (MQ=255)
                                                                                     aGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAgg                                                                  <  1:1446620/150‑1 (MQ=255)
                                                                                      ggATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAggg                                                                 <  1:6561193/150‑1 (MQ=255)
                                                                                         ttACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTTATTGCGGATTCGGTGGGGAGTGCAGGGACt                                                              >  2:3586482/1‑150 (MQ=11)
                                                                                         ttACTCAGCGCCGCGCTCTTCGCCCTTCTGTTCGTTGCCGGCGGCATCGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACt                                                              <  2:4077705/150‑1 (MQ=17)
                                                                                          tACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACtt                                                             <  2:1558218/150‑1 (MQ=255)
                                                                                          tACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACtt                                                             <  1:5414158/150‑1 (MQ=255)
                                                                                           aCTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTGGTGAGTTCAGGGACgtt                                                            >  2:5936082/1‑150 (MQ=255)
                                                                                           aCTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGCTTGCGGATTCGGTGGGGAGGTCAGGGACgtt                                                            >  2:5935886/1‑150 (MQ=11)
                                                                                            cTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCTAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACtttt                                                           <  1:3086844/150‑1 (MQ=255)
                                                                                            cTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACtttt                                                           <  1:5895884/150‑1 (MQ=255)
                                                                                             tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCACGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCAACCGTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCGTTAGGGGTAGGGGCCTTTTg                                                          >  2:2975620/1‑150 (MQ=14)
                                                                                                    tGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGa                                                   <  1:2435210/150‑1 (MQ=255)
                                                                                                      cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTTATTGCGGATTCGGTTGGGAGTTCAGGGACTTTTGAAAGTGAtg                                                 >  2:3690334/1‑150 (MQ=255)
                                                                                                      cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtg                                                 <  1:1160387/150‑1 (MQ=255)
                                                                                                        cgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCg                                                                                                                                                 >  2:6284744/1‑52 (MQ=35)
                                                                                                           ttCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggtggt                                            <  1:5304894/150‑1 (MQ=255)
                                                                                                            tCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggtggtg                                           <  1:3055892/150‑1 (MQ=255)
                                                                                                            tCGCCCATCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCAAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggtggtg                                           <  1:4241142/150‑1 (MQ=11)
                                                                                                              gCCCTTCGGGTCGTTGCCTTCGGCAACGCTCTCTCGCTGCCGCTAGAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATTGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGATAGTGATGGTGGTggg                                         <  1:3251979/150‑1 (MQ=14)
                                                                                                              gCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGATTATTTGATTGCGGATTCGCTTGAGAGTTCAGGTACTTTTGAAAGTGATGGTGGTggg                                         <  2:1276244/150‑1 (MQ=255)
                                                                                                              gCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTTATTGCGGATTAGTTTGTGAGTTCAGGGACTTTTGAAAGTGATGGTGGTggg                                         >  2:4797033/1‑150 (MQ=255)
                                                                                                               cgcTTCGGGTCGTTGCGTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTgggg                                        <  1:2022050/148‑1 (MQ=255)
                                                                                                               cccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTgggg                                        <  1:745349/150‑1 (MQ=255)
                                                                                                                ccgtCGGGGCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTggggg                                       <  1:3471866/147‑1 (MQ=255)
                                                                                                                ccTTCGGGGCGTTGGGTGCGGCTACGCTCTCTCGCGGTCGCTCGAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCCGGGTCTTTTGAAAGTGATGGTGGTggggg                                       <  1:3470791/150‑1 (MQ=11)
                                                                                                                 cTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGTTTGCGGAGTCGCTGGGGAGTTCAGGGACTTTTGAAAGGGATGGTGGGGGGGGa                                      >  2:3762132/1‑150 (MQ=14)
                                                                                                                  ttCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCTAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGaa                                     <  1:6394488/150‑1 (MQ=255)
                                                                                                                  ttCGGGTCGTTGCCTGCGGCAAAGCTCTCACGCTGACGCTCGAGCCGAACCATAGTCGAAGCTTCTCATCCTTCCCCGCATTGGCAGAATATTTTATTGCGGATTAGGTGGTGAGTTCAGGGACTTTTGAAAATTAAGGGGGTggggcgt                                     >  2:3350529/1‑147 (MQ=11)
                                                                                                                   tCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCTAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAg                                    <  1:856251/150‑1 (MQ=255)
                                                                                                                    gagagcGTTGCCGCAGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCCGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGTAAGGGTTGGTGGTGGGGGACgg                                   >  1:4447824/6‑150 (MQ=17)
                                                                                                                      ggTCGTTGCCTGCGGCAAAGCTCTAACCCTGAAGCTCGAGCCGAACCTAAGTCGAAGCGTCTCAAACTTACCCCCATGGGCAGAAAAATTTATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTTGTGGGGGGAGGAt                                 >  2:739661/1‑150 (MQ=14)
                                                                                                                           ttGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTc                            >  2:588970/1‑150 (MQ=255)
                                                                                                                            tGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCa                           <  2:2215543/150‑1 (MQ=255)
                                                                                                                              ccTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGTGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCCTAAga                         >  2:2359835/1‑149 (MQ=255)
                                                                                                                               cTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgcg                        >  2:5439298/1‑150 (MQ=255)
                                                                                                                                  cGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTg                     <  2:5490154/150‑1 (MQ=255)
                                                                                                                                       aCGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTTATGGTGGTGGGGGAAGGATTACTCAGCGCTgcgcgc                >  1:5903377/1‑150 (MQ=255)
                                                                                                                                       aCGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCATCCCCGCATGGGCACAAGATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTAATCAGCGCGgcgcgc                >  1:1725950/1‑150 (MQ=17)
                                                                                                                                       aCGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTATCATCCATCCCCCCATGGGCAGAATATTTTATTTCGGATTAGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACAAAGCGCTgcgcgc                >  1:257565/1‑150 (MQ=17)
                                                                                                                                          ctctctCGCTGACGCTCGAGCCGAACCTTCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTc             <  2:5633202/150‑1 (MQ=255)
                                                                                                                                           tctctcGCTGTCGCTCTAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGTTTATTTGATTGCGGATTCGCTTGAGTTTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCg            <  2:2666503/150‑1 (MQ=14)
                                                                                                                                           tctctcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCTCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCg            <  2:4014458/150‑1 (MQ=255)
                                                                                                                                           tctctcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCg            <  2:1424739/150‑1 (MQ=255)
                                                                                                                                           tctctcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTGGAGAGTTCAGGGACTTTGGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCGGCGCGCTTCg            >  2:6406881/1‑150 (MQ=255)
                                                                                                                                           tctctcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTACCCGCATGGGCAGAATATTTGATTGCGGATTAGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGAGCGCTTCg            >  1:4370279/1‑150 (MQ=255)
                                                                                                                                           tctcacGCTGCCGCTCTAGACGAACCGTCGTCGAAGCTTCTCATCCGTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTTAGATTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCg            <  2:4483034/150‑1 (MQ=14)
                                                                                                                                             tctcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGATTTTTTGATTGCGGATTCGATTGAGAGTTCAGGGACTTATGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGcc          <  2:6504186/150‑1 (MQ=255)
                                                                                                                                                    gACGCTCGAGCCGAACCTTCGTCGACGCTTCTCCTCCGTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCggg   <  1:3352435/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGT  >  NC_000913/1286816‑1287114

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: