Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,282,421 A→G 8.3% 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,4210AG 8.3% 218.1 / 4.7 96intergenic (+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 (33/51);  new base G (5/2);  total (41/55)
Fisher's exact test for biased strand distribution p-value = 1.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGG  >  CP009273/1282282‑1282568
                                                                                                                                           |                                                                                                                                                   
cccGCTTGGGCAGAAGATTTGAGTGCGGATTCGTTTGAGAATTCGGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAg                                                                                                                                             <  2:2778266/148‑1 (MQ=255)
    cTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGa                                                                                                                                         >  2:1162086/1‑148 (MQ=255)
    cTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGa                                                                                                                                         >  2:3412995/1‑148 (MQ=255)
      tGGGGAGATTATTTGCTTGGGGCTTAGTTTGAGTATTGGGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAAc                                                                                                                                       <  2:2537798/148‑1 (MQ=14)
         gCAGAATATATGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCtt                                                                                                                                    <  2:2594597/148‑1 (MQ=255)
          cAGATTATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTa                                                                                                                                   <  1:1271140/148‑1 (MQ=255)
          cAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCttg                                                                                                                                   >  2:2336946/1‑147 (MQ=255)
          cAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTa                                                                                                                                   <  1:1241642/148‑1 (MQ=255)
          cAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTa                                                                                                                                   <  1:1247000/148‑1 (MQ=255)
          cAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTAGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTa                                                                                                                                   >  2:2336990/1‑148 (MQ=255)
           aGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAg                                                                                                                                  <  1:356112/148‑1 (MQ=255)
             aaTATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCCCAGTCGCTCGAGGCGAACCGTAGta                                                                                                                                >  2:2963333/1‑147 (MQ=14)
               tagtttaTTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGa                                                                                                                              <  2:2321693/142‑1 (MQ=255)
                atTTGTTTGCGGCTTCGTGTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGaa                                                                                                                             <  1:1337772/148‑1 (MQ=255)
                atTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGGCTGCGGCCACGCTCTCTCGCTGGCGCACCAGACGAACCTTAGTCGaa                                                                                                                             >  2:963281/1‑148 (MQ=14)
                atTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCCACGCACTCTCGCTGACGCGCGAGCCGAACCCTAGTCGag                                                                                                                             >  2:1263032/1‑147 (MQ=14)
                 tttGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTGCGTCGGCAAAGCTAAAACGCTGACGCAACAGCCCAACCCAAGTCGAAg                                                                                                                            >  1:3210205/1‑148 (MQ=14)
                 tttGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGGCTGCGGCAACCCTCTCTCGCTGAAGCTCGAGCCCAACCGTAGTCGaat                                                                                                                            >  1:3208083/1‑147 (MQ=14)
                 tttGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGGCTGCGGCAACCCACTCACGCGGACGCGCGAGCCGCACCGTTGTCGGAg                                                                                                                            >  1:2857973/1‑148 (MQ=14)
                 tttGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCCGCAACCCACTCACGCTGACGCTCGAGCCGAAACTTAGTCGAAg                                                                                                                            >  2:1692512/1‑148 (MQ=14)
                 tttGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCGTTGTCGAAg                                                                                                                            >  1:2857995/1‑148 (MQ=255)
                  ttGATTGCGGATTCGTATGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGc                                                                                                                           <  2:989296/148‑1 (MQ=255)
                    gATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCtt                                                                                                                         >  2:271298/1‑148 (MQ=255)
                     aTTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTtc                                                                                                                        <  1:2620109/148‑1 (MQ=255)
                     aTTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTtc                                                                                                                        <  1:1726660/148‑1 (MQ=255)
                     aTTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAAACCACTCTCGCTGACGCTCCAGCCGAACCATAGTAGAAGcata                                                                                                                        >  2:1192565/1‑145 (MQ=14)
                       tGCGGATTCGTTTGAGAATTCCTGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTtctc                                                                                                                      <  2:146609/148‑1 (MQ=255)
                       tGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACCCACTCTCGCTGTCGCGCGAGACGAACCATAGTCGAAGCtaaga                                                                                                                      >  2:1014688/1‑144 (MQ=14)
                         cGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCAt                                                                                                                    >  2:2138750/1‑148 (MQ=255)
                            aTTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAg                                                                                                                                             <  1:2880304/120‑1 (MQ=255)
                            aTTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAg                                                                                                                                             >  2:2880304/1‑120 (MQ=255)
                            aTTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCCAAGCGCTCGCGCGGTCGCGCGCGCCCTCCCGTAGTCGCACCATCTCATCCt                                                                                                                 >  1:2205281/1‑148 (MQ=14)
                              tCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCCACGCGCTCTCGGAGACGCGCGCGCCGAACCTGAGTCGATGCTTATCATCCTTc                                                                                                               >  2:945106/1‑148 (MQ=14)
                                  ttGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCGGAGGCAACCCACACTCGCTGACGCTCGAGCCGAACAATAGTCGAAGCTTATCATACAACCCCg                                                                                                           >  1:2019397/1‑148 (MQ=14)
                                  ttGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACCCTCTCTAGGTGAAGCGCCAGCCGAACCATAGTAGAACCTTGTCATCCATacccc                                                                                                           >  1:790748/1‑147 (MQ=14)
                                  ttGAGAAGTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCg                                                                                                           <  1:2245094/148‑1 (MQ=255)
                                     agaATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCAt                                                                                                        <  1:2989826/148‑1 (MQ=255)
                                      gtaTTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATg                                                                                                       <  1:2104257/146‑1 (MQ=255)
                                       aaTTGCTGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATgg                                                                                                      <  1:123117/148‑1 (MQ=255)
                                         ttCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGGCGGCGGCAACGCGCTCACGCGGAATCGGGAGCCGCACCTGTGTCGTAACTTGTGAAACTTCCCCGCATgggg                                                                                                    >  2:1503493/1‑147 (MQ=14)
                                              tgtgctttGACTGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAat                                                                                               <  1:2524925/143‑1 (MQ=255)
                                              gggCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAat                                                                                               <  2:1892630/148‑1 (MQ=39)
                                                 cTTTTTAAAGTGTTGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAtt                                                                                            <  1:3167719/148‑1 (MQ=255)
                                                 cTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAtt                                                                                            >  2:2643446/1‑148 (MQ=37)
                                                 cTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCACTCTCGCTGACGCACGAGCCGAACCTTAGTCGAAGCTTCTCATCCATCCCCCCAAGGGCAGAAAAtt                                                                                            >  2:2643349/1‑148 (MQ=11)
                                                 cTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAAAGCACTCTCGCTGACGCTCGAGCCCAACCGTAGTCCAAGCTTCTCATACATACCCGCAAGGGCAGAATAtt                                                                                            >  2:2643292/1‑148 (MQ=11)
                                                   tttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCTAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTg                                                                                          <  2:3412745/148‑1 (MQ=38)
                                                     tGATAGTGATGGTGGTGGGGGATGGATTACTCAGCGCTGCGCGCTTCGCGCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGAt                                                                                        <  1:2389389/148‑1 (MQ=255)
                                                       aaaGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACCCTCTCTCGCTGACGCTCGAGCCGAAACTTAGTAGAAGCTTCTCATCCATACCCCCATGGGCAGGATAtttcaatg                                                                                      >  2:1976030/1‑143 (MQ=11)
                                                          gTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCgg                                                                                   <  2:355998/148‑1 (MQ=39)
                                                             aTGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGGCGGCGGCCAACCTCTCAAGCTGACGCGCCCGCCGAACCATATTCGAACCATATCAAACATCCCCCCATGGGGAGCATATTTTATTGAGGAtt                                                                                >  1:2225267/1‑148 (MQ=2)
                                                              tggtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACg                                                                                                                                                   <  1:2216120/80‑1 (MQ=34)
                                                              tggtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACg                                                                                                                                                   >  2:2216120/1‑80 (MQ=34)
                                                                gtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGc                                                                             <  1:2138750/148‑1 (MQ=40)
                                                                 tggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCt                                                                            <  1:2469760/148‑1 (MQ=40)
                                                                     gggggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTgag                                                                        <  1:3033179/148‑1 (MQ=255)
                                                                     gggcgAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTgag                                                                        <  1:3033033/148‑1 (MQ=255)
                                                                           aGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCa                                                                  <  1:1162086/148‑1 (MQ=255)
                                                                           aGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCa                                                                  <  1:3412995/148‑1 (MQ=255)
                                                                            ggATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAg                                                                 <  1:1692512/148‑1 (MQ=255)
                                                                            ggATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACCCTCTCTCGCGGGCGGGCGAGCCGCACCGTAGTGGAAGCGTCTCATCCGTCCCCGCCTGGGCCGGATAGTTGATTGGGGGTTCGCGTGTGGGTTCAg                                                                 >  2:1850146/1‑148 (MQ=2)
                                                                             gATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAgg                                                                <  1:510221/148‑1 (MQ=255)
                                                                                tACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCGCTCTCGCTGACGCTCGAGCCGAAACGTAGTCGAAGCTTCTCATCCTTCCCCCCAGGGGCAGAATATTTTATTTCGGATTCGGTTGGGAGTTGAGGGag                                                             >  2:754812/1‑147 (MQ=14)
                                                                                 aCTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACt                                                            <  1:271298/148‑1 (MQ=255)
                                                                                  cTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACtt                                                           <  1:2083782/148‑1 (MQ=255)
                                                                                      gcgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGa                                                       <  1:2643292/148‑1 (MQ=255)
                                                                                      gcgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGa                                                       <  1:2643349/148‑1 (MQ=255)
                                                                                      gcgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGa                                                       <  1:2643446/148‑1 (MQ=255)
                                                                                         cTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTTGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAg                                                    <  1:409239/148‑1 (MQ=255)
                                                                                         cTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAg                                                    <  1:409339/148‑1 (MQ=255)
                                                                                         cTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAg                                                    <  1:409281/148‑1 (MQ=255)
                                                                                           gcgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTg                                                  <  1:2848272/148‑1 (MQ=255)
                                                                                                  tCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCCAGGCGAACCTTAGTCGAAGCTTCTCAACCATCCCCGCATGGGCCGAATATTTGTTTGCGGATTCGGTTGAGATTTAAGGGGCTTTTGAAAGTGTtggtgg                                           >  2:1938269/1‑148 (MQ=11)
                                                                                                    gCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGATTGAGAGTTCAGGGACTTTTGAAAGTGAtggtggtg                                         <  1:1089965/148‑1 (MQ=255)
                                                                                                      ccTTCGGGTCGTTGCCTGCGGCAACCCACTCTCGCTGACGCGCGAGCCGAACCATAGTAGAATCGTATAAACCATACCCGCCTGGGCCGAATATTTGATTTCGGATTAGCTTGTGAGTTGAGGGCCTTTTGAAAGTGAtggtgggtgg                                       >  2:1557271/1‑144 (MQ=9)
                                                                                                        ttCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTggggg                                     <  1:1192565/148‑1 (MQ=255)
                                                                                                        ttCGGGTCGTTGCCGGCGGCAACGCTCACGCGCTGACGCTCGAGCCGAACCATATTCGAAGCTTCTCATCCATCCCCGCCTGGGCCGAATAATTGTTTGCGGATTCGGTTGATAGTTAAGGGACTTTTGTAAGTTATGGTGGTggggg                                     >  2:3402701/1‑148 (MQ=255)
                                                                                                         tCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGa                                    <  1:2593143/148‑1 (MQ=255)
                                                                                                         tCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGa                                    <  1:2593108/148‑1 (MQ=255)
                                                                                                          cGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGaa                                   <  1:2783225/148‑1 (MQ=255)
                                                                                                            ggTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTTGTCGAAGCTTCTCAACCATCCCCCCCTGGGCAGAATATTTGATTGCGGTTTAGCTGGGGAGTTGAGGGACTTTTGTAAGTGATGGTGGTGGGGGGAgg                                 >  2:3495462/1‑148 (MQ=11)
                                                                                                            ggTCGTTGCCTGCGGCAACGCTCACTCGCTGACGCTCGAGCCCAACCTTAGTCGAAGCTTGTCATACTTACCCCCAAGGGCACAATATTTTTTTGCGGATAGCGTTGGTCGTTATGGGCCTTTTTAAAGTGTTGGGGGTGGGGGAAgc                                 >  2:3495554/1‑147 (MQ=2)
                                                                                                              tCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAt                               <  1:1938269/148‑1 (MQ=255)
                                                                                                              tCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAt                               <  1:1938262/148‑1 (MQ=255)
                                                                                                                  tGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGGTGGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACt                           >  2:2077500/1‑148 (MQ=255)
                                                                                                                   gCCTGCGGCAACGCTCTCTCGCTGGCGCTCTAGACGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTc                          <  1:1593872/148‑1 (MQ=255)
                                                                                                                   gCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAttcgaa                          <  1:1557271/148‑5 (MQ=255)
                                                                                                                       gCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgcg                      <  2:3210205/148‑1 (MQ=255)
                                                                                                                       gCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCAAGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgcg                      <  2:3208083/148‑1 (MQ=255)
                                                                                                                        cGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACGTTATTCGAAGCTTCTCATCCTTCCCGGCTTTGTCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCTTTTACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgcgc                     <  1:2967287/148‑1 (MQ=14)
                                                                                                                             aCGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTgcgc                >  1:536606/1‑148 (MQ=255)
                                                                                                                              cGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTgcgcg               <  2:3304151/148‑1 (MQ=255)
                                                                                                                                 tctctcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTGGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGTACGGAATAGGCAGCTCTgcgcgcat            >  2:1054735/1‑146 (MQ=17)
                                                                                                                                      cGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCCGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGGGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGgcc       >  1:5599/1‑148 (MQ=255)
                                                                                                                                      cGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTAGTCAGAGCTGCGAGCTTCGcgc       >  2:472785/1‑146 (MQ=255)
                                                                                                                                         tGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGTGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTc    >  1:3400182/1‑148 (MQ=255)
                                                                                                                                           aCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGGTGGGGAGTTCAGAGAGGTATGTCGGGGTGGGGGGTGCGGGAAGGAGTACTCGGCGCTGCGCGCGGGGCCCTTCgg  >  2:905492/1‑148 (MQ=9)
                                                                                                                                           aCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCgg  >  1:2470795/1‑148 (MQ=255)
                                                                                                                                           |                                                                                                                                                   
CCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGG  >  CP009273/1282282‑1282568

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: