Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,209,162 C→A 100% I9I (ATC→ATA yohO → putative membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,209,1620CA98.3% 734.4 / 3.2 229I9I (ATC→ATAyohOputative membrane protein
Reads supporting (aligned to +/- strand):  ref base C (4/0);  new base A (142/83);  total (146/83)
Fisher's exact test for biased strand distribution p-value = 2.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 229 total aligned reads displayed.

CAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATCGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACG  >  CP009273/2209017‑2209306
                                                                                                                                                 |                                                                                                                                                
  gcagcaACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGccc                                                                                                                                              >  2:1969410/1‑148 (MQ=255)
         cTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCcctgttcct                                                                                                                                       >  1:3124644/1‑148 (MQ=255)
          tCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATGGGGGTCATAGCcctgttcctg                                                                                                                                      <  1:2582639/148‑1 (MQ=255)
            aCCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCcctgttcctgtt                                                                                                                                    >  1:702740/1‑148 (MQ=255)
            aCCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCcctgttcctgtt                                                                                                                                    >  2:2730129/1‑148 (MQ=255)
             ccAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGttt                                                                                                                                   >  2:1186891/1‑148 (MQ=255)
              cAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTa                                                                                                                                  >  2:2173462/1‑148 (MQ=255)
               aGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTAt                                                                                                                                 <  2:1819073/148‑1 (MQ=255)
               aGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATGGGGGTCATAGCCCTGTTCCTGTTTAt                                                                                                                                 >  1:3535961/1‑148 (MQ=255)
                     cAAGAACAAGTCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGtt                                                                                                                           >  1:520156/1‑148 (MQ=255)
                      aaGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGAGTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTa                                                                                                                          <  1:2340131/148‑1 (MQ=255)
                       aGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAg                                                                                                                         >  2:3178930/1‑148 (MQ=255)
                            aaGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGa                                                                                                                    >  1:2394351/1‑148 (MQ=255)
                             aGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGaa                                                                                                                   <  2:2204652/148‑1 (MQ=255)
                             aGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGaa                                                                                                                   >  1:1833384/1‑148 (MQ=255)
                              ggCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAAt                                                                                                                  <  1:143248/148‑1 (MQ=255)
                              ggCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAAt                                                                                                                  >  2:670567/1‑148 (MQ=255)
                               gCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAAtt                                                                                                                 >  2:3044558/1‑148 (MQ=255)
                                ccATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATtg                                                                                                                >  1:652551/1‑148 (MQ=255)
                                 cATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATtgg                                                                                                               >  1:375367/1‑148 (MQ=255)
                                  aTGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATtggt                                                                                                              <  1:2323049/148‑1 (MQ=255)
                                    gATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATtggtgg                                                                                                            >  2:410473/1‑148 (MQ=255)
                                      tatTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCg                                                                                                          >  1:2739383/1‑148 (MQ=255)
                                       atTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGt                                                                                                         <  2:1289163/148‑1 (MQ=255)
                                       atTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGt                                                                                                         <  2:1144821/148‑1 (MQ=255)
                                        ttttCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTc                                                                                                        <  2:2951181/148‑1 (MQ=255)
                                             ccACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCt                                                                                                   <  1:1582745/148‑1 (MQ=255)
                                                  cGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAg                                                                                              <  2:1012353/148‑1 (MQ=255)
                                                  cGGCATTTAATGATTAGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATCGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAg                                                                                              >  1:2059855/1‑148 (MQ=255)
                                                     cATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGtt                                                                                           >  1:1721713/1‑148 (MQ=255)
                                                      aTTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTta                                                                                          >  1:1490735/1‑148 (MQ=255)
                                                       tttAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTtat                                                                                         >  2:2629243/1‑148 (MQ=255)
                                                        ttAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTtata                                                                                        >  2:708105/1‑148 (MQ=255)
                                                        ttAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTtata                                                                                        <  1:1272581/148‑1 (MQ=255)
                                                        ttAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTtata                                                                                        >  1:2539629/1‑148 (MQ=255)
                                                          aaTGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATAcc                                                                                      <  1:2980806/148‑1 (MQ=255)
                                                            tGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCtt                                                                                    >  1:2389006/1‑148 (MQ=255)
                                                            tGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCtt                                                                                    >  1:1447508/1‑148 (MQ=255)
                                                             gATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCttt                                                                                   <  2:520156/148‑1 (MQ=255)
                                                             gATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGGTATACCttt                                                                                   >  1:2358684/1‑148 (MQ=255)
                                                                tCGTTAACAAATGCGCTTTACTGTACACTCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTAtt                                                                                <  2:777092/148‑1 (MQ=255)
                                                                tCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTAtt                                                                                <  2:2564953/148‑1 (MQ=255)
                                                                tCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTAtt                                                                                >  1:131544/1‑148 (MQ=255)
                                                                tCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTAtt                                                                                >  1:1675885/1‑148 (MQ=255)
                                                                tCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTAtt                                                                                <  2:2316430/148‑1 (MQ=255)
                                                                 cGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTAttt                                                                               >  1:2681211/1‑148 (MQ=255)
                                                                 cGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTAttt                                                                               >  1:2470869/1‑148 (MQ=255)
                                                                   ttAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTg                                                                             <  2:1125744/148‑1 (MQ=255)
                                                                      aCAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGt                                                                          >  1:267000/1‑148 (MQ=255)
                                                                       cAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTg                                                                         <  2:723473/148‑1 (MQ=255)
                                                                       cAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTg                                                                         >  2:953135/1‑148 (MQ=255)
                                                                       cAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTg                                                                         >  2:1040743/1‑148 (MQ=255)
                                                                       cAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTg                                                                         <  1:1037297/148‑1 (MQ=255)
                                                                        aaaTGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCATGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGc                                                                        >  2:3418924/1‑148 (MQ=255)
                                                                         aaTGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCt                                                                       >  2:3511490/1‑148 (MQ=255)
                                                                          aTGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTg                                                                      >  2:1183278/1‑148 (MQ=255)
                                                                          aTGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTg                                                                      >  2:1975768/1‑148 (MQ=255)
                                                                            gcgcTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGc                                                                    <  1:1248024/148‑1 (MQ=255)
                                                                                 ttACCGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAgc                                                               >  1:903329/1‑148 (MQ=255)
                                                                                   aCTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAgcgc                                                             >  2:1869875/1‑148 (MQ=255)
                                                                                           aTCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCa                                                     >  2:3005937/1‑148 (MQ=255)
                                                                                             ccTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGc                                                   <  2:1746352/148‑1 (MQ=255)
                                                                                             ccTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGc                                                   <  1:2233379/148‑1 (MQ=255)
                                                                                              cTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGcc                                                  >  1:332497/1‑148 (MQ=255)
                                                                                                     tAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCaaa                                           >  2:1574623/1‑148 (MQ=255)
                                                                                                     tAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCaaa                                           >  2:2942934/1‑148 (MQ=255)
                                                                                                     tAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCaaa                                           >  2:2621322/1‑148 (MQ=255)
                                                                                                     tAACATAAGGAGTGCATTAATCGCCTAGCTAAAATTGGGGTCATAGCCCTGTACCTGGTTCTGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTATATTTTGCGTGATGGCTAAGCGCCTGCACCAGCCTTTCaaa                                           >  2:149812/1‑148 (MQ=255)
                                                                                                        cATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAg                                        >  1:3341565/1‑148 (MQ=255)
                                                                                                        cATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAg                                        >  1:107851/1‑148 (MQ=255)
                                                                                                          tAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGc                                      >  1:1435172/1‑148 (MQ=255)
                                                                                                             ggAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGt                                   >  2:2749010/1‑148 (MQ=255)
                                                                                                              gAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTc                                  >  1:338724/1‑148 (MQ=255)
                                                                                                               aGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCt                                 >  1:2751720/1‑148 (MQ=255)
                                                                                                                gTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTg                                <  1:3457178/148‑1 (MQ=255)
                                                                                                                gTGCATTATGCGCATAGCTAAAATGGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTg                                <  2:3535961/148‑1 (MQ=255)
                                                                                                                   cATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCga                             <  2:2756210/148‑1 (MQ=255)
                                                                                                                    aTTATGCGCATAGCTAAAATTGGGGTCATCGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCgat                            >  2:266432/1‑148 (MQ=255)
                                                                                                                    aTTATGCGCATAGCTAAAATTGGGGTCATCGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCgat                            >  2:2658989/1‑148 (MQ=255)
                                                                                                                      tATGCGCATAGCTACAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCgatga                          <  2:41729/148‑1 (MQ=255)
                                                                                                                         gcgcATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATcg                       >  1:1013491/1‑148 (MQ=255)
                                                                                                                         gcgcATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATcg                       >  1:1063630/1‑148 (MQ=255)
                                                                                                                         gcgcATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATcg                       <  1:1768821/148‑1 (MQ=255)
                                                                                                                         gcgcATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATcg                       >  1:1782981/1‑148 (MQ=255)
                                                                                                                          cgcATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATcgc                      >  1:304056/1‑148 (MQ=255)
                                                                                                                          cgcATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATcgc                      <  1:2173452/148‑1 (MQ=255)
                                                                                                                          cgcATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATcgc                      >  1:308239/1‑148 (MQ=255)
                                                                                                                             aTAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATcgccgc                   <  1:1969410/148‑1 (MQ=255)
                                                                                                                             aTAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATcgccgc                   <  2:2539629/148‑1 (MQ=255)
                                                                                                                              tAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCgccgcc                  <  1:2484331/148‑1 (MQ=255)
                                                                                                                                   aaaaTTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGc             <  2:2346444/148‑1 (MQ=255)
                                                                                                                                      aTTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGccac          >  2:3013072/1‑148 (MQ=255)
                                                                                                                                       ttGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGccacc         <  2:2389006/148‑1 (MQ=255)
                                                                                                                                       ttGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGccacc         >  2:352300/1‑148 (MQ=255)
                                                                                                                                       ttGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGccacc         <  2:297275/148‑1 (MQ=255)
                                                                                                                                        tGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGccacca        >  2:1676162/1‑148 (MQ=255)
                                                                                                                                          gggTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGt      <  2:598694/148‑1 (MQ=255)
                                                                                                                                          gggTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGt      <  1:3279867/148‑1 (MQ=255)
                                                                                                                                             tCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAAc   >  2:2756179/1‑148 (MQ=255)
                                                                                                                                              cATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACg  <  1:889954/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATCGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACG  >  CP009273/2209017‑2209306

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: