Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,286,612 A→T 16.4% intergenic (+86/+454) 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,6120AT16.4% 149.6 / 2.5 60intergenic (+86/+454)narI/rttRnitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript
Reads supporting (aligned to +/- strand):  ref base A (16/33);  new base T (9/2);  total (25/35)
Fisher's exact test for biased strand distribution p-value = 5.11e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGATGGTGTGGCGTTTATCTTCCGTCTGCACCTGGTGCTGGGGATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAA  >  NC_000913/1286379‑1286840
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                       
cGATGGTGTGGCGTTTATCTTCCGTCTGCACCTGGTGCTGGGGATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATTGGCACGGGTGGGGGGGGGTTGTCGGGTGCGGTTCCCTTTTCCgg                                                                                                                                                                                                                                   >  2:644300/1‑238 (MQ=255)
    ggttttgtgTTTTTCTTCCGTCTGCCCCTGGTGCTGGGGATGACGTTATTCCTGCGGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTa                                                                                                                                                                                                                               <  2:554757/231‑1 (MQ=255)
        tgGCGTTTATCTTCCGTCTGCACCTGGTGCTGGGGATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTTTGGGCACGGGGCGGGGGGAGTGTTCGGTTGCGCTTCGCTTATCCGACCTCCAgg                                                                                                                                                                                                                           >  1:19203/1‑239 (MQ=255)
         ggCGTTTATCTTCCGTCTGCACCTGGTGCTGGGGATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAggg                                                                                                                                                                                                                          <  1:644300/239‑1 (MQ=255)
          gCGTTTATCTTCCGTCTGCCCCTGGTGCTGGGGATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAggg                                                                                                                                                                                                                          <  1:729527/238‑1 (MQ=255)
               tATCTTCCGTCTGCACCTGGTGCTGGGGATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGTGCTGGGTGAGTTGTCGGGTGCGCCTCCCTTTTCCGGCCTTCAGGgggggag                                                                                                                                                                                                                    >  1:252787/1‑238 (MQ=255)
               tATCTTCCGTCTGCACCTGGTGCTGGGGATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGACGGGGTTTTTTTATGGGCACGGTGCGGGGGGAGTTTTCGGAAGCGCTTCGCTTATACGACCCAAAGGggcggga                                                                                                                                                                                                                    >  1:170335/1‑237 (MQ=255)
                aTCTTCCGTCTGCACCTGGTGCTGGGGATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGTTGCGCTTCCCTTTTCCGACCTACAGGGGAGGatt                                                                                                                                                                                                                   >  1:554757/1‑238 (MQ=255)
                        tCTGCACCTGGTGCTGGGGATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTTGAGGATCTGACACGTACGTACCAGCTGGTGCGCGCCCGTCACTAAGCGCATTTTAGTTCACATAGACCCGGCTTCGGCGGGGGGTTTTTCATGGTCCGGTCCGGGGGGAGGGGTGGGATGCGGTTAGCTTTACCGACCGACGAGGGGGGATAGTgaaagg                                                                                                                                                                                                           >  2:605451/1‑234 (MQ=255)
                          tGCACCTGGTGCTGGGGATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGGTGGGCTTCGCTTTTCCGACCTCAAGGGGGCGGTTTTTTAGGccc                                                                                                                                                                                                         >  1:846128/1‑239 (MQ=255)
                             cccgtgTGCTGGGGATTGCGTTACTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCgg                                                                                                                                                                                                       <  1:488988/233‑1 (MQ=255)
                                   tgttGGGGGTGGCGTTTTTCCTGCGTTTCCCGTTTTCGCGTCTGATACACATCTGGGGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGGGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACg                                                                                                                                                                                                 <  2:756076/235‑1 (MQ=255)
                                             gACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACAGTGCGGGGTGAGTTGTCGGATGAGCTTCGCTTTTCCGACCTACAGGGGGGGTGATTGTGGGGCCGGGAAACGGGGTGcccaag                                                                                                                                                                                     >  2:441219/1‑235 (MQ=255)
                                                gTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGgga                                                                                                                                                                                                                        <  1:777293/202‑1 (MQ=255)
                                                          gtgTTCCCGTTCTCGGGTTTTATACACATCGGGAGGGTACCGGTGGAGTATCGGACACGTAAGTACCAGCGGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCt                                                                                                                                                                       <  2:511813/238‑1 (MQ=255)
                                                          ctgTTCCCGTTTCCGCGTTTGGTTCCCATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCt                                                                                                                                                                       <  2:222010/239‑1 (MQ=255)
                                                          ctgTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATTGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCCCTTTTCCGACCTACAGGGGGGGATATTGTAGGCCCGGTAAACCGGGGGCCATCGGGCAGAATTAGc                                                                                                                                                                        >  1:23639/1‑238 (MQ=255)
                                                           tgTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTCGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCt                                                                                                                                                                       <  1:826164/238‑1 (MQ=255)
                                                            gTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCg                                                                                                                                                                     <  1:778441/239‑1 (MQ=255)
                                                                       cgcgTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACCGGGGAGGATATTGTAGGCCCGGGAAACGTGGGGCCTTCGGGCCAAATTTGCTCGGGGGAACAggg                                                                                                                                                          >  1:696775/1‑238 (MQ=255)
                                                                            agagacagcaTCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGa                                                                                                                                                     <  1:52031/231‑1 (MQ=255)
                                                                                     ctttGGGGCGTTCCGGTGGAGTATCTGAACCGTAAGTACCAGCTGGTGGGCGCTCGTCACTAATGGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAg                                                                                                                                            <  2:751763/236‑1 (MQ=255)
                                                                                     aTCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAAc                                                                                                                                                                                                  <  1:207874/187‑1 (MQ=255)
                                                                                     aTCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCt                                                                                                                                                                       <  1:195136/212‑1 (MQ=255)
                                                                                               tACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGGGGGTTGTCGGATTCGCTTCGCTTATCCGCCCCCACGGGGGGGGTATTGTAGGCCCGGGGAAAGTGGGGCCCATCGGACAAAATTGCCCGAGGTAAACGGTTTTAATTTTTTAAGgaagggaagg                                                                                                                                  >  2:232938/1‑239 (MQ=255)
                                                                                               tACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTCTTTTTATCGGCACGGTGCGGGGTGCGTTGTCGGGGGTGCTTCGCTTTTCCGACCTACCGGGGGGGGTTTTGTAGGCGCGCGCAACGCGGTGCCATCGGGCACAATTAGCCCGGGGGAACAGGTTTTTTTTATTTGAAGAAGGGagag                                                                                                                                  >  2:777848/1‑236 (MQ=255)
                                                                                               tACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGGTTTTTTATGGGCACGGTGCGGGGGGAGTTGTCGGGTGCGCTTCGCTTATCCGGCCTACAGGGGAGGATATTGTAGGCCCCGTAAACGTGGGGCCAACGGGCAAAATTAGATCGGGGTAACAGGGTTTGATTATTTAAGaaatggtcct                                                                                                                                  >  2:52031/1‑232 (MQ=255)
                                                                                               tACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAAATTTAGTTCACCTAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGCGAGTTTTCGGCTGGGCCTCGCCTATCCGGACCAACAGGGAGGATATTTTAGGCCCGGTCAACGTGGTGGCACTCGGCGAAATTTACCGCGGGTAAACGGTTTTGTTTTTTTAAGGAAGCggggg                                                                                                                                  >  2:731739/1‑235 (MQ=255)
                                                                                               tACCGGGGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGAtg                                                                                                                                   <  1:884447/238‑1 (MQ=255)
                                                                                                aCCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCt                                                                                                                                                                       <  1:191758/201‑1 (MQ=255)
                                                                                                  cGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGAtggtgg                                                                                                                               <  1:232938/239‑1 (MQ=255)
                                                                                                   gggtgggTATTTGACACGTAAGTTCCAGCTGGTGCGCGCTAGTCACTAAGGGAATTTTAGTTGGCATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGAtggtggt                                                                                                                              <  2:23639/233‑1 (MQ=255)
                                                                                                         gtttttttcccgttatttccAGCTGGTGCGCGCTTGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTggggg                                                                                                                         <  2:393131/220‑1 (MQ=255)
                                                                                                         gtttctgcaccGTTAGTACCCGCTGGTGCGCGCTCGGTATTAAGGGAATTTTAGTTCACATAGGCCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTggggg                                                                                                                         <  2:252787/228‑1 (MQ=255)
                                                                                                         gTATCTGACAGGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTggggg                                                                                                                         <  1:631088/238‑1 (MQ=255)
                                                                                                             cTGACACGTAAGTCCCAGCTGGGGCGCTGGCGTCCCTAAGTGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGa                                                                                                                    <  1:6885/239‑1 (MQ=255)
                                                                                                                    gTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGGTTTGATTATTTAAGGAAGCGCTGGGGGTGGGGGAAGGGTTACTCa                                                                                                             >  1:576086/1‑239 (MQ=255)
                                                                                                                                   tGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCCCGGTGCGGGGGGAGTTGTCGGGTGCGCTTCGCTTTTCCGCCCTCCAGGGGGGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCACAATTAGACCGGGGTAACCGGTTTTGATTATTTAAAGAAGGGGTGGTGGTGGGGGGAGGAGTTATCAGCGCcgcgcggtcg                                                                                               >  2:209370/1‑234 (MQ=255)
                                                                                                                                         gtCGTCCCTAAGCGAATTTTAGTTCCCATAGGCCCTGATTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTc                                                                                        <  1:605451/238‑1 (MQ=255)
                                                                                                                                                      gAATTTTAGTTCACATGGACCCTTCTTCGGCGGGTTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGcc                                                                             <  1:506148/237‑1 (MQ=255)
                                                                                                                                                               gtCACATAGGCCCTGCTTCGGCGGGGTTTTTTTAGGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGc                                                                  <  2:696775/238‑1 (MQ=255)
                                                                                                                                                                            tGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGGTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACCCTCTCTCGCTGgc                                                      >  1:512163/1‑238 (MQ=255)
                                                                                                                                                                             gCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGa                                                                                                                                                     <  2:521045/142‑1 (MQ=255)
                                                                                                                                                                             gCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGa                                                                                                                                                     >  1:521045/1‑142 (MQ=255)
                                                                                                                                                                             gCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGGGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTCGCTGCGGCaaacaccctccgcctgcg                                                     >  1:761560/1‑222 (MQ=39)
                                                                                                                                                                                  ggcggGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAg                                               <  1:209370/239‑1 (MQ=255)
                                                                                                                                                                                         gTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGGTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTTCGGCCACGCCCTCTCCCTGGCCCCCAAGTCGaa                                          >  1:6513/1‑238 (MQ=255)
                                                                                                                                                                                          tttttttATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAAc                                         <  1:557829/238‑1 (MQ=255)
                                                                                                                                                                                                taggggCACGGTGGGGGGGGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTc                                  <  1:441219/236‑1 (MQ=255)
                                                                                                                                                                                                                    aGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGa                                                                                                                                     <  1:858761/120‑1 (MQ=255)
                                                                                                                                                                                                                    aGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGa                                                                                                                                     >  2:858761/1‑120 (MQ=255)
                                                                                                                                                                                                                       tggcGGATGCGCTTGTCTTTTGCGACCTACAGGGGAGGATTTTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCGCGGAGTAACTGGTTTGGATTATTGAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGc            <  2:691956/235‑1 (MQ=21)
                                                                                                                                                                                                                       tGTCGGGTTGGCTTCGCTTTTCCGACCGACAGGGGAGGGTTTTGTTGGCCCGGTAAACGGGGTGGCAGCGGGAAAAATTGGCTCGGAGTAACGGGTTTTGATTTTTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCCGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGc            <  2:585404/238‑1 (MQ=18)
                                                                                                                                                                                                                           ggtctcgcttcgtttATCCGGCGTAGAGGGGAGGATTTGGTAGGCCTGGTTAACGTTGTGCCATCGGGCAATATTAGCTCGGAGTAACAGGTTTGTATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTggg       <  2:250809/234‑1 (MQ=21)
                                                                                                                                                                                                                              tGCGCTTGGCTTATCCGACCTAAGGGGGGTGATTTTGGAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCa     <  1:908161/238‑1 (MQ=255)
                                                                                                                                                                                                                              tGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCGGGTAAAGGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCa     <  1:410381/238‑1 (MQ=255)
                                                                                                                                                                                                                              tGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGGACCTTGGTCGAAACTTCTCATCCTTCCCCGCTTGGGCAg    >  1:585404/1‑239 (MQ=255)
                                                                                                                                                                                                                              tGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAg    >  1:133795/1‑239 (MQ=255)
                                                                                                                                                                                                                                cgcttcgcttATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTTCTCGCTTCGCCCTTCGGGTCGTAGCCTGCGGCACCGCTCTCTCGCCGGCGCTCAAGTCGACCCTTGGTCGAAGCTTCCCCTCCTTcccccgtgggggggaa  >  2:908161/1‑227 (MQ=25)
                                                                                                                                                                                                                                cgcttcgcttATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGTCCGTTGCCTGCGGCAACGCCCTCTCGCTGGCGCTCAAGTCGAACCTTGGCCGAACCTTTCCTTCCTTCCCCCCTTTGGCGGag  >  1:202456/1‑238 (MQ=25)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                       
CGATGGTGTGGCGTTTATCTTCCGTCTGCACCTGGTGCTGGGGATGACGTTATTCCTGCTGTTCCCGTTCTCGCGTCTGATACACATCTGGAGCGTACCGGTGGAGTATCTGACACGTAAGTACCAGCTGGTGCGCGCTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGTCGGATGCGCTTCGCTTATCCGACCTACAGGGGAGGATATTGTAGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAA  >  NC_000913/1286379‑1286840

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