Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,189,766 C→A 12.7% T112K (ACG→AAG)  yjaG → DUF416 domain protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,189,7660CA12.7% 1055.2 / 84.6 318T112K (ACG→AAG) yjaGDUF416 domain protein
Reads supporting (aligned to +/- strand):  ref base C (146/131);  new base A (27/14);  total (173/145)
Fisher's exact test for biased strand distribution p-value = 1.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 330 total aligned reads displayed.

GACGCAAAAGTAAATTTCGACAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATGGGACATTCAGTGGGAAATATTC  >  CP009273/4189618‑4189914
                                                                                                                                                    |                                                                                                                                                    
gACGCAAAAGTAAATTTCGACAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAc                                                                                                                                                      <  1:1430318/149‑1 (MQ=255)
gACGCAAAAGTAAATTTCGACAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAc                                                                                                                                                      >  1:4734738/1‑149 (MQ=255)
gACGCAAAAGTAAATTTCGACAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAc                                                                                                                                                      >  1:316817/1‑149 (MQ=255)
  cGCAAAAGTAAATTTCGACAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGc                                                                                                                                                    <  1:4470627/149‑1 (MQ=255)
    cAAAAGTAAATTTCGACAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTc                                                                                                                                                  <  1:766902/149‑1 (MQ=255)
     aaaaGTAAATTTCGACAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCg                                                                                                                                                 >  2:2903673/1‑149 (MQ=255)
       aaGTAAATTTCGACAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGaa                                                                                                                                               >  1:2066141/1‑149 (MQ=255)
         gTAAATTTCGACAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAaca                                                                                                                                             <  2:1504093/149‑1 (MQ=255)
               ttCGACAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGt                                                                                                                                       <  2:2033427/149‑1 (MQ=255)
               ttCGACAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGt                                                                                                                                       <  2:2981292/149‑1 (MQ=255)
                   aCAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGaa                                                                                                                                   >  2:728929/1‑149 (MQ=255)
                    cAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAg                                                                                                                                  >  2:4385/1‑149 (MQ=255)
                     aGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGt                                                                                                                                 >  1:3256354/1‑149 (MQ=255)
                      gCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTg                                                                                                                                <  1:1037973/149‑1 (MQ=255)
                        cAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAg                                                                                                                              <  2:3789122/149‑1 (MQ=255)
                         aaCTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGt                                                                                                                             >  2:3684513/1‑149 (MQ=255)
                         aaCTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAACACGCGGTGGAAGTGAGt                                                                                                                             <  2:351364/149‑1 (MQ=255)
                           cTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTaa                                                                                                                           >  1:1690577/1‑149 (MQ=255)
                            tGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAg                                                                                                                          <  1:4629224/149‑1 (MQ=255)
                            tGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAg                                                                                                                          >  1:3869878/1‑149 (MQ=255)
                             ggAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGa                                                                                                                         >  1:3315878/1‑149 (MQ=255)
                              gagaAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGAc                                                                                                                        >  2:2878998/1‑149 (MQ=255)
                              gagaAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGAc                                                                                                                        >  2:1742485/1‑149 (MQ=255)
                                 aaaTTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTc                                                                                                                     <  2:3351004/149‑1 (MQ=255)
                                       gaagaaGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCAc                                                                                                               >  1:3400944/1‑149 (MQ=255)
                                         agaagCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCCATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGa                                                                                                             >  2:2220436/1‑149 (MQ=255)
                                         agaagCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGa                                                                                                             >  2:3830778/1‑149 (MQ=255)
                                               gATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTg                                                                                                       >  1:4780859/1‑149 (MQ=255)
                                                 ttCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCg                                                                                                     >  2:4690275/1‑149 (MQ=255)
                                                   ccTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGAt                                                                                                   <  2:1288565/149‑1 (MQ=255)
                                                     ttCACCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGc                                                                                                 >  2:537828/1‑149 (MQ=255)
                                                      tCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCt                                                                                                <  2:1838151/149‑1 (MQ=255)
                                                      tCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCt                                                                                                <  2:4227551/149‑1 (MQ=255)
                                                         gCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGa                                                                                             >  1:1629850/1‑149 (MQ=255)
                                                          ccGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGaa                                                                                            >  1:2767181/1‑149 (MQ=255)
                                                           cgacgaTTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGaaa                                                                                           >  1:3605652/1‑149 (MQ=255)
                                                           cgacgaTTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGaaa                                                                                           <  1:1297593/149‑1 (MQ=255)
                                                            gacgaTTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAAt                                                                                          >  2:3441670/1‑149 (MQ=255)
                                                              cgaTTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGa                                                                                        >  1:423281/1‑149 (MQ=255)
                                                              cgaTTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGa                                                                                        <  2:1421627/149‑1 (MQ=255)
                                                               gaTTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGAc                                                                                       >  2:4129007/1‑149 (MQ=255)
                                                                aTTTCGATCTGTACGGCTTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGaat                                                                                      >  2:2304998/1‑147 (MQ=255)
                                                                aTTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACt                                                                                      >  2:394358/1‑149 (MQ=255)
                                                                 tttCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTc                                                                                     <  2:1885728/149‑1 (MQ=255)
                                                                 tttCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTc                                                                                     >  2:16881/1‑149 (MQ=255)
                                                                   tCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAg                                                                                   <  2:3306530/149‑1 (MQ=255)
                                                                    cGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAgg                                                                                  <  2:3265192/149‑1 (MQ=255)
                                                                    cGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAgg                                                                                  <  2:4520540/149‑1 (MQ=255)
                                                                     gATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGc                                                                                 <  1:1845230/149‑1 (MQ=255)
                                                                     gATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGc                                                                                 >  2:628558/1‑149 (MQ=255)
                                                                     gATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGc                                                                                 >  1:1157870/1‑149 (MQ=255)
                                                                       tCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTg                                                                               >  2:1463734/1‑149 (MQ=255)
                                                                       tCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTg                                                                               >  1:230921/1‑149 (MQ=255)
                                                                         tGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGt                                                                             <  1:2600243/149‑1 (MQ=255)
                                                                               gCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGaa                                                                       >  1:585466/1‑149 (MQ=255)
                                                                                  tttATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATg                                                                    >  2:1103491/1‑149 (MQ=255)
                                                                                    tATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAg                                                                  >  2:4304489/1‑149 (MQ=255)
                                                                                     aTCCGGCACTCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGc                                                                 <  2:2178621/149‑1 (MQ=255)
                                                                                        cGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGAt                                                              <  1:3125660/149‑1 (MQ=255)
                                                                                        cGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGAt                                                              <  1:2249681/149‑1 (MQ=255)
                                                                                          gCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGa                                                            >  1:4159188/1‑149 (MQ=255)
                                                                                          gCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGAGATGa                                                            >  2:2123851/1‑149 (MQ=255)
                                                                                             aTCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAaga                                                         >  2:983279/1‑149 (MQ=255)
                                                                                             aTCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTACGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAaga                                                         >  2:1254235/1‑149 (MQ=255)
                                                                                                gATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCt                                                      <  1:1006266/149‑1 (MQ=255)
                                                                                                 aTGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCtt                                                     >  1:868225/1‑149 (MQ=255)
                                                                                                    ccTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTaaa                                                  >  2:2646097/1‑149 (MQ=255)
                                                                                                     cTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAg                                                 >  1:213122/1‑149 (MQ=255)
                                                                                                       gCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGTCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGaa                                               >  2:1066798/1‑149 (MQ=255)
                                                                                                       gCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGaa                                               >  2:3206967/1‑149 (MQ=255)
                                                                                                         gTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGaaaa                                             >  1:1496893/1‑149 (MQ=255)
                                                                                                          tGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAAc                                            >  1:1864510/1‑149 (MQ=255)
                                                                                                          tGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAAc                                            >  2:866767/1‑149 (MQ=255)
                                                                                                          tGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAAc                                            >  1:406901/1‑149 (MQ=255)
                                                                                                            gCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAAccc                                          >  1:2412697/1‑149 (MQ=255)
                                                                                                             agTTAAGTTATCTGGTACATTCGAGTTTGAGTGGTGAAAAGCTCGAACACGAGGTGGAAGTGAGTAAGACCTCAATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCa                                         <  1:2842590/148‑1 (MQ=255)
                                                                                                              gTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAg                                        >  1:2807303/1‑149 (MQ=255)
                                                                                                                tAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCt                                      <  2:3887086/149‑1 (MQ=255)
                                                                                                                  aGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGt                                    >  1:3721401/1‑149 (MQ=255)
                                                                                                                   gTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTa                                   >  2:1623864/1‑149 (MQ=255)
                                                                                                                     gAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTaga                                 <  2:1943050/149‑1 (MQ=255)
                                                                                                                     gAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTaga                                 <  2:4483036/149‑1 (MQ=255)
                                                                                                                       aCTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGc                               <  2:906961/149‑1 (MQ=255)
                                                                                                                            tCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGaa                          <  1:3508817/149‑1 (MQ=255)
                                                                                                                              cATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATg                        >  1:2893880/1‑149 (MQ=255)
                                                                                                                              cATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATg                        <  1:2522186/149‑1 (MQ=255)
                                                                                                                              cATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATg                        <  2:2503861/149‑1 (MQ=255)
                                                                                                                                        tGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCTAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATGGGACATTCAg              <  2:2213036/149‑1 (MQ=255)
                                                                                                                                        tGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATGGGACATTCAg              <  1:4585202/149‑1 (MQ=255)
                                                                                                                                        tGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATGGGACATTCAg              <  1:993965/149‑1 (MQ=255)
                                                                                                                                         gAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATGGGACATTCAGt             >  2:126098/1‑149 (MQ=255)
                                                                                                                                            tggtgAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATGGGACATTCAGTggg          <  2:4734738/149‑1 (MQ=255)
                                                                                                                                            tggtgAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATGGGACATTCAGTggg          <  1:1115677/149‑1 (MQ=255)
                                                                                                                                            tggtgAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATGGGACATTCAGTggg          >  1:3038820/1‑149 (MQ=255)
                                                                                                                                             ggtgAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATGGGACATTCAGTGGGa         >  2:2240266/1‑149 (MQ=255)
                                                                                                                                               tgAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATGGGACATTCAGTGGGaaa       <  1:4385/149‑1 (MQ=255)
                                                                                                                                                  aaCGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATGGGACATTCAGTGGGAAatat    >  1:4343060/1‑149 (MQ=255)
                                                                                                                                                  aaCGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATGGGACATTCAGTGGGAAatat    <  1:175471/149‑1 (MQ=255)
                                                                                                                                                    cGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATGGGACATTCAGTGGGAAATATTc  <  2:4578288/149‑1 (MQ=255)
                                                                                                                                                    cGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATGGGACATTCAGTGGGAAATATTc  <  1:3660258/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GACGCAAAAGTAAATTTCGACAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATGGGACATTCAGTGGGAAATATTC  >  CP009273/4189618‑4189914

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: