Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,885,611 G→T 93.7% R628S (CGT→AGT)  yoaA ← putative ATP‑dependent helicase, DinG family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,885,6110GT93.7% 472.7 / 21.8 167R628S (CGT→AGT) yoaAputative ATP‑dependent helicase, DinG family
Reads supporting (aligned to +/- strand):  ref base G (6/5);  new base T (98/58);  total (104/63)
Fisher's exact test for biased strand distribution p-value = 7.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 169 total aligned reads displayed.

TCGTTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACGAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATCAGTCGCCCTACCCCCTG  >  CP009273/1885466‑1885758
                                                                                                                                                 |                                                                                                                                                   
 cgtTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAAc                                                                                                                                                  <  1:2969790/148‑1 (MQ=255)
 cgtTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAAc                                                                                                                                                  <  2:2974794/148‑1 (MQ=255)
 cgtTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGGGCTTAGATCAATAAAATGGCGCGCATCATACCATACTCCGTAACAAATTACCTGCAGGATGGTATCGCAAGGAAACTAAc                                                                                                                                                  >  2:2537734/1‑148 (MQ=255)
          gggCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGc                                                                                                                                         >  1:1686385/1‑148 (MQ=255)
           ggCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCa                                                                                                                                        <  1:2949761/148‑1 (MQ=255)
           ggCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCa                                                                                                                                        <  1:811117/148‑1 (MQ=255)
               aCAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGt                                                                                                                                    <  1:2726888/148‑1 (MQ=255)
                cAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTc                                                                                                                                   >  1:1623153/1‑148 (MQ=255)
                   agCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACg                                                                                                                                >  2:2364004/1‑148 (MQ=255)
                     cAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCCAGGAAACTAACCGCACGGGCAATGTCACggg                                                                                                                              >  2:2332574/1‑148 (MQ=255)
                     cAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACggg                                                                                                                              <  2:3130597/148‑1 (MQ=255)
                        gCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGc                                                                                                                           >  1:1604207/1‑148 (MQ=255)
                            ctGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTgg                                                                                                                       <  2:2199729/148‑1 (MQ=255)
                                cgcgGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGcgcg                                                                                                                   >  1:3034073/1‑148 (MQ=255)
                                cgcgGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACGCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACGAACCGCACGGGCAATGTCACGGGTGCGTGGcgcg                                                                                                                   <  2:420942/148‑1 (MQ=255)
                                     tATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCgcgggcgg                                                                                                              >  1:1013669/1‑148 (MQ=255)
                                      aTCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGggcggc                                                                                                             >  2:1539290/1‑148 (MQ=255)
                                       tCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCa                                                                                                            >  1:352647/1‑148 (MQ=255)
                                       tCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACGAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCa                                                                                                            >  2:2684574/1‑148 (MQ=255)
                                              ccAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGc                                                                                                     >  1:1955462/1‑148 (MQ=255)
                                              ccAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGc                                                                                                     <  2:546093/148‑1 (MQ=255)
                                              ccAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGc                                                                                                     >  1:1003483/1‑148 (MQ=255)
                                              ccAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGc                                                                                                     >  1:420884/1‑148 (MQ=255)
                                              ccAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGc                                                                                                     >  1:420735/1‑148 (MQ=255)
                                              ccAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGc                                                                                                     <  2:2619884/148‑1 (MQ=255)
                                                   aTTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGaa                                                                                                >  2:2756854/1‑148 (MQ=255)
                                                     tCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAAc                                                                                              >  2:345462/1‑148 (MQ=255)
                                                        cATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGAGCATCATACCATACTCCGTAACAAATTAACTCGAGTATGGTATCTCAAGGAAACTAACAGCACTGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTc                                                                                           >  1:200101/1‑148 (MQ=255)
                                                                  cGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAg                                                                                 <  1:1431098/148‑1 (MQ=255)
                                                                   gTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAgg                                                                                <  2:1380497/148‑1 (MQ=255)
                                                                    tGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGa                                                                               >  1:2699084/1‑148 (MQ=255)
                                                                    tGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACGAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGa                                                                               >  1:868544/1‑148 (MQ=255)
                                                                     gCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGAc                                                                              >  2:136305/1‑148 (MQ=255)
                                                                     gCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGAc                                                                              >  2:928202/1‑148 (MQ=255)
                                                                     gCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGAc                                                                              >  1:125108/1‑148 (MQ=255)
                                                                       ttAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGc                                                                            >  1:2005687/1‑148 (MQ=255)
                                                                       ttAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGc                                                                            >  2:1145634/1‑148 (MQ=255)
                                                                        tAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCa                                                                           >  2:2661548/1‑148 (MQ=255)
                                                                        tAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCa                                                                           >  2:2661544/1‑148 (MQ=255)
                                                                        tAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCa                                                                           >  2:255085/1‑148 (MQ=255)
                                                                          gATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATc                                                                         <  1:3282430/148‑1 (MQ=255)
                                                                          gATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATc                                                                         >  2:2898881/1‑148 (MQ=255)
                                                                          gATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATc                                                                         >  2:2766952/1‑148 (MQ=255)
                                                                          gATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATc                                                                         <  1:2504489/148‑1 (MQ=255)
                                                                           aTCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCa                                                                        >  1:2369064/1‑148 (MQ=255)
                                                                           aTCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCa                                                                        >  1:2369077/1‑148 (MQ=255)
                                                                            tCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATCTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCAc                                                                       >  2:2349352/1‑148 (MQ=255)
                                                                               aTAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAg                                                                    >  2:2180102/1‑148 (MQ=255)
                                                                                tAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGc                                                                   >  2:2378848/1‑148 (MQ=255)
                                                                                tAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGc                                                                   >  2:1678416/1‑148 (MQ=255)
                                                                                  aaaaGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCACGGAAACGAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCg                                                                 >  2:291106/1‑148 (MQ=255)
                                                                                  aaaaGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCg                                                                 <  2:2696120/148‑1 (MQ=255)
                                                                                    aaGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGAt                                                               >  1:3147488/1‑148 (MQ=255)
                                                                                    aaGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGAt                                                               >  2:221731/1‑148 (MQ=255)
                                                                                    aaGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGAt                                                               >  1:3147551/1‑148 (MQ=255)
                                                                                    aaGGCGCGCATCAGACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGAt                                                               <  1:2609808/148‑1 (MQ=255)
                                                                                     gggCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGAtt                                                              >  1:2978713/2‑148 (MQ=255)
                                                                                     aGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGAtt                                                              <  1:1210802/148‑1 (MQ=255)
                                                                                      ggCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTg                                                             >  1:3231654/1‑148 (MQ=255)
                                                                                       gcgcgcATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGt                                                            >  1:290825/1‑148 (MQ=255)
                                                                                            catcatACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACaa                                                       >  2:882881/1‑148 (MQ=255)
                                                                                            catcatACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACaa                                                       >  2:1878296/1‑148 (MQ=255)
                                                                                              tcatACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAAt                                                     <  1:2747273/148‑1 (MQ=255)
                                                                                              tcatACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAAt                                                     <  2:2941108/148‑1 (MQ=255)
                                                                                               cataccatacTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATc                                                    >  1:2103063/1‑148 (MQ=255)
                                                                                                 taccatacTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCAc                                                  <  2:104686/148‑1 (MQ=255)
                                                                                                 taccatacTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCAc                                                  <  2:1152699/148‑1 (MQ=255)
                                                                                                 taccatacTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCAc                                                  <  1:263721/148‑1 (MQ=255)
                                                                                                 taccatacTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCAc                                                  <  1:1840879/148‑1 (MQ=255)
                                                                                                   ccatacTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCa                                                >  1:1552537/1‑148 (MQ=255)
                                                                                                    catacTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAg                                               <  1:1702089/148‑1 (MQ=255)
                                                                                                     atacTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGc                                              >  1:480555/1‑148 (MQ=255)
                                                                                                               aCAAATTACCTGGAGGATGGTATCGCAAGGAAACGAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGAt                                    <  2:1527385/148‑1 (MQ=255)
                                                                                                                cAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGAtc                                   >  2:405708/1‑148 (MQ=255)
                                                                                                                cAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGAtc                                   >  2:406482/1‑148 (MQ=255)
                                                                                                                  aaTTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGAtcgt                                 >  2:383729/1‑148 (MQ=255)
                                                                                                                  aaTTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGAtcgt                                 <  1:2190118/148‑1 (MQ=255)
                                                                                                                  aaTTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGAtcgt                                 <  1:3296910/148‑1 (MQ=255)
                                                                                                                   aTTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGAtcgtc                                >  2:2797597/1‑148 (MQ=255)
                                                                                                                       ccTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCg                            <  1:669458/148‑1 (MQ=255)
                                                                                                                         tGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTc                          <  1:602801/148‑1 (MQ=255)
                                                                                                                           gaggaTGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTcgc                        >  1:426806/1‑148 (MQ=255)
                                                                                                                              gaTGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGaa                     >  1:915560/1‑148 (MQ=255)
                                                                                                                              gaTGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGaa                     <  2:1956807/148‑1 (MQ=255)
                                                                                                                              gaTGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGaa                     <  2:733082/148‑1 (MQ=255)
                                                                                                                              gaTGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGAGCGTCGGCGTCGCGaa                     >  2:1840797/1‑148 (MQ=255)
                                                                                                                                tGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATc                   <  1:1977079/148‑1 (MQ=255)
                                                                                                                                tGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATc                   <  1:2276098/148‑1 (MQ=255)
                                                                                                                                tGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATc                   >  1:1476038/1‑148 (MQ=255)
                                                                                                                                tGGTATCGCAAGGAAACGAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATc                   <  2:832147/148‑1 (MQ=255)
                                                                                                                                 ggTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATCa                  >  2:2434239/1‑148 (MQ=255)
                                                                                                                                 ggTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATCa                  >  2:2387017/1‑148 (MQ=255)
                                                                                                                                    aTCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATCAGTc               >  2:686495/1‑148 (MQ=255)
                                                                                                                                       gctaGGAAACTAACCGCACGGTCAATGTCACTGGTGCTTGGCGCGGGCTGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCTCAAATCACCAGCACGCCACGATCGTGGTCGTCGCGAATCTTTCGcc            <  2:2577375/145‑1 (MQ=255)
                                                                                                                                       gCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATCAGTCGcc            <  2:1323988/148‑1 (MQ=255)
                                                                                                                                       gCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATCAGTCGcc            <  2:1323975/148‑1 (MQ=255)
                                                                                                                                           ggAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATCAGTCGCCCTAc        >  2:2796512/1‑148 (MQ=255)
                                                                                                                                            gAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATCAGTCGCCCTAcc       >  2:2119994/1‑148 (MQ=255)
                                                                                                                                               aCGAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATCAGTCGCCCTAccccc    <  2:2834214/148‑1 (MQ=255)
                                                                                                                                                 gAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATCAGTCGCCCTACCCCCTg  <  2:608319/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TCGTTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACGAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATCAGTCGCCCTACCCCCTG  >  CP009273/1885466‑1885758

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