Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,322,609 C→A 90.8% W223L (TGG→TTG)  yfaQ ← tandem DUF2300 domain protein, putative host defense protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,322,6090CA90.8% 399.9 / 35.9 153W223L (TGG→TTG) yfaQtandem DUF2300 domain protein, putative host defense protein
Reads supporting (aligned to +/- strand):  ref base C (8/6);  new base A (71/68);  total (79/74)
Fisher's exact test for biased strand distribution p-value = 7.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

Only 100 of 158 total aligned reads displayed.

TGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCAC  >  CP009273/2322464‑2322751
                                                                                                                                                 |                                                                                                                                              
tgtgATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCaac                                                                                                                                              <  1:2241047/148‑1 (MQ=255)
tgtgATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCaac                                                                                                                                              >  1:1212287/1‑148 (MQ=255)
 gtgATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCaaca                                                                                                                                             <  2:1945292/148‑1 (MQ=255)
 gtgATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCaaca                                                                                                                                             <  2:275590/148‑1 (MQ=255)
    aTTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCg                                                                                                                                          <  2:497459/148‑1 (MQ=255)
     ttttCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCgg                                                                                                                                         >  2:3248341/1‑148 (MQ=255)
     ttttCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCgg                                                                                                                                         >  2:1282190/1‑148 (MQ=255)
     ttttCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCgg                                                                                                                                         >  1:2694473/1‑148 (MQ=255)
       ttcttcAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGgtg                                                                                                                                       <  1:1012278/148‑1 (MQ=255)
            cAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCa                                                                                                                                  >  2:277303/1‑148 (MQ=255)
             aGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCAt                                                                                                                                 <  2:3234580/148‑1 (MQ=255)
                ggATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCgg                                                                                                                              <  1:1531896/148‑1 (MQ=255)
                ggATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCgg                                                                                                                              <  2:803210/148‑1 (MQ=255)
                ggATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCgg                                                                                                                              <  1:201091/148‑1 (MQ=255)
                ggATAGCGGGCAAACAGCTCCACCTCTACACAGTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCgg                                                                                                                              <  2:803770/148‑1 (MQ=255)
                  aTAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGtt                                                                                                                            >  2:251738/1‑148 (MQ=255)
                   tAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTa                                                                                                                           <  1:1160023/148‑1 (MQ=255)
                   tAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTa                                                                                                                           <  2:1526718/148‑1 (MQ=255)
                     gCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAAc                                                                                                                         >  2:212087/1‑148 (MQ=255)
                        ggCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAAccag                                                                                                                      >  1:2111606/1‑148 (MQ=255)
                         gCAAACAGCTCCACCTCTACACATTGCCCACGGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAAccagc                                                                                                                     <  1:2989506/148‑1 (MQ=255)
                             aCAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACcagccagc                                                                                                                 >  2:3195472/1‑148 (MQ=255)
                             aCAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACcagccagc                                                                                                                 >  1:337562/1‑148 (MQ=255)
                               aGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCaa                                                                                                               >  1:3280568/1‑148 (MQ=255)
                                   ccACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAacc                                                                                                           >  1:821371/1‑148 (MQ=255)
                                    cACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAacca                                                                                                          <  2:1020212/148‑1 (MQ=255)
                                      ccTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAaccacc                                                                                                        >  2:3190274/1‑148 (MQ=255)
                                      ccTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAaccacc                                                                                                        <  1:41172/148‑1 (MQ=255)
                                       ctctACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCAccc                                                                                                       <  1:1003417/148‑1 (MQ=255)
                                         ctACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTg                                                                                                     <  1:1168866/148‑1 (MQ=255)
                                             acaTTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACgg                                                                                                 <  2:2004830/148‑1 (MQ=255)
                                              caTTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGa                                                                                                <  2:2463818/148‑1 (MQ=255)
                                              caTTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGa                                                                                                <  2:2463678/148‑1 (MQ=255)
                                                ttGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGAt                                                                                               >  1:2210260/1‑147 (MQ=255)
                                                    ccACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCt                                                                                          <  1:3019759/148‑1 (MQ=255)
                                                     cACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATCAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTg                                                                                         <  1:3235482/148‑1 (MQ=255)
                                                        tGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGttt                                                                                      >  2:2855650/1‑148 (MQ=255)
                                                        tGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGttt                                                                                      <  1:704693/148‑1 (MQ=255)
                                                        tGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGttt                                                                                      <  1:704660/148‑1 (MQ=255)
                                                          gCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTat                                                                                    >  2:2290855/1‑148 (MQ=255)
                                                             aCCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTatcat                                                                                 >  1:1021260/1‑148 (MQ=255)
                                                                  cGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCa                                                                            >  1:127374/1‑148 (MQ=255)
                                                                    tCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCaaa                                                                          <  1:3126000/148‑1 (MQ=255)
                                                                    tCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCaaa                                                                          >  2:3098155/1‑148 (MQ=255)
                                                                    tCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCaaa                                                                          <  2:2532446/148‑1 (MQ=255)
                                                                    tCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCaaa                                                                          <  1:302395/148‑1 (MQ=255)
                                                                             gCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCa                                                                 <  1:1043883/148‑1 (MQ=255)
                                                                             gCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCa                                                                 >  1:2870627/1‑148 (MQ=255)
                                                                             gCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCa                                                                 <  1:2696980/148‑1 (MQ=255)
                                                                             gCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCa                                                                 >  1:2870523/1‑148 (MQ=255)
                                                                                aaaCCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGt                                                              <  2:2686014/148‑1 (MQ=255)
                                                                                 aaCCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGtt                                                             >  1:2079303/1‑148 (MQ=255)
                                                                                  aCCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGttt                                                            >  2:559645/1‑148 (MQ=255)
                                                                                  aCCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGttt                                                            >  1:1688066/1‑148 (MQ=255)
                                                                                    ccGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTc                                                          <  1:157051/148‑1 (MQ=255)
                                                                                       aTTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCAcc                                                       >  1:2693536/1‑148 (MQ=255)
                                                                                       aTTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCAcc                                                       >  1:2091782/1‑148 (MQ=255)
                                                                                       aTTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCAcc                                                       <  1:44308/148‑1 (MQ=255)
                                                                                        ttAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCAccc                                                      >  1:2940517/1‑148 (MQ=255)
                                                                                          aaGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGt                                                    <  2:3167387/148‑1 (MQ=255)
                                                                                          aaGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCGGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGt                                                    >  1:325286/1‑148 (MQ=255)
                                                                                            gCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTAc                                                  >  1:2343208/1‑148 (MQ=255)
                                                                                            gCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTAc                                                  <  1:2292932/148‑1 (MQ=255)
                                                                                              aCCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCCCCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTAccc                                                >  2:3257301/1‑148 (MQ=255)
                                                                                               ccGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCg                                               >  2:687525/1‑148 (MQ=255)
                                                                                                  tCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACt                                            >  1:1684150/1‑148 (MQ=255)
                                                                                                     cATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTAcc                                         >  2:141113/1‑148 (MQ=255)
                                                                                                        aaCGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAg                                      >  1:186046/1‑148 (MQ=255)
                                                                                                        aaCGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAg                                      >  2:1750833/1‑148 (MQ=255)
                                                                                                        aaCGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAg                                      >  1:2282853/1‑148 (MQ=255)
                                                                                                         aCGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGc                                     <  1:818611/148‑1 (MQ=255)
                                                                                                         aCGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGc                                     <  2:2131359/148‑1 (MQ=255)
                                                                                                         aCGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGc                                     >  1:357097/1‑148 (MQ=255)
                                                                                                          cGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCa                                    <  1:778132/148‑1 (MQ=255)
                                                                                                          cGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCa                                    <  1:3195537/148‑1 (MQ=255)
                                                                                                          cGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCa                                    <  1:3195472/148‑1 (MQ=255)
                                                                                                            ggTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATg                                  >  1:1174217/1‑148 (MQ=255)
                                                                                                             gTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATgg                                 <  1:2048259/148‑1 (MQ=255)
                                                                                                                   aCCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGa                           >  2:727062/1‑148 (MQ=255)
                                                                                                                    ccGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGac                          >  2:924333/1‑148 (MQ=255)
                                                                                                                    ccGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGac                          >  2:1554782/1‑148 (MQ=255)
                                                                                                                    ccGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGac                          >  1:816261/1‑148 (MQ=255)
                                                                                                                    ccGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGac                          >  1:815862/1‑148 (MQ=255)
                                                                                                                     cGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGaca                         >  2:131051/1‑148 (MQ=255)
                                                                                                                     cGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGaca                         >  2:411136/1‑148 (MQ=255)
                                                                                                                       ttttGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACAcc                       >  2:2938306/1‑148 (MQ=255)
                                                                                                                       ttttGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACAcc                       >  2:2938305/1‑148 (MQ=255)
                                                                                                                        tttGCTGTACCCACTGCAAGTAACCAACATCGTTCTACAATCGGTTAACCAGCCAGCAAAACCACCATGACTGATATCCAGTTTATCAAAGCCAAACCAGCTCCAGGTTTTCACCCCTACCCGACTACCGACCATGGATGCGACACCg                      >  2:3228070/1‑148 (MQ=255)
                                                                                                                        tttGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCg                      <  1:251738/148‑1 (MQ=255)
                                                                                                                           gCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGAtt                   <  2:2822366/148‑1 (MQ=255)
                                                                                                                           gCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGAtt                   <  2:537552/148‑1 (MQ=255)
                                                                                                                            cTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGAttt                  <  1:2842992/148‑1 (MQ=255)
                                                                                                                                tGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGc              <  2:2956199/148‑1 (MQ=255)
                                                                                                                                 gCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGcc             <  2:1166587/148‑1 (MQ=255)
                                                                                                                                    aCTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCg          >  2:2252635/1‑148 (MQ=255)
                                                                                                                                     cTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGt         <  1:1633102/148‑1 (MQ=255)
                                                                                                                                     cTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGt         >  2:1415349/1‑148 (MQ=255)
                                                                                                                                      tGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTc        <  2:2196750/148‑1 (MQ=255)
                                                                                                                                       gCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTcc       >  1:2996613/1‑148 (MQ=255)
                                                                                                                                           gTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGca   <  2:1303894/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
TGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCAC  >  CP009273/2322464‑2322751

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: