Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,280,478 T→G 5.4% I25M (ATT→ATG nrfD → formate‑dependent nitrite reductase, membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,280,4780TG 5.4% 1203.1 / 46.0 342I25M (ATT→ATGnrfDformate‑dependent nitrite reductase, membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (174/148);  new base G (7/12);  total (182/160)
Fisher's exact test for biased strand distribution p-value = 1.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 349 total aligned reads displayed.

TACCGCTACAAGCTGGCGCTGGGAACCAAACCGAAGCTGTACCGCGTACCGTTTAAATACGGGGAGGTGAGCCAATGACGCAGACTTCCGCATTTCATTTTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACCACGCTGATTGTCGGGCCGGGCGCGGTGATCCTCGGTCTGTTGATCCTCGTCTTCCACC  >  CP009273/4280330‑4280626
                                                                                                                                                    |                                                                                                                                                    
tACCGCTACAAGCTGGCGCTGGGAACCAAACCGAAGCTGTACCGCGTACCGTTTAAATACGGGGAGGTGAGCCAATGACGCAGACTTCCGCATTTCATTTTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGAtt                                                                                                                                                      <  1:4604025/149‑1 (MQ=255)
tACCGCTACAAGCTGGCGCTGGGAACCAAACCGAAGCTGTACCGCGTACCGTTTAAATACGGGGAGGTGAGCCAATGACGCAGACTTCCGCATTTCATTTTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGAtt                                                                                                                                                      >  1:3550499/1‑149 (MQ=255)
   cGCTACAAGCTGGCGCTGGGAACCAAACCGAAGCTGTACCGCGTACCGTTTAAATACGGGGAGGTGAGCCAATGACGCAGACTTCCGCATTTCATTTTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGt                                                                                                                                                   <  1:215151/149‑1 (MQ=255)
    gCTACAAGCTGGCGCTGGGAACCAAACCGAAGCTGTACCGCGTACCGTTTAAATACGGGGAGGTGAGCCAATGACGCAGACTTCCGCATTTCATTTTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTa                                                                                                                                                  >  2:2168456/1‑149 (MQ=255)
        cAAGCTGGCGCTGGGAACCAAACCGAAGCTGTACCGCGTACCGTTTAAATACGGGGAGGTGAGCCAATGACGCAGACTTCCGCATTTCATTTTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTc                                                                                                                                              <  2:699551/149‑1 (MQ=255)
          aGCTGGCGCTGGGAACCAAACCGAAGCTGTACCGCGTACCGTTTAAATACGGGGAGGTGAGCCAATGACGCAGACTTCCGCCTTTCATTTTGAATCGCTGGTGTGGGACTGGCCGAGTGCCATCTACCTGTTTTTGATTGGTATTTctg                                                                                                                                            >  2:323852/1‑149 (MQ=255)
          aGCTGGCGCTGGGAACCAAACCGAAGCTGTACCGCGTACCGTTTAAATACGGGGAGGTGAGCCAATGACGCAGACTTCCGCATTTCATTTTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTctg                                                                                                                                            >  2:2144119/1‑149 (MQ=255)
              ggCGCTGGGAACCAAACCGAAGCTGTACCGCGTACCGTTTAAATACGGGGAGGTGAGCCAATGACGCAGACTTCCGCATTTCATTTTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGctgg                                                                                                                                        <  2:3900200/149‑1 (MQ=255)
                 gcTGGGAACCAACCCGAAGCTGTACCGCGTACCGTTTAAATACGGGGAGGTGAGCCAATGACGCAGACTTCCGCATTTCATTTTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGctggtct                                                                                                                                     <  2:981428/149‑1 (MQ=255)
                 gcTGGGAACCAAACCGAAGCTGTACCGCGTACCGTTTAAATACGGGGAGGTGAGCCAATGACGCAGACTTCCGCATTTCATTTTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGctggtct                                                                                                                                     >  1:2951653/1‑149 (MQ=255)
                 gcTGGGAACCAAACCGAAGCTGTACCGCGTACCGTTTAAATACGGGGAGGTGAGCCAATGACGCAGACTTCCGCATTTCATTTTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGctggtct                                                                                                                                     <  2:5140034/149‑1 (MQ=255)
                    gggAACCAAACCGAAGCTGTACCGCGTACCGTTTAAATACGGGGAGGTGAGCCAATGACGCAGACTTCCGCATTTCATTTTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGctggtctggt                                                                                                                                  >  1:4724614/1‑149 (MQ=255)
                     ggAACCAAACCGAAGCTGTACCGCGTACCGTTTAAATACGGGGAGGTGAGCCAATGACGCAGACTTCCGCATTTCATTTTGAATCGCTGGTGTGGGACTGGCGGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTg                                                                                                                                 <  1:4819735/149‑1 (MQ=255)
                      gAACCAAACCGAAGCTGTACCGCGTACCGTTTAAATACGGGGAGGTGAGCCAATGACGCAGACTTCCGCATTTCATTTTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGa                                                                                                                                <  2:1746283/149‑1 (MQ=255)
                       aaCCAAACCGAAGCTGTACCGCGTACCGTTTAATTACGGGGAGGTGAGCCAATGACGCAGACTTCCGCATTTCATTTTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGAc                                                                                                                               >  2:448391/1‑149 (MQ=255)
                       aaCCAAACCGAAGCTGTACCGCGTACCGTTTAAATACGGGGAGGTGAGCCAATGACGCAGACTTCCGCATTTCATTTTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGAc                                                                                                                               >  1:3284696/1‑149 (MQ=255)
                         ccAAACCGAAGCTGTACCGCGTACCGTTTAAATACGGGGAGGTGAGCCAATGACGCAGACTTCCGCATTTCATTTTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGc                                                                                                                             >  1:961748/1‑149 (MQ=255)
                          cAAACCGAAGCTGTACCGCGTACCGTTTAAATACGGGGAGGTGAGCCAATGACGCAGACTTCCGCATTTCATTTTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCt                                                                                                                            <  1:1478817/149‑1 (MQ=255)
                           aaaCCGAAGCTGTACCGCGTACCGTTTAAATACGGGGAGGTGAGCCAATGACGCAGACTTCCGCATTTCATTTTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTg                                                                                                                           >  2:4304284/1‑149 (MQ=255)
                           aaaCCGAAGCTGTACCGCGTACCGTTTAAATACGGGGAGGTGAGCCAATGACGCAGACTTCCGCATTTCATTTTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTg                                                                                                                           >  2:3205324/1‑149 (MQ=255)
                           aaaCCGAAGCTGTACCGCGTACCGTTTAAATACGGGGAGGTGAGCCAATGACGCAGACTTCCGCATTTCATTTTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTGTTTGAGTGGTATTTCTGCTGGTCTGGTGACGCTg                                                                                                                           >  2:3614968/1‑149 (MQ=255)
                           aaaCCGAAGCTGTACCGCGTACCGTTTAAATACGGGGAGGTGAGCCAATGACGCAGACTTCCGCATTTCATTTTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCGGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTg                                                                                                                           >  1:3465443/1‑149 (MQ=255)
                            aaCCGAAGCTGTACCGCGTACCGTTTAAATACGGGGAGGTGAGCCAATGACGCAGACTTCCGCATTTCATTTTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTgg                                                                                                                          <  1:2205056/149‑1 (MQ=255)
                                 aaGCTGTACCGCGTACCGTTTAAATACGGGGAGGTGAGCCAATGACGCAGACTTCCGCATTTCATTTTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTg                                                                                                                     >  1:1989813/1‑149 (MQ=255)
                                   gCTGTACCGCGTACCGTTTAAATCCGGGGAGGTGAGCCAATGACGCAGACTTCCGCATTTCATTTTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCt                                                                                                                   <  1:5114180/149‑1 (MQ=255)
                                      gTACCGCGTACCGTTTAAATACGGGGAGGTGAGCCAATGACGCAGACTTCCGCATTTCATTTTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGtt                                                                                                                >  1:1836141/1‑149 (MQ=255)
                                       tACCGCGTACCGTTTAAATACGGGGAGGTGAGCCAATGACGCAGACTTCCGCATTTCATTTTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTa                                                                                                               <  2:2883929/149‑1 (MQ=255)
                                          cgcgTACCGTTTAAATACGGGGAGGTGAGCCAATGACGCAGACTTCCGCATTTCATTTTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGt                                                                                                            <  2:620328/149‑1 (MQ=255)
                                             gTACCGTTTAAATACGGGGAGGTGAGCCAATGACGCAGACTTCCGCATTTCATTTTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGc                                                                                                         <  2:3371388/149‑1 (MQ=255)
                                             gTACCGTTTAAATACGGGGAGGTGAGCCAATGACGCAGACTTCCGCATTTCATTTTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGc                                                                                                         >  2:3046553/1‑149 (MQ=255)
                                              tACCGTTTAAATACGGGGAGGTGAGCCAATGACGCAGACTTCCGCATTTCATTTTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCt                                                                                                        <  1:102391/149‑1 (MQ=255)
                                               aCCGTTTAAATACGGGGAGGTGAGCCAATGACGCAGACTTCCGCATTTCATTTTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCtt                                                                                                       >  1:3201576/1‑149 (MQ=255)
                                                ccGTTTAAATACGGGGAGGTGAGCCAATGACGCAGACGTCCGCATTTCATTGTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATGGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTc                                                                                                      <  2:4660787/149‑1 (MQ=255)
                                                  gTTTAAATACGGGGAGGTGAGCCAATGACGCAGACTTCCGCATTTCATTTTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTa                                                                                                    >  1:3604475/1‑149 (MQ=255)
                                                    ttAAATACGGGGAGGTGAGCCAATGACGCAGACTTCCGCATTTCATTTTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATGGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTAcc                                                                                                  >  2:772175/1‑149 (MQ=255)
                                                      aaaTACGGGGAGGTGAGCCAATGACGCAGACTTCCGCATTTCATTTTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCg                                                                                                >  2:3100071/1‑149 (MQ=255)
                                                        aTACGGGGAGGTGAGCCAATGACGCAGACTTCCGCATTTCATTTTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCa                                                                                              <  2:853408/149‑1 (MQ=255)
                                                         tACGGGGAGGTGAGCCAATGACGCAGACTTCCGCATTTCATTTTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAg                                                                                             <  2:4628460/149‑1 (MQ=255)
                                                         tACGGGGAGGTGAGCCAATGACGCAGACTTCCGCATTTCATTTTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAg                                                                                             <  2:4193665/149‑1 (MQ=255)
                                                                    gAGCCAATGACGCAGACTTCCGCATTTCATTTTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGc                                                                                  <  1:576560/149‑1 (MQ=255)
                                                                               gCAGACTTCCGCATTTCATTTTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGt                                                                       <  1:2835388/149‑1 (MQ=255)
                                                                                  gACTTCCGCATTTCATTTTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATGGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTtgc                                                                    <  2:2706234/149‑1 (MQ=255)
                                                                                    cTTCCGCATTTCATTTTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTtgctg                                                                  >  2:4747332/1‑149 (MQ=255)
                                                                                     ttCCGCATTTCATTTTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTtgctgc                                                                 <  2:3135436/149‑1 (MQ=255)
                                                                                     ttCCGCATTTCATTTTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTtgctgc                                                                 >  1:4094682/1‑149 (MQ=255)
                                                                                     ttCCGCATTTCATTTTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTtgctgc                                                                 >  2:159907/1‑149 (MQ=255)
                                                                                        cGCATTTCATTTTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATGGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGca                                                              <  2:4450387/149‑1 (MQ=255)
                                                                                           atttcatttTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGcacca                                                           >  2:676037/1‑149 (MQ=255)
                                                                                           atttcatttTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGcacca                                                           >  2:2498007/1‑149 (MQ=255)
                                                                                            tttcatttTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGcaccac                                                          >  1:3416446/1‑149 (MQ=255)
                                                                                              tcatttTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACCACGc                                                        <  2:3275719/149‑1 (MQ=255)
                                                                                              tcatttTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACAACGc                                                        <  1:1602866/149‑1 (MQ=255)
                                                                                               catttTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACCACGCt                                                       >  1:3805092/1‑149 (MQ=255)
                                                                                               catttTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACCACGCt                                                       >  2:4073968/1‑149 (MQ=255)
                                                                                               catttTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACCACGCt                                                       >  2:3453598/1‑149 (MQ=255)
                                                                                               catttTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACCACGCt                                                       <  2:4767312/149‑1 (MQ=255)
                                                                                                 ttttGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACCACGCTGa                                                     >  2:4956348/1‑149 (MQ=255)
                                                                                                  tttGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACCACGCTGAt                                                    <  2:5075731/149‑1 (MQ=255)
                                                                                                     gactCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGGATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACCACGCTGATTGt                                                 <  2:3177021/146‑1 (MQ=255)
                                                                                                      aaTCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACGGTACGTTGCTGCGCACCACGCTGATTGTc                                                >  1:649589/1‑149 (MQ=255)
                                                                                                      aaTCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACCACGCTGATTGTc                                                >  1:2667371/1‑149 (MQ=255)
                                                                                                      aaTCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACCACGCTGATTGTc                                                <  1:4929471/149‑1 (MQ=255)
                                                                                                       aTCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACCACGCTGATTGTcg                                               <  1:2745287/149‑1 (MQ=255)
                                                                                                            tGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACCACGCTGATTGTcgggccg                                          >  2:3659070/1‑149 (MQ=255)
                                                                                                             ggTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACCACGCTGATTGTcgggccgg                                         <  1:4355382/149‑1 (MQ=255)
                                                                                                             ggTGTGGGACTGGCCGATTGCCATCTACCTGTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACCACGCTGATTGTcgggccggg                                        <  1:333382/149‑1 (MQ=255)
                                                                                                             ggTGTGGGACTAGCCGATTGCCAACAACCTGTAGTTGAGAGGGATGTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGGCGCTTCTACCCGCACGCGGGCGGTGCCGACAGTACGTTGCTGCGCACCACGCTGATACTcgggccgg                                         >  2:1171309/1‑149 (MQ=255)
                                                                                                               tgtgGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACCACGCTGATTGTcgggccgggc                                       >  2:849070/1‑149 (MQ=255)
                                                                                                                gtgGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACCACGCTGATTGTCGGGCCGGgcg                                      <  1:1114061/149‑1 (MQ=255)
                                                                                                                 tgGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACCACGCTGATTGTCGGGCCGGgcgc                                     >  1:2789291/1‑149 (MQ=255)
                                                                                                                  gggACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACCACGCTGATTGTCGGGCCGGGcgcg                                    >  2:3994771/1‑149 (MQ=255)
                                                                                                                   ggACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACCACGCTGATTGTCGGGCCGGGCGCgg                                   >  2:975345/1‑149 (MQ=255)
                                                                                                                     aCTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACCACGCTGATTGTCGGGCCGGGCGCGGTg                                 <  1:2436076/149‑1 (MQ=255)
                                                                                                                      cTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACCACGCTGATTGTCGGGCCGGGCGCGGTGa                                <  2:2842484/149‑1 (MQ=255)
                                                                                                                       tGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACCACGCTGATTGTCGGGCCGGGCGCGGTGAt                               <  1:2109614/149‑1 (MQ=255)
                                                                                                                        ggCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACCACGCTGATTGTCGGGCCGGGCGCGGTGATc                              >  2:2938763/1‑149 (MQ=255)
                                                                                                                        ggCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACCACGCTGATTGTCGGGCCGGGCGCGGTGATc                              >  2:4972756/1‑149 (MQ=255)
                                                                                                                          ccGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACCACGCTGATTGTCGGGCCGGGCGCGGTGATCCt                            >  1:140674/1‑149 (MQ=255)
                                                                                                                          ccGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACCACGCTGATTGTCGGGCCGGGCGCGGTGATCCt                            <  2:3547618/149‑1 (MQ=255)
                                                                                                                          ccGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACCACGCTGATTGTCGGGCCGGGCGCGGTGATCCt                            >  2:2329383/1‑149 (MQ=255)
                                                                                                                           cGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACCACGCTGATTGTCGGGCCGGGCGCGGTGATCCTc                           <  2:1459283/149‑1 (MQ=255)
                                                                                                                             aTTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACCACGCTGATTGTCGGGCCGGGCGCGGTGATCCTCgg                         >  1:901538/1‑149 (MQ=255)
                                                                                                                                gCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACCACGCTGATTGTCGGGCCGGGCGCGGTGATCCTCGGTCt                      <  1:2283899/149‑1 (MQ=255)
                                                                                                                                gCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACCACGCTGATTGTCGGGCCGGGCGCGGTGATCCTCGGTCt                      <  2:2071467/149‑1 (MQ=255)
                                                                                                                                 ccATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACCACGCTGATTGTCGGGCCGGGCGCGGTGATCCTCGGTCTg                     <  1:2888274/149‑1 (MQ=255)
                                                                                                                                  cATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACCACGCTGATTGTCGGGCCGGGCGCGGTGATCCTCGGTCTGt                    <  1:4023134/149‑1 (MQ=255)
                                                                                                                                  cATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACCACGCTGATTGTCGGGCCGGGCGCGGTGATCCTCGGTCTGt                    >  2:1240022/1‑149 (MQ=255)
                                                                                                                                   aTCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACCACGCTGATTGTCGGGCCGGGCGCGGTGATCCTCGGTCTGtt                   >  1:3036732/1‑149 (MQ=255)
                                                                                                                                   aTCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACCACGCTGATTGTCGGGCCGGGCGCGGTGATCCTCGGTCTGtt                   >  1:4287509/1‑149 (MQ=255)
                                                                                                                                   aTCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACCACGCTGATTGTCGGGCCGGGCGCGGTGATCCTCGGTCTGtt                   <  2:3795626/149‑1 (MQ=255)
                                                                                                                                    tCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACCACGCTGATTGTCGGGCCGGGCGCGGTGATCCTCGGTCTGTTg                  >  1:5017663/1‑149 (MQ=255)
                                                                                                                                       aCCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACCACGCTGATTGTCGGGCCGGGCGCGGTGATCCTCGGTCTGTTGATc               >  1:3860536/1‑149 (MQ=255)
                                                                                                                                         cTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACCACGCTGATTGTCGGGCCGGGCGCGGTGATCCTCGGTCTGTTGATCCt             <  1:1263054/149‑1 (MQ=255)
                                                                                                                                           gTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACCACGCTGATTGTCGGGCCGGGCGCGGTGATCCTCGGTCTGTTGATCCTCg           >  1:3703243/1‑149 (MQ=255)
                                                                                                                                           gTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACCACGCTGATTGTCGGGCCGGGCGCGGTGATCCTCGGTCTGTTGATCCTCg           >  2:1898775/1‑149 (MQ=255)
                                                                                                                                            tttttGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACCACGCTGATTGTCGGGCCGGGCGCGGTGATCCTCGGTCTGTTGATCCTCGt          <  1:793775/149‑1 (MQ=255)
                                                                                                                                             ttttGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACCACGCTGATTGTCGGGCCGGGCGCGGTGATCCTCGGTCTGTTGATCCTCGTc         <  1:3498541/149‑1 (MQ=255)
                                                                                                                                                 gATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACCACGCTGATTGTCGGGCCGGGCGCGGTGATCCTCGGTCTGTTGATCCTCGTCTTcc     <  1:2246292/149‑1 (MQ=255)
                                                                                                                                                 gATTGGGATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACCACGCTGATTGTCGGGCCGGGCGCGGTGATCCTCGGTCTGTTGATCCTCGTCTTcc     <  1:4556706/149‑1 (MQ=255)
                                                                                                                                                  aTTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACCACGCTGATTGTCGGGCCGGGCGCGGTGATCCTCGGTCTGTTGATCCTCGTCTTCCa    >  1:1746131/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TACCGCTACAAGCTGGCGCTGGGAACCAAACCGAAGCTGTACCGCGTACCGTTTAAATACGGGGAGGTGAGCCAATGACGCAGACTTCCGCATTTCATTTTGAATCGCTGGTGTGGGACTGGCCGATTGCCATCTACCTGTTTTTGATTGGTATTTCTGCTGGTCTGGTGACGCTGGCCGTGCTGTTACGTCGCTTCTACCCGCAGGCGGGCGGTGCAGACAGTACGTTGCTGCGCACCACGCTGATTGTCGGGCCGGGCGCGGTGATCCTCGGTCTGTTGATCCTCGTCTTCCACC  >  CP009273/4280330‑4280626

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