Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,189,766 C→A 46.4% 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,7660CA46.4% 199.9 / inf 293T112K (ACG→AAG) yjaGDUF416 domain protein
Reads supporting (aligned to +/- strand):  ref base C (61/96);  new base A (75/61);  total (136/157)
Fisher's exact test for biased strand distribution p-value = 6.83e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 299 total aligned reads displayed.

AGACGCAAAAGTAAATTTCGACAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATGGGACATTCAGTGGGAAATAT  >  CP009273/4189617‑4189912
                                                                                                                                                     |                                                                                                                                                  
aGACGCAAAAGTAAATTTCGACAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAc                                                                                                                                                    <  2:614201/150‑1 (MQ=255)
    gCAAAAGTAAATTTCGACAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTc                                                                                                                                                <  1:2585020/150‑1 (MQ=255)
      aaaaGTAAATTTCGACAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGa                                                                                                                                              >  1:1553862/1‑150 (MQ=255)
       aaaGTAAATTTCGACAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGaa                                                                                                                                             >  1:3916624/1‑150 (MQ=255)
       aaaGTAAATTTCGACAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGaa                                                                                                                                             >  1:2051622/1‑150 (MQ=255)
        aaGTAAATTTCGACAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAac                                                                                                                                            <  2:478723/150‑1 (MQ=255)
        aaGTAAATTTCGACAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAac                                                                                                                                            >  1:241452/1‑150 (MQ=255)
        aaGTAAATTTCGACAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAac                                                                                                                                            >  2:3374855/1‑150 (MQ=255)
         aGTAAATTTCGACAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAaca                                                                                                                                           >  1:3747202/1‑150 (MQ=255)
         aGTAAATTTCGACAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAaca                                                                                                                                           >  2:1113379/1‑150 (MQ=255)
         aGTAAATTTCGACAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAaca                                                                                                                                           >  1:941755/1‑150 (MQ=255)
          gTAAATTTCGACAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAacac                                                                                                                                          <  2:2418604/150‑1 (MQ=255)
             aaTTTCGACAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACAcgcg                                                                                                                                       >  1:632507/1‑150 (MQ=255)
              aTTTCGACAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCgg                                                                                                                                      <  2:218229/150‑1 (MQ=255)
               tttCGACAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAACACGCGGt                                                                                                                                     >  1:3015357/1‑150 (MQ=255)
                ttCGACAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTg                                                                                                                                    <  1:4145013/150‑1 (MQ=255)
                 tCGACAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTgg                                                                                                                                   <  2:779932/150‑1 (MQ=255)
                   gACAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAACACGCGGTGGaa                                                                                                                                 >  2:4489728/1‑150 (MQ=255)
                    aCAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAg                                                                                                                                >  1:118982/1‑150 (MQ=255)
                     cAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAACACGCGGTGGAAGt                                                                                                                               >  2:4006099/1‑150 (MQ=255)
                       gCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGa                                                                                                                             >  2:2994992/1‑150 (MQ=255)
                         cAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGTCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAACACGCGTTGGAAGTGAGt                                                                                                                           >  1:2424272/1‑150 (MQ=255)
                             tGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGa                                                                                                                       <  1:1796033/150‑1 (MQ=255)
                               gagaAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGAcc                                                                                                                     >  2:3539149/1‑150 (MQ=255)
                                agaAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCt                                                                                                                    <  2:1007335/150‑1 (MQ=255)
                                  aaaTTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCTcc                                                                                                                  <  2:1301075/150‑1 (MQ=255)
                                    aTTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCAt                                                                                                                <  1:1521963/150‑1 (MQ=255)
                                          agaagCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGAc                                                                                                          >  2:3955234/1‑150 (MQ=255)
                                            aagCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCg                                                                                                        <  2:331352/150‑1 (MQ=255)
                                             agCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGt                                                                                                       >  1:4395575/1‑150 (MQ=255)
                                              gCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGtt                                                                                                      <  2:1459277/150‑1 (MQ=255)
                                               cGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTg                                                                                                     >  1:866195/1‑150 (MQ=255)
                                                      ttCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCt                                                                                              >  1:1071508/1‑150 (MQ=255)
                                                       tCAGCCGACGATTTCGAACTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTg                                                                                             <  2:2370992/150‑1 (MQ=255)
                                                        cAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTgg                                                                                            <  2:4490016/150‑1 (MQ=255)
                                                          gCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGaa                                                                                          <  1:3314746/150‑1 (MQ=255)
                                                          gCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGaa                                                                                          >  1:3626821/1‑150 (MQ=255)
                                                          gCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGaa                                                                                          >  2:3678930/1‑150 (MQ=255)
                                                           ccGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGaaa                                                                                         >  2:310043/1‑150 (MQ=255)
                                                           ccGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGaaa                                                                                         >  2:2608503/1‑150 (MQ=255)
                                                            cgacgaTTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAAt                                                                                        >  2:4235305/1‑150 (MQ=255)
                                                             gacgaTTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATg                                                                                       <  2:456676/150‑1 (MQ=255)
                                                              acgaTTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGa                                                                                      <  1:1322661/150‑1 (MQ=255)
                                                                     cGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGc                                                                               >  2:3401813/1‑150 (MQ=255)
                                                                     cGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGc                                                                               >  2:4138491/1‑150 (MQ=255)
                                                                     cGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTGTGAGTGGTGAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGc                                                                               >  2:1981502/1‑150 (MQ=255)
                                                                      gATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCt                                                                              <  1:801069/150‑1 (MQ=255)
                                                                       aTCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTg                                                                             <  1:536235/150‑1 (MQ=255)
                                                                          tGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTc                                                                          <  1:3263906/150‑1 (MQ=255)
                                                                           gTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTcg                                                                         <  2:3495440/150‑1 (MQ=255)
                                                                            tACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTcgc                                                                        <  2:325869/150‑1 (MQ=255)
                                                                            tACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTcgc                                                                        >  2:267003/1‑150 (MQ=255)
                                                                             aCGGCGTGTATCCGGCAATCGATGCCTGCGTGGCGGTAAGTGAACTGGTCCATTCGCGTGTGAGTGGGGAAAAGCTCGAACACGCGGTGGACGTGAGTACGTCCTCCATCACGACCGGTGCGATGCTGGCAATGACTCTGGCGGGGcgcg                                                                       >  2:3680064/1‑150 (MQ=255)
                                                                                  gTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATg                                                                  >  1:2786299/1‑150 (MQ=255)
                                                                                     tATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGc                                                               <  1:978376/150‑1 (MQ=255)
                                                                                     tATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGc                                                               <  2:1417236/150‑1 (MQ=255)
                                                                                     tATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGc                                                               >  1:700245/1‑150 (MQ=255)
                                                                                         cGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATg                                                           >  2:2750206/1‑150 (MQ=255)
                                                                                          ggCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGa                                                          >  1:1372674/1‑150 (MQ=255)
                                                                                            cAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAag                                                        <  2:688157/150‑1 (MQ=255)
                                                                                            cAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAag                                                        <  2:2729843/150‑1 (MQ=255)
                                                                                                cGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCt                                                    <  2:869599/150‑1 (MQ=255)
                                                                                                cGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCt                                                    <  2:2330699/150‑1 (MQ=255)
                                                                                                 gATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCtt                                                   >  2:1634983/1‑150 (MQ=255)
                                                                                                 gATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCtt                                                   >  2:1838468/1‑150 (MQ=255)
                                                                                                 gATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCtt                                                   <  1:2132683/150‑1 (MQ=255)
                                                                                                 gATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCtt                                                   <  1:1313710/150‑1 (MQ=255)
                                                                                                 gATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCtt                                                   <  1:4030522/150‑1 (MQ=255)
                                                                                                   tGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTaa                                                 >  1:2254541/1‑150 (MQ=255)
                                                                                                   tGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTaa                                                 <  1:3478316/150‑1 (MQ=255)
                                                                                                    gCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTaaa                                                <  1:3526873/150‑1 (MQ=255)
                                                                                                    gCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTaaa                                                <  1:1002194/150‑1 (MQ=255)
                                                                                                    gCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTaaa                                                <  1:3280964/150‑1 (MQ=255)
                                                                                                       tGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGaa                                             >  2:2093859/1‑150 (MQ=255)
                                                                                                         cGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGaaaa                                           >  1:3675114/1‑150 (MQ=255)
                                                                                                                 tAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTg                                   <  1:3374855/150‑1 (MQ=255)
                                                                                                                   aGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTa                                 >  1:3216053/1‑150 (MQ=255)
                                                                                                                    gTGAACTGGTCCTTTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGATAACCCAGCTGTag                                <  2:1363806/150‑1 (MQ=255)
                                                                                                                    gTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTag                                <  1:3093775/150‑1 (MQ=255)
                                                                                                                    gTGAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTag                                <  1:3093741/150‑1 (MQ=255)
                                                                                                                      gAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTagag                              <  1:1292441/150‑1 (MQ=255)
                                                                                                                      gAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTagag                              <  2:3525830/150‑1 (MQ=255)
                                                                                                                      gAACTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTagag                              >  2:1670482/1‑150 (MQ=255)
                                                                                                                         cTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCaa                           <  2:4456321/150‑1 (MQ=255)
                                                                                                                         cTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCaa                           >  1:3419647/1‑150 (MQ=255)
                                                                                                                         cTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCaa                           >  1:1146322/1‑150 (MQ=255)
                                                                                                                         cTGGTCCATTCGCGTTTGAGTGGTGAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCaa                           >  2:1779052/1‑150 (MQ=255)
                                                                                                                          tGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAg                          <  1:2550374/150‑1 (MQ=255)
                                                                                                                          tGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAg                          >  1:2217152/1‑150 (MQ=255)
                                                                                                                          tGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAg                          <  1:1849999/150‑1 (MQ=255)
                                                                                                                              ccATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATg                      <  2:3747202/150‑1 (MQ=255)
                                                                                                                                aTTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATggg                    >  2:1190878/1‑150 (MQ=255)
                                                                                                                                   cgcgTTTGAGTGGTGAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATGGGACa                 >  1:3861617/1‑150 (MQ=255)
                                                                                                                                     cgTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATGGGACAtt               <  1:4235305/150‑1 (MQ=255)
                                                                                                                                      gTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATGGGACATTc              >  1:1574951/1‑150 (MQ=255)
                                                                                                                                        ttGAGTGGTGAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCCGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATGGGACATTCAg            <  2:2219062/150‑1 (MQ=255)
                                                                                                                                           aGTGGTGAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATGGGACATTCAGTgg         <  2:962646/150‑1 (MQ=255)
                                                                                                                                              ggtgAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATGGGACATTCAGTGGGaa      >  2:3813585/1‑150 (MQ=255)
                                                                                                                                               gtgAAAAGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATGGGACATTCAGTGGGaaa     <  1:724216/150‑1 (MQ=255)
                                                                                                                                                  aaaCGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATGGGACATTCAGTGGGAAatat  >  1:3340151/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
AGACGCAAAAGTAAATTTCGACAGCCAACTGGAGAAATTTGAAGAAGCGATTCCTTCAGCCGACGATTTCGATCTGTACGGCGTTTATCCGGCAATCGATGCCTGCGTGGCGTTAAGTGAACTGGTCCATTCGCGTTTGAGTGGTGAAACGCTCGAACACGCGGTGGAAGTGAGTAAGACCTCCATCACGACCGTTGCGATGCTGGAAATGACTCAGGCTGGTCGCGAAATGAGCGATGAAGAGCTTAAAGAAAACCCAGCTGTAGAGCAAGAATGGGACATTCAGTGGGAAATAT  >  CP009273/4189617‑4189912

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: