Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,286,955 T→A 6.8% intergenic (+429/+111) 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,9550TA 6.8% 560.1 / 5.6 199intergenic (+429/+111)narI/rttRnitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript
Reads supporting (aligned to +/- strand):  ref base T (53/122);  new base A (15/2);  total (75/124)
Fisher's exact test for biased strand distribution p-value = 4.02e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 222 total aligned reads displayed.

GTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCG  >  NC_000913/1286806‑1287104
                                                                                                                                                      |                                                                                                                                                     
gTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGct                                                                                                                                                       >  1:5335422/1‑150 (MQ=255)
gTCGAAGCTTCTCATCCTTCCCCGATTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGct                                                                                                                                                       <  2:6126774/150‑1 (MQ=255)
   gAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGctcgc                                                                                                                                                    >  1:1389341/1‑148 (MQ=255)
    aaGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGctctct                                                                                                                                                   <  2:1314379/150‑1 (MQ=255)
    aaGCTTCTCATCCTTCCCCGCGTGGGCAGAATATTTTATTTGGGATTCGTTTGAGATTTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGctctct                                                                                                                                                   <  2:4633142/150‑1 (MQ=11)
    aaGCTTCTCATCCGTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGctctct                                                                                                                                                   <  2:5879351/150‑1 (MQ=255)
     aGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCtctctc                                                                                                                                                  <  1:1731187/150‑1 (MQ=255)
     aGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCtctctc                                                                                                                                                  <  1:254349/150‑1 (MQ=255)
          ctccgccTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTATGAGAATTCCGGGGCTTTTGCAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGa                                                                                                                                             <  1:3521872/145‑1 (MQ=255)
          ctcATCCTTCCCCGCTTGGGCAGAATAATTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGa                                                                                                                                             <  2:2973262/150‑1 (MQ=255)
             aTCCTTCCCCGCTTGGGCATATTATTTGATTTCGTATTCGTTTGAGAATTCCGGTGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGc                                                                                                                                          <  2:5242031/150‑1 (MQ=14)
             aTCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGGTGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACCCACTAACGCAGACGc                                                                                                                                          >  1:5178587/1‑150 (MQ=14)
                 ttCCGCGCTTGGGCAGTATATTTGTTTGCGTATTCGTTTGAGTTTTACGTGTATTTTGCAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGa                                                                                                                                      <  2:1237447/150‑1 (MQ=14)
                 ttCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGa                                                                                                                                      >  2:6078444/1‑150 (MQ=255)
                 ttCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTAGCAAGCGGCAACGCTCACACGCTGACGCTcca                                                                                                                                      >  1:4945414/1‑148 (MQ=14)
                 gtCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGa                                                                                                                                      <  2:5596869/149‑1 (MQ=255)
                  tCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTTGGGGAAGGATTACTCAGCGCAGCGCGCTTCGCCCTTCTGGTC‑CTTTCCTGCGGCAAAGCCCTCTCGCAGACGCTCGac                                                                                                                                     >  1:1655361/1‑149 (MQ=14)
                   ccccGCTTGGGCAGACTATTTGATTGCGGATTCGTTTGAGTATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGc                                                                                                                                    <  1:4823010/150‑1 (MQ=255)
                       gCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGaa                                                                                                                                <  2:4854323/150‑1 (MQ=255)
                             gCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAg                                                                                                                          <  2:3880914/150‑1 (MQ=255)
                             gCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTAGCCCTTCGGGTC‑GTGGGCAGCGGCCACGCGCACTCCCGGGCGCGCGAGGCGTAccgtgg                                                                                                                          >  1:5567920/1‑146 (MQ=11)
                               ataaTATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTAGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCAGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCCAAACATAGTc                                                                                                                        >  2:3580700/3‑150 (MQ=14)
                               aGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGCGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTc                                                                                                                        >  2:4085552/1‑150 (MQ=255)
                                  atatTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGaa                                                                                                                     <  1:4402146/150‑1 (MQ=255)
                                   tatTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAg                                                                                                                    <  1:3340524/150‑1 (MQ=255)
                                         agtGCGGATTCGTGTGAGCAGTCCGGGGCGTTTGAACGTGATGGTGGGGGGGGACGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTGGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTGAGTCGAAGCTtctc                                                                                                              <  1:4277665/148‑1 (MQ=14)
                                         aTTGCGGATTCGTTTGAGCATTCCGGGGCTTTTGAAAGAGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTtctc                                                                                                              <  2:3744110/150‑1 (MQ=255)
                                                agtCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTc                                                                                                       <  2:4151270/148‑1 (MQ=255)
                                                  tCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCAGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTacc                                                                                                     >  2:1986129/1‑150 (MQ=255)
                                                  gcGTTTGAGAATTCCGGGGCTTCTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTccc                                                                                                     <  2:4247071/149‑1 (MQ=255)
                                                        gagaATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATg                                                                                               <  2:223403/150‑1 (MQ=255)
                                                        gagaATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACCCACTCACGCAGACGCACCAGCCGAACCGAAGTCGAAGCTTAACAAACATCCCCGCaag                                                                                               >  1:1900104/1‑148 (MQ=14)
                                                        gagaATTCCGGGGCAATAGAAACTGAAGCTGGTGGCGGAAGGATTACTCAGCGCTGCGAGCTTCGCCCTACGGCTA‑TTTTCCTGCCCCAAATCACTCTCTCAGACCATCCAGCCGAACCTTACTCGAAGCTTCTCATCCTTCCCCGCATg                                                                                               <  2:5012497/150‑1 (MQ=9)
                                                         agtaTTGCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATgg                                                                                              <  1:674185/147‑1 (MQ=255)
                                                         agaATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCCACGCACACGCGCTGACGCTCGAGCCGAACCGTAGTCGAAGCTTCTCATACATCCCCCCCTgg                                                                                              >  1:4677235/1‑150 (MQ=14)
                                                         agaATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATgg                                                                                              <  2:5347678/150‑1 (MQ=255)
                                                               ggttgtgttttGTAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGaa                                                                                        <  1:4142683/143‑1 (MQ=255)
                                                                   ggCTTTTGTAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAtt                                                                                    <  2:3411443/150‑1 (MQ=38)
                                                                    gCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAttt                                                                                   <  2:2838541/150‑1 (MQ=38)
                                                                    gCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCGCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAttt                                                                                   >  2:5634422/1‑150 (MQ=38)
                                                                      ttttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGa                                                                                 <  1:4599503/150‑1 (MQ=38)
                                                                         tGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTTCCAGCGGCAACGCACTCTCGCTGACGCACGAGCCGAAACATAGTAGAAGCTTCTAATCCTTCCCCGCATGGGCAGCATAATTTATTg                                                                              >  1:6464305/1‑150 (MQ=11)
                                                                           aaaGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACCCACACTCGCAGACGCTCGAGCCGAACCGTAGTCGAAGCTTCTCAACCATCCCCGCATGGGCCGCATATATGATTGCg                                                                            >  2:2612871/1‑150 (MQ=14)
                                                                            aaGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCgg                                                                           <  2:3591090/150‑1 (MQ=255)
                                                                                  tggtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGc                                                                     <  1:747472/150‑1 (MQ=40)
                                                                                  gttcgttggGGGACTGATTACTCTGCGCTGCCCG‑GTCGGCCGTCCGGTCTGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCTAGCCGAACCTTCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGc                                                                     <  2:6466024/144‑1 (MQ=11)
                                                                                      ggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTga                                                                 <  1:5969206/150‑1 (MQ=255)
                                                                                      ggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTga                                                                 <  1:924733/150‑1 (MQ=255)
                                                                                      ggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTga                                                                 <  1:3987442/150‑1 (MQ=255)
                                                                                       gtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTgag                                                                <  1:5488705/150‑1 (MQ=255)
                                                                                        tggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTgaga                                                               <  2:5335422/150‑1 (MQ=255)
                                                                                        tggGGGAAGGATGACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTgaga                                                               <  1:155628/150‑1 (MQ=255)
                                                                                         gggggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGTATATTTGATTGCGGATTCGCTTGagag                                                              <  1:2768332/150‑1 (MQ=255)
                                                                                           gggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGtt                                                            <  2:1655357/150‑1 (MQ=255)
                                                                                              aaGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCTAGTCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAg                                                         <  1:755363/150‑1 (MQ=255)
                                                                                              aaGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAg                                                         <  1:2469983/150‑1 (MQ=255)
                                                                                              aaGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAg                                                         <  2:4493356/150‑1 (MQ=255)
                                                                                              aaGGATTACTCAGCGCTGCGCGCTTCGCCCGTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAg                                                         <  1:2355337/150‑1 (MQ=255)
                                                                                               aGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAgg                                                        <  1:1446620/150‑1 (MQ=255)
                                                                                                ggATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAggg                                                       <  1:6561193/150‑1 (MQ=255)
                                                                                                 gATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGa                                                      <  2:3835547/150‑1 (MQ=255)
                                                                                                   ttACTCAGCGCCGCGCTCTTCGCCCTTCTGTTC‑GTTGCCGGCGGCATCGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACt                                                    <  2:4077705/150‑1 (MQ=17)
                                                                                                    tACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCTAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACtt                                                   <  1:4560607/150‑1 (MQ=255)
                                                                                                    tACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACtt                                                   <  2:1558218/150‑1 (MQ=255)
                                                                                                    tACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCGTAGTCGAAGCTTCTCATACTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACtt                                                   >  2:208577/1‑150 (MQ=255)
                                                                                                     aCTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTGGTGAGTTCAGGGACgtt                                                  >  2:5936082/1‑150 (MQ=255)
                                                                                                     aCTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGCTTGCGGATTCGGTGGGGAGGTCAGGGACgtt                                                  >  2:5935886/1‑150 (MQ=11)
                                                                                                     aCTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACttt                                                  <  2:2411459/150‑1 (MQ=255)
                                                                                                      cTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCTAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACtttt                                                 <  1:3086844/150‑1 (MQ=255)
                                                                                                      cTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACtttt                                                 <  1:5895884/150‑1 (MQ=255)
                                                                                                       tCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCATCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGGGAGTTCAGGGACTTTTg                                                >  2:3647735/1‑150 (MQ=255)
                                                                                                       tCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCACGCTGACGCTCGAGCCGAACCTTAGTCGAAGCGTCTCATCCTTCCCCGCATGGGCAGAATATTTTATTGCGGATTCGCTGGTGAGTTCAGGGACTTTTg                                                >  2:641913/1‑150 (MQ=14)
                                                                                                             cTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAATTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTg                                          >  2:3947015/1‑150 (MQ=255)
                                                                                                              tGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGa                                         <  1:2435210/150‑1 (MQ=255)
                                                                                                                gtcgcTTCGCGCTTCTGTTC‑GTTTCTTGCGGCAACGCTCTCTCGCTGTCGCTCTAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtg                                       <  2:6436781/148‑1 (MQ=14)
                                                                                                                cgcgcTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtg                                       <  1:1160387/150‑1 (MQ=255)
                                                                                                                  cgcTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCg                                                                                                                                       >  2:6284744/1‑52 (MQ=35)
                                                                                                                  cgcTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCg                                                                                                                                       <  1:6284744/52‑1 (MQ=35)
                                                                                                                     ttCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggtggt                                  <  1:5304894/150‑1 (MQ=255)
                                                                                                                      tCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggtggtg                                 <  1:3055892/150‑1 (MQ=255)
                                                                                                                      tCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTATTGAAAGTGAtggtggtg                                 <  1:3964062/150‑1 (MQ=255)
                                                                                                                      tCGCCCATCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCAAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggtggtg                                 <  1:4241142/150‑1 (MQ=11)
                                                                                                                        gCCCTTCGGGTC‑GTTGCCTTCGGCAACGCTCTCTCGCTGCCGCTAGAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATTGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGATAGTGATGGTGGTggg                               <  1:3251979/150‑1 (MQ=14)
                                                                                                                        gCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTTATTGCGGATTAGTTTGTGAGTTCAGGGACTTTTGAAAGTGATGGTGGTggg                               >  2:4797033/1‑150 (MQ=255)
                                                                                                                        gCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTggg                               <  2:1188697/150‑1 (MQ=255)
                                                                                                                         cgcTTCGGGTC‑GTTGCGTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTgggg                              <  1:2022050/148‑1 (MQ=255)
                                                                                                                          ccTTCGGGGC‑GTTGGGTGCGGCTACGCTCTCTCGCGGTCGCTCGAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCCGGGTCTTTTGAAAGTGATGGTGGTggggg                             <  1:3470791/150‑1 (MQ=11)
                                                                                                                            ttCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCTAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGaa                           <  1:6394488/150‑1 (MQ=255)
                                                                                                                            ttCGGGTC‑GTTGCCTGCGGCAAAGCTCTCACGCTGACGCTCGAGCCGAACCATAGTCGAAGCTTCTCATCCTTCCCCGCATTGGCAGAATATTTTATTGCGGATTAGGTGGTGAGTTCAGGGACTTTTGAAAATTAAGGGGGTggggcgt                           >  2:3350529/1‑147 (MQ=11)
                                                                                                                                ggTC‑GTTGCCTGCGGCAAAGCTCTAACCCTGAAGCTCGAGCCGAACCTAAGTCGAAGCGTCTCAAACTTACCCCCATGGGCAGAAAAATTTATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTTGTGGGGGGAGGAt                       >  2:739661/1‑150 (MQ=14)
                                                                                                                                  tC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTa                     <  1:1197580/150‑1 (MQ=255)
                                                                                                                                          cTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgcg              >  2:5439298/1‑150 (MQ=255)
                                                                                                                                             cGGCAACGCTCTCTCGCTGCCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTg           <  2:5490201/150‑1 (MQ=255)
                                                                                                                                               gCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGTGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTgcg         <  2:2740055/150‑1 (MQ=255)
                                                                                                                                                  aCGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTTATGGTGGTGGGGGAAGGATTACTCAGCGCTgcgcgc      >  1:5903377/1‑150 (MQ=255)
                                                                                                                                                     ctctctCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAGTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTc   <  1:363194/150‑1 (MQ=255)
                                                                                                                                                      tctctcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCTCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCg  <  2:4014458/150‑1 (MQ=255)
                                                                                                                                                      tctctcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCg  <  2:1424739/150‑1 (MQ=255)
                                                                                                                                                      tctctcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCGTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCg  >  2:3133629/1‑150 (MQ=255)
                                                                                                                                                      tctctcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTACCCGCATGGGCAGAATATTTGATTGCGGATTAGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGAGCGCTTCg  >  1:4370279/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTC‑GTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCG  >  NC_000913/1286806‑1287104

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: ‑