Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 288,249 T→A 100% T137T (ACA→ACT)  yagK ← CP4‑6 prophage; conserved protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273288,2490TA99.0% 968.2 / ‑0.6 291T137T (ACA→ACT) yagKCP4‑6 prophage; conserved protein
Reads supporting (aligned to +/- strand):  ref base T (3/0);  new base A (126/162);  total (129/162)
Fisher's exact test for biased strand distribution p-value = 8.60e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01
Rejected as polymorphism: E-value score below prediction cutoff.
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 292 total aligned reads displayed.

GGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAACGGT  >  CP009273/288104‑288394
                                                                                                                                                 |                                                                                                                                                 
  tAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGat                                                                                                                                              <  1:822203/149‑1 (MQ=255)
    gTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGatca                                                                                                                                            >  1:1492201/1‑149 (MQ=255)
        aGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATccc                                                                                                                                        <  2:3998531/149‑1 (MQ=255)
        aGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATccc                                                                                                                                        <  2:1698365/149‑1 (MQ=255)
         gTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCct                                                                                                                                       >  1:1002305/1‑149 (MQ=255)
            gggTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTa                                                                                                                                    >  1:2924432/1‑149 (MQ=255)
             ggTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTaa                                                                                                                                   <  2:3572171/149‑1 (MQ=255)
             ggTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTaa                                                                                                                                   >  1:3225134/1‑149 (MQ=255)
              gTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAg                                                                                                                                  >  1:2968416/1‑149 (MQ=255)
              gTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAg                                                                                                                                  <  2:672238/149‑1 (MQ=255)
                    aGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATc                                                                                                                            >  2:3500920/1‑149 (MQ=255)
                     gcgcCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCg                                                                                                                           >  1:331058/1‑149 (MQ=255)
                        ccTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTcc                                                                                                                        >  2:2063490/1‑149 (MQ=255)
                         cTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTATCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCt                                                                                                                       <  2:370699/149‑1 (MQ=255)
                         cTGGTAGTTCTGCTGGAAGTCAGTACTATGGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCt                                                                                                                       <  1:3248058/149‑1 (MQ=255)
                            gTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGAt                                                                                                                    <  2:2467940/149‑1 (MQ=255)
                            gTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGAt                                                                                                                    <  2:400489/149‑1 (MQ=255)
                                ttCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTa                                                                                                                <  2:2330710/149‑1 (MQ=255)
                                ttCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTa                                                                                                                >  2:1348797/1‑149 (MQ=255)
                                 tCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTaa                                                                                                               <  1:273086/149‑1 (MQ=255)
                                 tCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTaa                                                                                                               >  2:455203/1‑149 (MQ=255)
                                  ctgctgGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAAt                                                                                                              <  1:2391558/149‑1 (MQ=255)
                                     ctgGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATccc                                                                                                           <  2:1210054/149‑1 (MQ=255)
                                        gAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGa                                                                                                        >  2:103618/1‑149 (MQ=255)
                                          aGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAgg                                                                                                      <  2:1575135/149‑1 (MQ=255)
                                          aGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAgg                                                                                                      <  1:3507291/149‑1 (MQ=255)
                                           gTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAggt                                                                                                     >  1:2249340/1‑149 (MQ=255)
                                               gTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGgta                                                                                                 <  2:1025022/149‑1 (MQ=255)
                                                    aTTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAg                                                                                            <  1:213403/149‑1 (MQ=255)
                                                    aTTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAg                                                                                            <  1:212343/149‑1 (MQ=255)
                                                      tGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCg                                                                                          <  2:3225134/149‑1 (MQ=255)
                                                        gtgtCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTc                                                                                        <  2:1164299/149‑1 (MQ=255)
                                                         tgtCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCt                                                                                       <  1:1614976/149‑1 (MQ=255)
                                                         tgtCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCt                                                                                       >  1:570192/1‑149 (MQ=255)
                                                           tCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCttt                                                                                     <  2:3047075/149‑1 (MQ=255)
                                                            cTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTttg                                                                                    <  2:54506/149‑1 (MQ=255)
                                                            cTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTttg                                                                                    <  1:62089/149‑1 (MQ=255)
                                                             tAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTttgt                                                                                   >  1:1508786/1‑149 (MQ=255)
                                                             tAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTttgt                                                                                   <  2:623346/149‑1 (MQ=255)
                                                                cacaTACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGa                                                                                <  1:1713722/149‑1 (MQ=255)
                                                                 acaTACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGaa                                                                               >  1:1910995/1‑149 (MQ=255)
                                                                  caTACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGaac                                                                              <  1:732701/149‑1 (MQ=255)
                                                                  caTACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGaac                                                                              <  1:2685669/149‑1 (MQ=255)
                                                                    tACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGaacaa                                                                            >  1:1672/1‑149 (MQ=255)
                                                                    tACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGaacaa                                                                            <  2:2823204/149‑1 (MQ=255)
                                                                    tACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCAGCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGaacaa                                                                            <  1:1744851/149‑1 (MQ=255)
                                                                       ttACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAg                                                                         >  1:961457/1‑149 (MQ=255)
                                                                         aCAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGc                                                                       <  2:601361/149‑1 (MQ=255)
                                                                                  aGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGAt                                                              <  1:3429129/149‑1 (MQ=255)
                                                                                  aGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGAt                                                              >  1:85333/1‑149 (MQ=255)
                                                                                  aGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGAt                                                              <  2:56435/149‑1 (MQ=255)
                                                                                   ggAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATa                                                             <  2:961457/149‑1 (MQ=255)
                                                                                   ggAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATa                                                             <  2:1609954/149‑1 (MQ=255)
                                                                                      aaaTGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTg                                                          >  2:1176883/1‑149 (MQ=255)
                                                                                       aaTGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTgg                                                         >  1:3758074/1‑149 (MQ=255)
                                                                                         tGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAATGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCa                                                       <  2:2862665/149‑1 (MQ=255)
                                                                                          ggACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCAt                                                      <  1:352624/149‑1 (MQ=255)
                                                                                          ggACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCAt                                                      <  2:3341721/149‑1 (MQ=255)
                                                                                           gaccagaccaGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCAtt                                                     <  1:1840381/149‑1 (MQ=255)
                                                                                              cagaccagGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTAc                                                  <  1:1348797/149‑1 (MQ=255)
                                                                                               agaccagGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTAcc                                                 >  1:4015873/1‑149 (MQ=255)
                                                                                                gaccagGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCt                                                <  1:3531886/149‑1 (MQ=255)
                                                                                                    agGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGAtt                                            >  1:988711/1‑149 (MQ=255)
                                                                                                     gggTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTc                                           <  2:182459/149‑1 (MQ=255)
                                                                                                     gggTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTc                                           >  1:1281636/1‑149 (MQ=255)
                                                                                                     gggTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTc                                           <  2:2226412/149‑1 (MQ=255)
                                                                                                      ggTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCa                                          >  2:295792/1‑149 (MQ=255)
                                                                                                        tGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCACATGTGATAGTGGCATTTACCTGATTCAGa                                        >  1:2762597/1‑149 (MQ=255)
                                                                                                         gATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGaa                                       <  1:1079526/149‑1 (MQ=255)
                                                                                                          aTCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAAt                                      >  1:3436014/1‑149 (MQ=255)
                                                                                                          aTCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAAt                                      >  2:3795380/1‑149 (MQ=255)
                                                                                                            cGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATAc                                    <  1:3500920/149‑1 (MQ=255)
                                                                                                              tctctctcCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATActc                                  >  1:1983790/1‑149 (MQ=255)
                                                                                                              tctctctcCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATActc                                  >  1:4069823/1‑149 (MQ=255)
                                                                                                               ctctctcCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATActct                                 <  1:3731148/149‑1 (MQ=255)
                                                                                                                 ctctcCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCttt                               <  2:425007/149‑1 (MQ=255)
                                                                                                                 ctctcCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCttt                               <  2:3271095/149‑1 (MQ=255)
                                                                                                                  tctcCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTc                              >  1:1314761/1‑149 (MQ=255)
                                                                                                                    tcCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGc                            >  2:1408038/1‑149 (MQ=255)
                                                                                                                    tcCAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGc                            >  2:3161249/1‑149 (MQ=255)
                                                                                                                     ccAGTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGcc                           >  1:364367/1‑149 (MQ=255)
                                                                                                                        gTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCat                        >  2:1926976/1‑149 (MQ=255)
                                                                                                                        gTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCat                        >  1:498189/1‑149 (MQ=255)
                                                                                                                        gTCTCAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCat                        <  2:1500871/149‑1 (MQ=255)
                                                                                                                         tctcAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCata                       <  1:3084106/149‑1 (MQ=255)
                                                                                                                         tctcAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCata                       <  2:2632888/149‑1 (MQ=255)
                                                                                                                            cAGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATataa                    <  1:3839839/149‑1 (MQ=255)
                                                                                                                             aGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATataat                   <  2:2699209/149‑1 (MQ=255)
                                                                                                                             aGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATataat                   <  2:1527907/149‑1 (MQ=255)
                                                                                                                             aGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATataat                   <  1:3141021/149‑1 (MQ=255)
                                                                                                                             aGGGCGCTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATataat                   >  1:3551147/1‑149 (MQ=255)
                                                                                                                                 cgcTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATaataaat               <  2:3737828/149‑1 (MQ=255)
                                                                                                                                 cgcTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATaataaat               <  2:351917/149‑1 (MQ=255)
                                                                                                                                  gcTGTACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATaataaata              >  1:3114126/1‑149 (MQ=255)
                                                                                                                                      tACCAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAgg          <  2:2180943/149‑1 (MQ=255)
                                                                                                                                        ccAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGc        <  2:2906642/149‑1 (MQ=255)
                                                                                                                                        ccAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGc        <  1:863502/149‑1 (MQ=255)
                                                                                                                                        ccAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGc        <  2:847719/149‑1 (MQ=255)
                                                                                                                                         cAGGCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCa       <  1:3692888/149‑1 (MQ=255)
                                                                                                                                           ggCTCCAGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAAc     <  1:254859/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATGTGATAGTGGCATTTACCTGATTCAGAATACTCTTTCGCCCATATAATAAATAAAGGGCAACGGT  >  CP009273/288104‑288394

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