Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,287,007 G→A 6.6% intergenic (+481/+59) 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,287,0070GA 6.6% 211.7 / 4.9 77intergenic (+481/+59)narI/rttRnitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript
Reads supporting (aligned to +/- strand):  ref base G (42/28);  new base A (2/3);  total (45/32)
Fisher's exact test for biased strand distribution p-value = 6.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCG  >  NC_000913/1286859‑1287154
                                                                                                                                                    |                                                                                                                                                   
tttGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCGAGAGACGAACCGTGGTCGAAGCTTCAAATCCATCCCCGc                                                                                                                                                    >  2:2543569/1‑150 (MQ=11)
     gaATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCCACCCACTCACGCAGAAGCACCAGCCGAACCAAAGTCGCAGCTTATCATACTTCCCCGCATggg                                                                                                                                               >  1:715852/1‑150 (MQ=14)
            ggggCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGaaaa                                                                                                                                        >  2:584514/1‑148 (MQ=39)
                   ttgtctgtgtTGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCACATGGGCAGAATATTTGAtt                                                                                                                                 <  1:2463726/144‑1 (MQ=255)
                   tggtctgtgtTGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGAtt                                                                                                                                 <  1:1803673/144‑1 (MQ=255)
                         gTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGAt                                                                                                                           <  1:2060798/150‑1 (MQ=37)
                          tGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGAtt                                                                                                                          <  2:2464217/150‑1 (MQ=39)
                             tggtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGTGTCGTTGCCTGCGGCAACGCTATCACGCTGACGCTCGAGCCGTACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGAttagg                                                                                                                       >  2:1669931/1‑147 (MQ=14)
                                  gtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCACATGGGCAGAATATTTGAtt                                                                                                                                 >  2:424873/1‑135 (MQ=21)
                                  gtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCACATGGGCAGAATATTTGAtt                                                                                                                                 >  2:424947/1‑135 (MQ=255)
                                  gtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCACATGGGCAGAATATTTGAtt                                                                                                                                 <  1:424947/135‑1 (MQ=255)
                                  gtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCACATGGGCAGAATATTTGAtt                                                                                                                                 <  1:424873/135‑1 (MQ=255)
                                   tggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTgaga                                                                                                                 <  2:522458/150‑1 (MQ=255)
                                        gAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGGTGGGGAGTTCa                                                                                                            >  2:359641/1‑150 (MQ=255)
                                        gAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCa                                                                                                            >  2:2459511/1‑150 (MQ=255)
                                         tatgATTACTCCGCGCCGCGCGATTCGCCCTTCGGGTGGTTGCCTGCGGCAACGCTCTCCGGCCGTCGCACGAGCCGAACCTTCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAg                                                                                                           <  2:893840/147‑1 (MQ=14)
                                           ggATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAATCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAggg                                                                                                         <  2:1734822/150‑1 (MQ=255)
                                            gATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCCGACGATTTGATTGCGGATTCGCTTGAGAGTTCAGGGa                                                                                                        >  2:2258431/1‑150 (MQ=255)
                                                 cTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACtttt                                                                                                   <  1:2069173/150‑1 (MQ=255)
                                                 cTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACtttt                                                                                                   <  1:1082087/150‑1 (MQ=255)
                                                  tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAGCCGTAGTCGAAGCGTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTg                                                                                                  >  2:213485/1‑150 (MQ=255)
                                                  tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGGAGCTTCTAATCCTTCCCCGCATGGGCACAATATTTGATTGCGGATGCGCTTGGGAGTTCAGGGACTTTTg                                                                                                  >  2:790169/1‑150 (MQ=14)
                                                    aGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCATTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGaa                                                                                                <  1:741139/150‑1 (MQ=255)
                                                     gcgctgcgctCTTCGCACTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCAGAACCTTAGTCGAAGCTTCTCATCCTTCCCCTCATGGGCAGAATATTTGATTGCGGATTCGCTTTAGAGTTCAGGGACTTTTGaaa                                                                                               <  1:1752236/150‑1 (MQ=17)
                                                         tGCGCGCTTAGCCCTAGGTGTCGTTGCCTGCGGCATCGCTCTCTCGCTGACGCTCGTGCGGATCGTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGa                                                                                           <  1:1704248/150‑1 (MQ=14)
                                                           cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtg                                                                                         <  1:2260998/150‑1 (MQ=255)
                                                           cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtg                                                                                         <  1:2260933/150‑1 (MQ=255)
                                                           cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtg                                                                                         <  2:1537352/150‑1 (MQ=255)
                                                                ttCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCg                                                                                                                        <  1:1596268/114‑1 (MQ=255)
                                                                ttCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCg                                                                                                                        >  2:1596268/1‑114 (MQ=255)
                                                                ttCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggtggt                                                                                    >  2:2094050/1‑150 (MQ=255)
                                                                ttCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTAGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAATTTATTGCGGATTCGGTTGGGAGTTCAGGGACTTTTGAAAGTGAtggttgg                                                                                    >  2:804024/1‑147 (MQ=14)
                                                                 tCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCTTCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggtggtg                                                                                   <  2:147742/150‑1 (MQ=255)
                                                                 tCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggtggtg                                                                                   <  1:1035033/150‑1 (MQ=255)
                                                                       ttCGGGTCGTTGCGTGCGGCCCCGCTCTCTCGCTGCCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGCTTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGaa                                                                             <  1:774705/150‑1 (MQ=14)
                                                                           gagcGTTGCCGCAGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCGTACCCGCAAGGGCCGAAGATTTGTTTGCGGATTCGGTTGTGAGTTGAGGGACTTTTGTAAGTGATGGTGGTGGGGGTAggct                                                                         >  1:662841/4‑148 (MQ=255)
                                                                           gagcGTTGCCGCAGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCGTAGTCGAAGCTTCTCAAACATCCCCGCAGGGGGAGAGTATTTGATTGCGGGTTAGCTTGAGAGTTGAGGGACTTTTGTGAGTGTTGGTGGGGGGGGGAgggt                                                                         >  1:662839/4‑148 (MQ=255)
                                                                            gTCGTTGCCTGCGGCAACTCTCTCTCGCTGTCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGATTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAtt                                                                        <  2:2055361/150‑1 (MQ=255)
                                                                            gTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAtt                                                                        <  2:2055767/150‑1 (MQ=255)
                                                                              cGTTGCCGCAGGCAACGCTCTCTAGCTGACGCTCGAGCCGAACCATAGTAGAAGCATATCATACATCCCCGCATGGGCAGAAAATTTGATTGCGGATTCGCTTGAGAGTTAAGGGAATTTTGAAAGTGATGGTGGTGGGGGAAGGATTaa                                                                      >  1:1788919/1‑149 (MQ=17)
                                                                               gTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACt                                                                     <  2:1266094/150‑1 (MQ=255)
                                                                                ttGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGGACCTTAGTCGAAGCTTCTCATCCTTACCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACGTTTGTAAGTGATGGTGGTGGGGGAAGGATTACTc                                                                    >  1:58775/1‑150 (MQ=255)
                                                                                ttGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGTGAGTTCAGGGACTTTTGTAAGTGATGGTGGTGGGGGCAGGATTACTc                                                                    >  1:1361240/1‑150 (MQ=255)
                                                                                ttGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCGTAGTCGAAGCTTCTCATCCGTCCCCCCAAGGGCCGAATATTTGATTGCGGCTTCGCGTGGGAGGTCAGGGCCGGTTGAAAGTGGTGGTGGGGGGGGACGGATTACTc                                                                    >  1:730144/1‑150 (MQ=14)
                                                                                ttGCCTGCGGCAACGCTCTCTCGCTGACGCTCCAGCCGAACCTTAGTCGAACCTTCTCATCCTTCCCCGCAAGGGCAGAAAATAACAATGCGGATTCCCTTCAAAGTTCAGGGACATTAGAAAGTGATGGTGGTGGGGGAAGGATTACTc                                                                    >  2:2499453/1‑150 (MQ=14)
                                                                                   ccTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGTTTGGGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAttcgtccga                                                                 >  2:1221287/1‑143 (MQ=255)
                                                                                    cTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgcg                                                                >  2:1745524/1‑150 (MQ=255)
                                                                                         gCAACGCTCTCTCGCTGCCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTgcg                                                           <  2:2192391/150‑1 (MQ=255)
                                                                                                    tcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCt                                                >  2:1604634/1‑150 (MQ=255)
                                                                                                        tgtctcTCCAGCCGAACCTCTGTCGAAGCGTATCCTCCTTCACCGCTTGGGCATATTTTTTGATTGCGGATTCGCTTGAGATTTCTGTGACTTTTGTAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCgg                                            <  2:2066416/145‑1 (MQ=11)
                                                                                                        tggcGCTCGAGCCGAACCTTAGTCGAAGCGTCTCACCCTTCCCCGCATGGGCAGAATCTTTGATTGCGGATTCGCTTGAGCGTTGTGTGACTTTTGAAAGGGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCgg                                            <  2:781525/147‑1 (MQ=14)
                                                                                                            gCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACACAGAGCTGAGCCCTTCTCGCTTCGGGtag                                        >  2:179036/1‑148 (MQ=17)
                                                                                                             cTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGt                                       >  2:2092012/1‑150 (MQ=255)
                                                                                                             cTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGt                                       >  2:2299127/1‑150 (MQ=255)
                                                                                                              tCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGATAGTGATGGTGGTGGGGGTAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGGCGtt                                      >  2:566455/1‑150 (MQ=255)
                                                                                                                gAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGTGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGATCTGAGACCTTCGCCCGTCGGGTCGATGc                                    >  2:1718472/1‑150 (MQ=14)
                                                                                                                gAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGATGc                                    >  2:1449179/1‑150 (MQ=255)
                                                                                                                  gCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGccg                                  >  1:880745/1‑149 (MQ=255)
                                                                                                                  gCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGGGGTGGGGGACGGATTACTCAGCGCTGCGCGCTTCGCCCAGCGGGTCGTGGCCt                                  >  2:2448723/1‑150 (MQ=14)
                                                                                                                    cGACCCTTGGTCGAAGCTTCGCCTCGTTCCACGCATGGGCAGCTTTTTTGATTTGGGATTCGCTTGAGCTTTCAGGTACTTTTGCAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGc                                <  2:1626029/150‑1 (MQ=11)
                                                                                                                    cGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGATGCCTGc                                >  1:2386406/1‑150 (MQ=255)
                                                                                                                    cGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGAGCGCTTCGCCCTTCGGGTCGtgtccagc                                >  1:1105349/1‑143 (MQ=255)
                                                                                                                       aCCTTAGTCGAATCTTCTAATCCTTCCCCGCATGGGCAGACTATCTGCTTGCTGTTTCGCTCGAGAGTTCAGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGc                             <  2:715852/150‑1 (MQ=14)
                                                                                                                         cTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCaa                           >  1:2360641/1‑150 (MQ=255)
                                                                                                                                    cctctcATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGc                <  2:688160/148‑1 (MQ=255)
                                                                                                                                      tctcATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTAGTCAGCGCTGCGCGCATCGCCCTTCGTGTCGTTTCGGGCGGCAAAGGTCTCTCggag              >  2:2345735/1‑147 (MQ=14)
                                                                                                                                      tctcATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTg              >  2:261775/1‑150 (MQ=255)
                                                                                                                                      tctcATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTg              >  2:2372531/1‑150 (MQ=255)
                                                                                                                                      tctcATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCGCTCTCGCTg              >  1:2032610/1‑150 (MQ=255)
                                                                                                                                       ctcATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAttcg                                                                      <  1:1924226/93‑3 (MQ=255)
                                                                                                                                       ctcATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAttcg                                                                      >  2:1924226/1‑91 (MQ=255)
                                                                                                                                        tcATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGgc            >  2:463270/1‑150 (MQ=255)
                                                                                                                                         cATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATTGTGGTGGGGGAAGGATTACTCTGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGgcg           <  1:213485/150‑1 (MQ=255)
                                                                                                                                              ttCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCACACTCGCAGGCGCACGa      >  2:2414424/1‑150 (MQ=255)
                                                                                                                                              ttCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCGTCCGGTCGTGGCAAGCGGCAACGCACACTCGCACGCGCTCGa      >  2:1416134/1‑150 (MQ=14)
                                                                                                                                                  ccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCg  >  1:1642871/1‑150 (MQ=255)
                                                                                                                                                  ccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCGGGCGCTCGAGTCg  >  1:1668474/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCG  >  NC_000913/1286859‑1287154

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: