Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,402,850 G→T 89.6% L2L (CTG→CTT)  yjfC → ATP‑Grasp family ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,402,8500GT89.6% 702.3 / 79.7 259L2L (CTG→CTT) yjfCATP‑Grasp family ATPase
Reads supporting (aligned to +/- strand):  ref base G (16/11);  new base T (106/126);  total (122/137)
Fisher's exact test for biased strand distribution p-value = 2.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 260 total aligned reads displayed.

TGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTGAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGAGCAGATCGAAA  >  CP009273/4402704‑4402995
                                                                                                                                                  |                                                                                                                                                 
 ggTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGac                                                                                                                                               <  2:3114351/150‑1 (MQ=255)
  gTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGaca                                                                                                                                              <  2:28138/150‑1 (MQ=255)
   tGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGacac                                                                                                                                             >  2:4087909/1‑150 (MQ=255)
    gATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGAcaca                                                                                                                                            >  2:44477/1‑150 (MQ=255)
      ttACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCAGTGGGGAGGCTAATCATGCTTAGACACAAc                                                                                                                                          >  1:44899/1‑150 (MQ=255)
       tACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACg                                                                                                                                         <  2:3792528/150‑1 (MQ=255)
        aCTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGt                                                                                                                                        >  2:11415/1‑150 (MQ=255)
          tCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTc                                                                                                                                      <  2:2413926/150‑1 (MQ=255)
          tCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTGAGACACAACGTTc                                                                                                                                      <  2:511467/150‑1 (MQ=255)
            cTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCt                                                                                                                                    >  1:3187888/1‑150 (MQ=255)
            cTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCt                                                                                                                                    >  2:2894544/1‑150 (MQ=255)
              ggCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCtgt                                                                                                                                  <  1:4071677/150‑1 (MQ=255)
                cgcTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGc                                                                                                                                <  1:3606716/150‑1 (MQ=255)
                   tCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGAc                                                                                                                             >  1:132714/1‑150 (MQ=255)
                     cGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACgg                                                                                                                           <  2:3683862/150‑1 (MQ=255)
                        ttCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGAt                                                                                                                        <  2:2037407/150‑1 (MQ=255)
                        ttCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGAt                                                                                                                        >  2:2715246/1‑150 (MQ=255)
                         tCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATc                                                                                                                       >  2:970788/1‑150 (MQ=255)
                          cTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCt                                                                                                                      >  1:3968866/1‑150 (MQ=255)
                             gCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGa                                                                                                                   >  2:2549723/1‑150 (MQ=255)
                             gCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGa                                                                                                                   <  1:3759449/150‑1 (MQ=255)
                              cAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTGAGACACAACGTTCCTGTGCGACGGGATCTGGAc                                                                                                                  >  1:2281994/1‑150 (MQ=255)
                               aaaaaaGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGAcc                                                                                                                 <  2:683679/150‑1 (MQ=255)
                                aaaaaGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTGAGACACAACGTTCCTGTGCGACGGGATCTGGACCa                                                                                                                >  2:4063308/1‑150 (MQ=255)
                                aaaaaGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATGGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCa                                                                                                                >  2:3383095/1‑150 (MQ=255)
                                 aaaaGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAg                                                                                                               >  1:1235425/1‑150 (MQ=255)
                                 aaaaGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTGAGACACAACGTTCCTGTGCGACGGGATCTGGACCAg                                                                                                               >  2:3148962/1‑150 (MQ=255)
                                   aaGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGAt                                                                                                             <  1:2535785/150‑1 (MQ=255)
                                   aaGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGAt                                                                                                             >  1:3631348/1‑150 (MQ=255)
                                     gagTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATcg                                                                                                           <  2:3704004/150‑1 (MQ=255)
                                     gagTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATcg                                                                                                           >  1:3642618/1‑150 (MQ=255)
                                       gTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATcgcc                                                                                                         <  1:867478/150‑1 (MQ=255)
                                       gTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATcgcc                                                                                                         <  1:541900/150‑1 (MQ=255)
                                           tACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTg                                                                                                     <  1:2030063/150‑1 (MQ=255)
                                             ctctTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGAc                                                                                                   <  1:1731936/150‑1 (MQ=255)
                                              tctTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACa                                                                                                  <  2:869512/150‑1 (MQ=255)
                                              tctTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCAGTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACa                                                                                                  >  1:812585/1‑150 (MQ=255)
                                               ctTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACaa                                                                                                 >  1:1919854/1‑150 (MQ=255)
                                                ttCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATTCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAAc                                                                                                >  1:1215583/1‑150 (MQ=255)
                                                ttCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTGAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAAc                                                                                                <  1:3783506/150‑1 (MQ=255)
                                                  cGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACgg                                                                                              >  2:3569186/1‑150 (MQ=255)
                                                  cGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACgg                                                                                              <  1:1407219/150‑1 (MQ=255)
                                                  cGGGCGGCTTTACCACCAAAAAAGCGTCAACCGTTACTCGCGGCTGCGATGTTCGTTCTTCCAGCGCCCGTGGGCATAGGGGAGTCTAATCAAGCGTAGACACAACGTTCCTGAGCGACGGGATCTGGACCAGATCGCCGCTGACAAccg                                                                                              >  2:2223699/1‑148 (MQ=255)
                                                   gggCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTGAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGt                                                                                             >  1:4205332/1‑150 (MQ=255)
                                                     gcggcTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGttt                                                                                           <  1:3211068/150‑1 (MQ=255)
                                                      cggcTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTc                                                                                          <  1:3911039/150‑1 (MQ=255)
                                                      cggcTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTc                                                                                          >  2:1391174/1‑150 (MQ=255)
                                                      cggcTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTc                                                                                          <  1:2349376/150‑1 (MQ=255)
                                                      cggcTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTc                                                                                          <  2:3892503/150‑1 (MQ=255)
                                                      cggcTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTc                                                                                          <  2:890388/150‑1 (MQ=255)
                                                         cTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGAc                                                                                       <  2:3815422/150‑1 (MQ=255)
                                                          tttACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACACCGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACt                                                                                      <  2:1160430/150‑1 (MQ=255)
                                                           ttACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACtt                                                                                     <  1:4609270/150‑1 (MQ=255)
                                                                accaAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACTTTCata                                                                                >  2:2931771/1‑150 (MQ=255)
                                                                  caAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACTTTCATatc                                                                              >  2:2485569/1‑150 (MQ=255)
                                                                   aaaaaaGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACTTTCATatca                                                                             >  2:2682722/1‑150 (MQ=255)
                                                                   aaaaaaGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACTTTCATatca                                                                             <  1:3899037/150‑1 (MQ=255)
                                                                    aaaaaGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTGAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACTTTCATatcat                                                                            >  1:292068/1‑150 (MQ=255)
                                                                     aaaaGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACTTTCATatcatc                                                                           <  1:149484/150‑1 (MQ=255)
                                                                     aaaaGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTGAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACTTTCATatcatc                                                                           <  1:2692973/150‑1 (MQ=255)
                                                                       aaGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACTTTCATATCATCGa                                                                         >  2:620182/1‑150 (MQ=255)
                                                                       aaGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTGAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACTTTCATATCATCGa                                                                         <  1:3861117/150‑1 (MQ=255)
                                                                          cGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGATCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACTTTAATATCATCGACaa                                                                      >  1:2760664/1‑150 (MQ=255)
                                                                           gTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACTTTCATATCATCGACAAt                                                                     <  1:2077091/150‑1 (MQ=255)
                                                                                ccGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTGAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACTTTCATATCATCGACAATGAAAt                                                                >  2:4125571/1‑150 (MQ=255)
                                                                                 cGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACTTTCATATCATCGACAATGAAATc                                                               <  2:1042361/150‑1 (MQ=255)
                                                                                   tttCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACTTTCATATCATCGACAATGAAATCTa                                                             <  2:3097265/150‑1 (MQ=255)
                                                                                     tctcGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACTTTCATATCATCGACAATGAAATCTAtt                                                           <  2:4079996/150‑1 (MQ=255)
                                                                                      ctcGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACGTTCATATCATCGACACAGAAATCTAttt                                                          >  1:1927979/1‑149 (MQ=255)
                                                                                              ggcTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACTTTCATATCATCGACAATGAAATCTATTGGGATgaga                                                  >  2:4387667/1‑150 (MQ=255)
                                                                                                 tATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTc                                               <  2:2460029/150‑1 (MQ=255)
                                                                                                 tATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCAAGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTc                                               >  1:1585381/1‑150 (MQ=255)
                                                                                                  aTGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCg                                              >  1:2680964/1‑150 (MQ=255)
                                                                                                   tGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCgg                                             <  1:1012178/150‑1 (MQ=255)
                                                                                                   tGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCgg                                             <  2:3298549/150‑1 (MQ=255)
                                                                                                   tGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCgg                                             <  1:4619202/150‑1 (MQ=255)
                                                                                                      tCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCt                                          <  1:2186172/150‑1 (MQ=255)
                                                                                                        gTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTa                                        <  1:4367835/150‑1 (MQ=255)
                                                                                                            ttcCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGt                                    >  1:1150555/1‑150 (MQ=255)
                                                                                                                aGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTa                                <  1:4407020/150‑1 (MQ=255)
                                                                                                                 gcgcCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTAc                               <  1:910331/150‑1 (MQ=255)
                                                                                                                  cgcCCGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTAct                              >  1:2830321/1‑150 (MQ=255)
                                                                                                                     ccGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTg                           >  1:3386571/1‑150 (MQ=255)
                                                                                                                      cGTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTgc                          >  2:4301194/1‑150 (MQ=255)
                                                                                                                       gTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTgcg                         >  1:2241391/1‑150 (MQ=255)
                                                                                                                       gTGGGCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTgcg                         >  2:1220588/1‑150 (MQ=255)
                                                                                                                       gTGGGCAGTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTgcg                         >  2:3498215/1‑150 (MQ=255)
                                                                                                                         gggCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGcc                       <  2:432942/150‑1 (MQ=255)
                                                                                                                           gCATTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAg                     >  1:3509497/1‑150 (MQ=255)
                                                                                                                             aTTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGAt                   <  2:1520847/150‑1 (MQ=255)
                                                                                                                             aTTGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGAt                   <  2:2696116/150‑1 (MQ=255)
                                                                                                                              ttGGGGAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGAtt                  <  2:417247/150‑1 (MQ=255)
                                                                                                                                ggggAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGa                >  2:1795976/1‑150 (MQ=255)
                                                                                                                                  ggAGGCTAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAag              <  1:540092/150‑1 (MQ=255)
                                                                                                                                        tAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGAGCAGa        >  1:3427832/1‑150 (MQ=255)
                                                                                                                                        tAATCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGAGCAGa        <  1:1026726/150‑1 (MQ=255)
                                                                                                                                         aaTCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGAGCAGAt       <  1:3111362/150‑1 (MQ=255)
                                                                                                                                           tCATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGAGCAGATCg     <  1:4286529/150‑1 (MQ=255)
                                                                                                                                            cATGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGAGCAGATCGa    >  1:65958/1‑150 (MQ=255)
                                                                                                                                              tGCTTAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGAGCAGATCGaaa  <  2:4601763/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGTGGGCATTGGGGAGGCTAATCATGCTGAGACACAACGTTCCTGTGCGACGGGATCTGGACCAGATCGCCGCTGACAACGGTTTCGACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGAGCAGATCGAAA  >  CP009273/4402704‑4402995

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: ‑