Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,286,983 A→T 7.9% intergenic (+457/+83) 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,9830AT 7.9% 755.6 / 2.8 284intergenic (+457/+83)narI/rttRnitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript
Reads supporting (aligned to +/- strand):  ref base A (100/141);  new base T (8/17);  total (120/164)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 288 total aligned reads displayed.

GGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG  >  NC_000913/1286834‑1287131
                                                                                                                                                     |                                                                                                                                                    
ggCAGAATATTTTATTGCGGCTTCGTTTGAGAATTCCGGGGCTTTAGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCAGACGCTCGAGCCGAACCTTa                                                                                                                                                      <  1:3559597/150‑1 (MQ=255)
  cAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGTGGCTTTTGATAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGt                                                                                                                                                    <  1:6311544/150‑1 (MQ=255)
     aaTATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAACGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGa                                                                                                                                                 <  2:10293843/150‑1 (MQ=255)
          ttGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCtt                                                                                                                                            <  1:388104/150‑1 (MQ=255)
           tGATTGCGGATTCGTTTGAGAATTCGGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTtc                                                                                                                                           <  1:7480432/150‑1 (MQ=255)
              ttGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCa                                                                                                                                        >  2:6372787/1‑150 (MQ=255)
              ttGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCCACGCACTCACGCTGACGCTCGAGCCGACACAATGTCGAAGCTTATCa                                                                                                                                        >  1:10104955/1‑150 (MQ=14)
               tGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCAt                                                                                                                                       <  1:4337398/150‑1 (MQ=255)
                     ttCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGACCCTTAGTCGAAGCTTCTCatccatcc                                                                                                                                 >  1:6612386/1‑150 (MQ=255)
                     ttCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGCGGGGGAAGGATTACTCAGCGCTGCGCGATTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCACTCGCGCTGGCGCACGAGCCGTACCTTGGTCGAAGCTTCTCATCCTTcc                                                                                                                                 >  2:2440740/1‑150 (MQ=11)
                     ggcGTTTGAGTATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTcc                                                                                                                                 <  2:6727896/148‑1 (MQ=255)
                         tttGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACCCTCACACGCAGACCCGCCACCCCAACCATAGTCGAAGCATCTCAACCATaccccc                                                                                                                             >  1:8195498/1‑148 (MQ=14)
                          ttGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCa                                                                                                                            >  2:5393357/1‑150 (MQ=255)
                              gaATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATggg                                                                                                                        <  1:9242238/150‑1 (MQ=255)
                               aaTTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTTGTCGAAGCTTCTCATCCATCCCCGCATGGGc                                                                                                                       >  2:5471297/1‑150 (MQ=255)
                                 ttCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAg                                                                                                                     <  1:1918217/150‑1 (MQ=255)
                                    cGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCCCATCGCCCATCGGGTCGATGCCGGCGGCCACGCGCTCACGCTGACGCACGAGCCGAACCGTCGTAGAAGCTACACATCCTTCCCCGCAGAGGCAGAat                                                                                                                  >  2:5646374/1‑150 (MQ=11)
                                       gtcTTTTGTAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAtt                                                                                                               <  2:8275219/148‑1 (MQ=37)
                                           tttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGCCGCTAGAGCCGAACCGTTGTCGAAGCTTAAAATCCATCCCCGCAACGGCACAATAattaaa                                                                                                           >  2:159253/1‑147 (MQ=11)
                                             tGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTGGCCTGCCGCAACGCTCTCTCGCGGACGCGCGAGCCGAACCGTAGTCGAAGCATCTCAACCTACCCCGCAAGGGCAGGATATATGaatt                                                                                                         >  1:3103772/1‑147 (MQ=11)
                                              gAAAGTGATGGTGGTGGGGGAAGGATTACTAAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTCACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTACCCGCAAGGGCAGAATATTTTAATGc                                                                                                        >  2:9228694/1‑150 (MQ=11)
                                                 aGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGa                                                                                                     <  2:7148568/150‑1 (MQ=39)
                                                 aGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGa                                                                                                     <  2:5415260/150‑1 (MQ=39)
                                                      tggtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCCCTGACGCACGAGCCGAACCTTAGTCGAAGCATATCATCCTTCCCCGCATGTGAAGAATATTTGATTTCGGAttagg                                                                                                >  2:6062828/1‑147 (MQ=14)
                                                        gtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCtt                                                                                              <  1:7725594/150‑1 (MQ=255)
                                                        gtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCGTAGTCGAAGCTTCTCATCCGTCCCCGCATGGGCAGACTATTTGATTGCGGATTCGCtt                                                                                              >  2:7450968/1‑150 (MQ=255)
                                                         tggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTg                                                                                             >  2:6824610/1‑150 (MQ=38)
                                                           gtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGc                                                                                                <  1:9604917/145‑1 (MQ=40)
                                                              ggggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGATTGTGAGt                                                                                        >  2:8510144/1‑150 (MQ=255)
                                                               gggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGtt                                                                                       <  2:1222141/150‑1 (MQ=255)
                                                               gggAAGGATTACTCAGCGCTGCGCGCGTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGtt                                                                                       <  1:5411599/150‑1 (MQ=255)
                                                                ggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCGTGGGCGAAGCTTCGCATCCTTCCCCGCATGGGCCGCCTCTTTGATTGCGGATTCGCTTGGGAGTTc                                                                                      >  2:8145294/1‑150 (MQ=11)
                                                                  aaGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAg                                                                                    <  1:6824610/150‑1 (MQ=255)
                                                                  aaGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAg                                                                                    <  1:6118664/150‑1 (MQ=255)
                                                                   aGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAgg                                                                                   <  1:3192409/150‑1 (MQ=255)
                                                                   aGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAgg                                                                                   <  1:8561756/150‑1 (MQ=255)
                                                                        tACTCAGCGCTGCGCGCGTCGCCCGCCGGGTCGTTGCCTGCGGCAACGCCCTCTCTCTGACGCTCGAGACCAACCTTAGTCGAAGCTTCTCATCCTTCCCCTCATGGGCATAATTTTTGATTGCGGATTCGCTTGATAGTTCAGGGACtt                                                                              <  1:5689603/150‑1 (MQ=14)
                                                                         aCTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACttt                                                                             <  2:3121315/150‑1 (MQ=255)
                                                                          cTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACtttt                                                                            <  1:9376776/150‑1 (MQ=255)
                                                                          cTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCGCTCTCGCTGTCGCTCGAGCCGAACCGTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGGGGGTGAGTTCAGGGACtttt                                                                            >  2:3048701/1‑150 (MQ=11)
                                                                           tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCATAGTCGAAGCTTCTCAACCTTGCCCGCATGGGCAGAATATTTGGTTGCGGATTCGCTTGGGAGTTCAGGGACGTTTg                                                                           >  2:9614124/1‑150 (MQ=14)
                                                                             aGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGaa                                                                         <  1:1463760/150‑1 (MQ=255)
                                                                             aGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGaa                                                                         <  1:2627674/150‑1 (MQ=255)
                                                                             aGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCTCGCTCGCTGACGCGCGAGCCGAACCGTAGTCGAAGCGTCTCATCCTTCCCCGCCGGGGCAGAATATTTGATTGCGGCTTCGCGTGAGAGTTCAGGGCCGTTTGaa                                                                         >  2:4729795/1‑150 (MQ=14)
                                                                              gcgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGaaa                                                                        <  1:2137471/150‑1 (MQ=255)
                                                                               cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAg                                                                       <  2:4418967/150‑1 (MQ=255)
                                                                               cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAg                                                                       <  2:9085336/150‑1 (MQ=255)
                                                                               cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAg                                                                       <  2:1320037/150‑1 (MQ=255)
                                                                                  tGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCGTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCCGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACGTTTGAAAGTGa                                                                    >  2:7178844/1‑150 (MQ=255)
                                                                                   gcgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAt                                                                   <  1:2816941/150‑1 (MQ=255)
                                                                                    cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtg                                                                  <  1:10276371/150‑1 (MQ=255)
                                                                                          tCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggtggtg                                                            <  1:7393194/150‑1 (MQ=255)
                                                                                           cGCCCTTCGGGTCGGTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggtggtgg                                                           <  1:966491/150‑1 (MQ=255)
                                                                                            gCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCGTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTggg                                                          <  2:4772069/150‑1 (MQ=255)
                                                                                            gCCCGTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTggg                                                          <  1:2059799/150‑1 (MQ=255)
                                                                                             cccTTCGGGTCGTTGCCGGCGGCAACGCTCTCTCGCGGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCAAACTTACCCGCATGGGCAGAAAATTTGATTGCGGTTTCGCTTGTGAGGTCAGGGACTTTTGAAAGTGATGGTGGTgggg                                                         >  2:3747582/1‑150 (MQ=14)
                                                                                                ttCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGaa                                                      <  2:1790285/150‑1 (MQ=255)
                                                                                                 tCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGACGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGTGTTCATGTACTTTTGAAAGCGATGGTGGTGGGGGAAg                                                     <  1:27254/150‑1 (MQ=14)
                                                                                                  cGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAgg                                                    <  1:3417559/150‑1 (MQ=255)
                                                                                                  cGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAgg                                                    <  1:10027135/150‑1 (MQ=255)
                                                                                                   gggTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATTGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGa                                                   <  1:4439272/150‑1 (MQ=255)
                                                                                                   gggTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATTTTTGATTGCGGATTCGATTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAgta                                                   <  1:3744955/150‑3 (MQ=255)
                                                                                                    ggTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAt                                                  <  2:2030885/150‑1 (MQ=255)
                                                                                                      tctttGCCTGCGGCAACGCTCTATCTCTGACGATCTAGCAGAACCTAACTAGATTATTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTa                                                <  1:1609849/147‑1 (MQ=17)
                                                                                                      tCGTTGCCTGCGGCAACGCTCTCTCGCTGACTCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTATTGAAAGTGATGGTGGTGGGGGAAGGATTa                                                <  1:7994966/150‑1 (MQ=255)
                                                                                                      tCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTa                                                >  2:7324754/1‑150 (MQ=255)
                                                                                                      tCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTa                                                <  1:7356854/150‑1 (MQ=255)
                                                                                                       cGTTGCCGCAGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAttt                                                                                                              <  2:9358004/87‑1 (MQ=14)
                                                                                                         ttGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGCAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCTGTTACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTc                                             <  1:182791/150‑1 (MQ=14)
                                                                                                            cgatctgCAACGCTCTCTCGCTGTCGCTCTAGTCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgc                                          <  2:8063695/144‑1 (MQ=14)
                                                                                                            ccgcaggCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGGTTGGGAGTTCAGGGACTTTTGAAAGGGTGGGTGGTGGGGCGAGGATTACACAgc                                          >  2:8837180/6‑150 (MQ=14)
                                                                                                            ccTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAAGGTGGTGGGGGAAGGATTACTCata                                          >  2:6560456/1‑148 (MQ=255)
                                                                                                               gCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCATCTCATCCATCCCCCCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAATTGATGGTGGTGGGGGAAGGATTACTAAGCGCt                                       >  1:7106387/1‑150 (MQ=17)
                                                                                                                  gCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTAGTCAGCGCTtcg                                    >  2:1531257/1‑150 (MQ=255)
                                                                                                                    aaCGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAAGATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTgcgcg                                  >  2:520463/1‑150 (MQ=255)
                                                                                                                     aCGCTCTCTCGCTGTCGCTCGAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTgcgcgc                                 <  2:7041116/150‑1 (MQ=255)
                                                                                                                       gCTCTCTCGCTGACGCGCGCGCCGAACCTTAGTCGAAGCTTCTCCCCCGTCCCCGCTTGGGCAGACTATTAGATTGCGGATTAGATTGAGTTTTCAGGGACTTTTGCAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCtt                               <  2:2597289/150‑1 (MQ=14)
                                                                                                                        ctgtgtcGCTTTCGCTCGAGTCGAACCTTAGTCTAAGCTTCTCATCCGTCCCCGCTTGGGCATATTTTTTGATTGCGGATTCGCTTGAGAGTTCCGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTc                              <  1:10205133/145‑1 (MQ=11)
                                                                                                                        ctctctCGCTGCCTCTCCAGCCGCACCTTCGTCGAAGGTTCTCATCCTTCCCCGCCTGGGCAGATTTTATGATTGCGGATTCGCTTGAGAGTTGTGTGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTc                              <  1:9928541/150‑1 (MQ=14)
                                                                                                                        ctctctCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTc                              <  2:5009188/150‑1 (MQ=255)
                                                                                                                        ctctctCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTATCATCCTTCCCCCCAAGGGCAGAATATTTGATTGCGGATTCGATTGTGAGTTCAGGGACTTTTGAAATTGATGGTGGTGGGGGAAGGATTACAAACCGCTGCGCGCtac                              >  1:5265813/1‑148 (MQ=17)
                                                                                                                         tctctcGCTGTCGCTCTAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGATTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCg                             <  1:5070130/150‑1 (MQ=255)
                                                                                                                         gctctcGCTGACGCTCGAGCCGATCCTTAGTCGAAGGTTCTCATCCTTCCCCGCATGGGCAGACTAATTGATTGCGGCTTCGATTGAGCGTTCAGGGACTTTTGTAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCg                             <  1:3255424/149‑1 (MQ=17)
                                                                                                                            ctcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGccc                          >  1:9179588/1‑150 (MQ=255)
                                                                                                                                cTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCg                      <  1:8119731/150‑1 (MQ=255)
                                                                                                                                cTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCg                      <  1:8119324/150‑1 (MQ=255)
                                                                                                                                cTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGAGCTGCGCGCTTCGCCCTTCg                      >  2:4962024/1‑150 (MQ=255)
                                                                                                                                  ctcGCTCTAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCggg                    <  2:4019419/148‑1 (MQ=255)
                                                                                                                                   aCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGt                   <  2:7848824/150‑1 (MQ=255)
                                                                                                                                       tCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGtt               >  1:8854177/1‑150 (MQ=255)
                                                                                                                                       tCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGtt               <  1:4729843/150‑1 (MQ=255)
                                                                                                                                        cgcgCCGAACCTTTGTCGTCGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTg              <  2:1755991/147‑1 (MQ=14)
                                                                                                                                        cGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTg              >  1:7899324/1‑150 (MQ=255)
                                                                                                                                        cGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTGGAGAGTTCAGGGACTTTTCAAACGGAGGGGCGTCGGGGAAGGAGAACTCAGCGCTGCGCGCTTCGCCCTGCGGGTCggag              >  2:9122080/1‑147 (MQ=17)
                                                                                                                                              gAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTTCGCGCTTCGCCCTTCGGGTCGTTGCCTGCg        >  2:97983/1‑150 (MQ=255)
                                                                                                                                               caCCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGAGTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCgg       <  2:8195498/149‑1 (MQ=255)
                                                                                                                                                  cTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAg                                                                                    >  1:5768820/1‑70 (MQ=255)
                                                                                                                                                  cTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAg                                                                                    <  2:5768820/70‑1 (MQ=255)
                                                                                                                                                  cTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGa                                                                    >  2:5741033/1‑86 (MQ=255)
                                                                                                                                                  cTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCaa    >  2:3549163/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG  >  NC_000913/1286834‑1287131

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: