Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,282,421 A→T 13.8% intergenic (+439/+101) narI → / ← rttR nitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,282,4210AT13.8% 301.9 / 8.8 149intergenic (+439/+101)narI/rttRnitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript
Reads supporting (aligned to +/- strand):  ref base A (48/70);  new base T (9/11);  total (64/85)
Fisher's exact test for biased strand distribution p-value = 8.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 150 total aligned reads displayed.

TCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGG  >  CP009273/1282276‑1282569
                                                                                                                                                 |                                                                                                                                                    
tCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGACGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCGCACTCGCTGTCGc                                                                                                                                                   >  1:2328185/1‑149 (MQ=255)
tCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGc                                                                                                                                                   >  2:46493/1‑149 (MQ=255)
 ccTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCCACGCTCGCTCGCTTACGCt                                                                                                                                                  >  1:3113694/1‑149 (MQ=255)
 ccTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCt                                                                                                                                                  <  1:92286/149‑1 (MQ=255)
    tCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGa                                                                                                                                               <  2:762748/149‑1 (MQ=255)
    tCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCACTCGCTGTCGCTCgc                                                                                                                                               >  2:4077586/1‑148 (MQ=255)
    tCCCCGCTGGGGCAGAAGATTTGATTGCGGATTCGTTTGAGAAGTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGa                                                                                                                                               <  1:3535909/149‑1 (MQ=255)
      cccGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTGCTGCGGCAACGCTCTCTAACTGACGCACGAGc                                                                                                                                             >  1:1394612/1‑149 (MQ=14)
          cTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGaa                                                                                                                                         <  1:2811671/149‑1 (MQ=255)
           ttGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAAc                                                                                                                                        <  2:4385271/149‑1 (MQ=255)
            tGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCCCCCCAATCACGCAGACGCTCGCGCCGAAcc                                                                                                                                       >  1:2875454/1‑149 (MQ=14)
             gggCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCt                                                                                                                                      <  1:1258643/149‑1 (MQ=255)
              ggCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGGCGGAGGCAACGCACACTCGCAGAAGCTCGAGCCGAAccat                                                                                                                                     >  1:199781/1‑147 (MQ=14)
                 aGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGt                                                                                                                                  <  1:3581245/149‑1 (MQ=255)
                  tattatTTGATTTCGGATTAGTTTGAGTATTCAGGGGATTTTGATATTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTAGGGTCGTTGCCTGCGGCATCGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTc                                                                                                                                 <  1:901690/146‑1 (MQ=14)
                  gAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAAAGCTCTCTCGCTGTCGCTCGAGCCGAACCATTGTc                                                                                                                                 >  2:1922997/1‑149 (MQ=255)
                  gAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCCAACCACACTCGCAGACGCGCGACCCCAACCAAAGTc                                                                                                                                 >  1:1427685/1‑149 (MQ=14)
                      atTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAg                                                                                                                             <  1:46493/149‑1 (MQ=255)
                      atTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACGCTCACTCGCAGACGCACGACCCGAACCAAAGTCCAAg                                                                                                                             >  2:2978662/1‑149 (MQ=14)
                       tttGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTAGCTGTCGCTCGCGCCGAACCGTAGTCGAAGc                                                                                                                            >  2:215987/1‑149 (MQ=255)
                       tttGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACGCTCACGAGCAGACGCACGAGACGAAACATAGTCGaatc                                                                                                                            >  2:175414/1‑147 (MQ=14)
                           aTTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGTAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTtct                                                                                                                        <  2:3471841/149‑1 (MQ=255)
                             tGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCATAGTCGAAGCTTCTCa                                                                                                                      >  2:3699251/1‑149 (MQ=255)
                             tGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGCCGCTCGAGCCCCACCGTAGTCGAAGCTTATCa                                                                                                                      >  2:3980588/1‑149 (MQ=14)
                             tGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCCGCAACCCTCAAACGCTGGCCCGCGAGTCCAACCTTGGTCGTAGcatatac                                                                                                                      >  2:3699254/1‑143 (MQ=11)
                                ggATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCGCGAGCCGAACCTTTGTCGAATCTTCTCAAcc                                                                                                                   >  2:3206572/1‑149 (MQ=255)
                                ggATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACCCTCTCTCGCTGACGCTCGATCCGAACCATAGTCGAAGCTTCTCATcc                                                                                                                   >  2:3420684/1‑149 (MQ=255)
                                 gATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCCACCCTCACACGCTGACGCGCGAGCCGAACCATAATAGAAGCCGCTCaaaca                                                                                                                  >  1:1099620/1‑145 (MQ=14)
                                  aTTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGGCAGCGGCAACACACAAACGCAGACGCACGCGCCGAAACAAAGTCGCAGCTTcacaaaaaa                                                                                                                 >  1:669948/1‑144 (MQ=14)
                                  aTTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCCACGCGCACTCGCGGACGCGCGCCCCGAACCGTAGTCGGAGCGTGTCAtacat                                                                                                                 >  1:2454612/1‑145 (MQ=14)
                                     cGTTTGAGTATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTccc                                                                                                              <  1:2385751/149‑1 (MQ=255)
                                     cGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTccc                                                                                                              <  1:4674323/149‑1 (MQ=255)
                                     cGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTccc                                                                                                              <  1:2052203/149‑1 (MQ=255)
                                     cGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTccc                                                                                                              <  1:2052986/149‑1 (MQ=255)
                                          gagttttCCGGGTCTTTTGTAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCAt                                                                                                         <  1:1409998/144‑1 (MQ=255)
                                           agaATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCATCGGGTCGTTGCCAGCGGCAACGCGCACTCGCGGGCCCTCGAGGCGAACCTAAGGCGAAGCTTCCCATCCTTACCCCCaac                                                                                                        >  2:3654548/1‑147 (MQ=11)
                                            gttttCGGGGTCTTTTGTAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATgg                                                                                                       <  2:3000439/146‑1 (MQ=255)
                                            gaATTCCTGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATgg                                                                                                       <  1:2061008/149‑1 (MQ=255)
                                            gaATTCCGGGGCTTTTGACAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATgg                                                                                                       <  1:3969799/149‑1 (MQ=255)
                                                 ccGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCTCTCTCGCTGACGCACGAGCCGAACCTTAGTCGAAGCTGCTCAACCTTCCCCGCCTGGGccgt                                                                                                  >  2:2544016/1‑146 (MQ=14)
                                                   ggggCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTTGTAGAATCATCTCAAACATCCCCGCATGGGCAGAat                                                                                                >  2:621387/1‑149 (MQ=11)
                                                   ggggCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACCCTCTCTCGCAGCCGCTCGAGACGAACCGTAGTCGAAGCTTCTCAACCTTCCCCCCAAGGGCACAat                                                                                                >  2:3791479/1‑149 (MQ=11)
                                                     ggCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAatat                                                                                              <  1:4196390/149‑1 (MQ=38)
                                                     ggCTTTTGAAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAatat                                                                                              <  1:290241/149‑1 (MQ=38)
                                                      gCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAtt                                                                                             <  1:215987/149‑1 (MQ=17)
                                                       cTTTTGAAAGTGATGGTGGTGGGGGACGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAttt                                                                                            <  1:3968336/149‑1 (MQ=37)
                                                        ttttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCACTCGCTGACGCTCGAGCCGAACCATAGTCGAAGCATCTCAACCTTACCAGCAAAGGCAGAATAGTTg                                                                                           >  2:256247/1‑149 (MQ=11)
                                                        ttttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAAAGCTCACACGCAGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCCCATGGGCAGAAAAttat                                                                                           >  2:3301779/1‑147 (MQ=11)
                                                         atctaaaGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGCCGCTCGAGCCGAACCTTCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAtttaa                                                                                          <  1:271000/145‑3 (MQ=2)
                                                           tgtaaGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGAtt                                                                                        <  1:3121150/146‑1 (MQ=255)
                                                           tGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCTCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGAtt                                                                                        <  2:2428622/149‑1 (MQ=255)
                                                              aaGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACCCTCTCTCGCTGACGCTCGAGCCGAACCATAGTCGAAGCTTCTCATCCTTCCCCGCAAGGGCAGAAAATTTAAtagat                                                                                     >  2:271000/1‑145 (MQ=11)
                                                                gTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCATACCCGCATGGGCAGAATATTTGATTGAGGa                                                                                   >  2:2631307/1‑149 (MQ=255)
                                                                 cgATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGAt                                                                                  <  1:631435/148‑1 (MQ=21)
                                                                 cgATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGAt                                                                                  <  1:631434/148‑1 (MQ=39)
                                                                  gATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCTAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGAtt                                                                                 <  2:2542841/149‑1 (MQ=255)
                                                                  gATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACCCTCACTAGCTGACGCTCGAGCCCAACCATAGTCGAATAATCTAAAACATACCCCCAAGGGCAGAATAATTTATTTAGGAtt                                                                                 >  2:3517441/1‑149 (MQ=11)
                                                                   aTGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGACGAACCTTAGTCGAAGCTTCACATCCATCCCCGCATGGGCCGAAAATGTGATTGCAGAtac                                                                                >  2:4636404/1‑147 (MQ=11)
                                                                     ggtggtggGGGCAGGATTACGCCGCGCTGCGCGCTTCGCCCTTCGGGTCGGTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCCTCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGc                                                                              <  2:3908660/149‑1 (MQ=14)
                                                                      gtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCt                                                                             <  1:664123/149‑1 (MQ=40)
                                                                         gtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTga                                                                          <  1:179558/149‑1 (MQ=255)
                                                                           gggggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCGCGCACGCTGGCGCTCGAGGCGAACCATAGTCGAAGCTTCTCATCCATACCCGCCGGGGCCGAATATTTTATTGCGGTTTCGCtgggga                                                                        >  2:899218/1‑144 (MQ=1)
                                                                             gggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCTAGACGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGt                                                                      <  1:921173/149‑1 (MQ=11)
                                                                             gggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCTCACTCGCAGTCGCTCGAGCCGCACCTTAGTCGAAGCTTCTCATCCATCCCCGCGGGGGCCGAATAATTGATTGCGGGTTCGCGGGAGAGt                                                                      >  2:1149131/1‑149 (MQ=14)
                                                                              ggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGtt                                                                     <  2:255472/149‑1 (MQ=255)
                                                                               gAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTc                                                                    <  1:535977/149‑1 (MQ=255)
                                                                                aaGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCa                                                                   <  1:2289430/149‑1 (MQ=255)
                                                                                 aGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGGCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTTATTGCGGATTCGGTTGGGAGTTGAg                                                                  >  2:312792/1‑149 (MQ=11)
                                                                                  ggATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAgg                                                                 <  1:550102/149‑1 (MQ=255)
                                                                                   gATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCTCACTCGCTGACGCTCGATCCGAACCTGAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCCGAAGATTTGATTGCGGCTTCGCTGGAGAGTTCAggg                                                                >  2:3115881/1‑149 (MQ=14)
                                                                                       aCTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGGGGCAACGCAATCACGCGGACGCGCGAGCCGAACCTTAGTCGCAGCTTCGCATACTTCCCCCCCAGGGCAGAAGATTTTATTGCGGGTTCGCGTGGGAGTTCAGGGCCtt                                                            >  2:1336862/1‑149 (MQ=14)
                                                                                              gcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAg                                                     <  1:4241630/149‑1 (MQ=255)
                                                                                               cTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGt                                                    <  1:860731/149‑1 (MQ=255)
                                                                                               cTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGt                                                    <  2:3280746/149‑1 (MQ=255)
                                                                                                  cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAt                                                 <  1:2170881/149‑1 (MQ=255)
                                                                                                       ttCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCCAGTCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggtgg                                            <  2:2789007/149‑1 (MQ=11)
                                                                                                        tCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggtggt                                           <  1:2761894/149‑1 (MQ=255)
                                                                                                           cccTTCGGGTCGTTGCCTGGGGCAACGCTCACACGCGTTCGCTCGTGCCGAACCTTAGTGGAAGCTTATCAACCTTACCCGCAGGGGCCGAAAATTTGATTGCGGGATCGCTTTTGAGTTGAGGGACTTTTGAATGTGTTGgggggggg                                        >  2:4400174/1‑149 (MQ=9)
                                                                                                            ccTTCGGGTCGTTGCCGGCGGCAACCCTCACACGCTGACGCTCGAGCCGAACCATAGTCGAAGCATCTCAAACTTCCCCGCAAGGGCACAATATTTTATTGCGGATTCGCTTGTGAGTTCAGGGACTTTTGAAAGTGATGGAGGTgggg                                       >  2:1390115/1‑149 (MQ=14)
                                                                                                              ttCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGa                                     <  1:1456389/149‑1 (MQ=255)
                                                                                                                cGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAg                                   <  2:2811544/149‑1 (MQ=255)
                                                                                                                 gggTCGTTGCCTTCGGCAACGCTCTCTCGCTGACGCGCGAGCCGAACCTTAGTCGAATCTTCTCATACATACCCGCAGGGGCAGTATATTTTATTGCGGATTGGGTTGGGAGTGAAGGGACGGGTGTAAGTGTTGGTGGTGGGGGGAgg                                  >  2:3074960/1‑149 (MQ=255)
                                                                                                                 gggTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGACCCTGAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGAGTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAgg                                  <  2:1440593/149‑1 (MQ=255)
                                                                                                                  ggTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCCAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGa                                 <  2:1644869/149‑1 (MQ=255)
                                                                                                                   gTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAt                                <  2:500763/149‑1 (MQ=255)
                                                                                                                      gTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGACTATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTAc                             <  1:1668296/149‑1 (MQ=255)
                                                                                                                      gTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTACCCGCAGGGGCCGAATATTTTATTGCGGTTTGGGTTGTGAGATCAGGGACTTTTGCAAGTGTTGGTGGGGGGGGAAGGAttgg                             >  1:546259/1‑147 (MQ=14)
                                                                                                                       ttGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGACTATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGACGGATTagt                            >  2:2467679/1‑147 (MQ=255)
                                                                                                                           cTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGATTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCCtacgc                        >  2:105991/1‑145 (MQ=255)
                                                                                                                           cTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTGAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAAGATGTGAGTGCGGCTTCGCTTGAGAGGTCAGGGACGTGTGTATGGGAGGGTGGTGGGGGACGGAGTAATCAgg                        >  2:1356323/1‑148 (MQ=17)
                                                                                                                               ggCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTGTTGAAAGTGATCGTCGTGGGGGAAGGAATACTCAGCGCTg                    >  1:566737/1‑149 (MQ=255)
                                                                                                                                 cAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAACTGATGGGGGTGGGGGAAGGATTACTcagcgcagcg                  >  1:2910066/1‑149 (MQ=255)
                                                                                                                                  aaCGCTCTCTCGCTGACGCTCGTGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGATTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTgcgc                 <  2:2144559/149‑1 (MQ=255)
                                                                                                                                       tctctcGCTGACGCTCGAGGCGATCCTTCTTCGAAGCTTCTGAGCGTTCCCCGCTTGGGCAGATTATTTGATTGCGTATTAGCTTGAGAGTTCCGGTACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTc            <  2:1768429/149‑1 (MQ=11)
                                                                                                                                       tctctcGCTGACGCTCGAGGCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCTTTGGCAGAATATTTGATTGCGGATTCGCTTGAGATTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTc            <  2:1768519/149‑1 (MQ=255)
                                                                                                                                            cGCTGTCGCTCTAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCt       <  2:185014/149‑1 (MQ=255)
                                                                                                                                              cTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCCGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTc     >  1:4559291/1‑149 (MQ=255)
                                                                                                                                              cTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTc     >  1:4559749/1‑149 (MQ=255)
                                                                                                                                               tGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCg    >  2:3979195/1‑149 (MQ=255)
                                                                                                                                                 aCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCggg  >  2:4047657/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGG  >  CP009273/1282276‑1282569

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: