Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,297,485 G→C 100% G44G (GGC→GGG alsK ← D‑allose kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,297,4850GC97.7% 859.1 / 8.2 213G44G (GGC→GGGalsKD‑allose kinase
Reads supporting (aligned to +/- strand):  ref base G (1/4);  new base C (112/96);  total (113/100)
Fisher's exact test for biased strand distribution p-value = 1.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-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 213 total aligned reads displayed.

CTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATGAGGCTGTACCTC  >  CP009273/4297343‑4297629
                                                                                                                                              |                                                                                                                                                
 tGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGa                                                                                                                                            >  2:1443854/1‑148 (MQ=255)
 tGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGa                                                                                                                                            >  1:1502703/1‑148 (MQ=255)
    tAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCa                                                                                                                                         <  1:1447292/148‑1 (MQ=255)
      aCGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAAt                                                                                                                                       >  2:375495/1‑148 (MQ=255)
       cGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATg                                                                                                                                      >  2:3474385/1‑148 (MQ=255)
         gCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGAc                                                                                                                                    >  1:2664578/1‑148 (MQ=255)
         gCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGAc                                                                                                                                    >  1:2667237/1‑148 (MQ=255)
            ggTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTc                                                                                                                                 <  1:2523087/148‑1 (MQ=255)
             gTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCt                                                                                                                                >  1:3178036/1‑148 (MQ=255)
                   cGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTc                                                                                                                          <  2:2067924/148‑1 (MQ=255)
                     tAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCg                                                                                                                        >  2:3131843/1‑148 (MQ=255)
                     tAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCg                                                                                                                        <  1:1044607/148‑1 (MQ=255)
                       gAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCt                                                                                                                      >  2:1638722/1‑148 (MQ=255)
                         aataatGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCttt                                                                                                                    >  2:3243236/1‑148 (MQ=255)
                         aataatGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCttt                                                                                                                    <  1:3260975/148‑1 (MQ=255)
                            aatGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCtttttt                                                                                                                 >  2:3499218/1‑148 (MQ=255)
                            aatGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCtttttt                                                                                                                 <  2:3362206/148‑1 (MQ=255)
                                gTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGc                                                                                                             >  2:2243197/1‑148 (MQ=255)
                                 tGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCa                                                                                                            <  1:3401839/148‑1 (MQ=255)
                                 tGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCa                                                                                                            <  2:1737341/148‑1 (MQ=255)
                                    gTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAgtg                                                                                                         <  1:1747178/148‑1 (MQ=255)
                                     tttATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAgtgt                                                                                                        <  2:731962/148‑1 (MQ=255)
                                     tttATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAgtgt                                                                                                        >  2:2460586/1‑148 (MQ=255)
                                      ttATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTa                                                                                                       >  2:2239463/1‑148 (MQ=255)
                                      ttATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTa                                                                                                       >  2:528351/1‑148 (MQ=255)
                                       tATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAg                                                                                                      >  1:3372967/1‑148 (MQ=255)
                                         tCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCg                                                                                                    >  1:2383528/1‑148 (MQ=255)
                                         tCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCg                                                                                                    >  1:320639/1‑148 (MQ=255)
                                          cTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGt                                                                                                   <  2:2712578/148‑1 (MQ=255)
                                           tttACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGtt                                                                                                  >  1:525732/1‑148 (MQ=255)
                                           tttACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGtt                                                                                                  <  1:1607973/148‑1 (MQ=255)
                                              aCTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCa                                                                                               >  2:1334430/1‑148 (MQ=255)
                                               cTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCAc                                                                                              <  2:670146/148‑1 (MQ=255)
                                                  aCCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCtt                                                                                           >  1:2440650/1‑148 (MQ=255)
                                                  aCCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCtt                                                                                           >  1:2265773/1‑148 (MQ=255)
                                                   ccAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTc                                                                                          <  1:2609630/148‑1 (MQ=255)
                                                   ccAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTc                                                                                          <  2:2508706/148‑1 (MQ=255)
                                                   ccAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTc                                                                                          <  1:585773/148‑1 (MQ=255)
                                                    cAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTct                                                                                         >  2:2777069/1‑148 (MQ=255)
                                                    cAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTct                                                                                         >  2:1024862/1‑148 (MQ=255)
                                                     aGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTctg                                                                                        <  1:2903617/148‑1 (MQ=255)
                                                     aGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTctg                                                                                        <  1:2087935/148‑1 (MQ=255)
                                                          cGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTc                                                                                   <  2:1246790/148‑1 (MQ=255)
                                                           ggAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTcc                                                                                  >  1:2601942/1‑148 (MQ=255)
                                                            gAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCg                                                                                 >  1:1664920/1‑148 (MQ=255)
                                                             aaaTCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGc                                                                                >  2:883072/1‑148 (MQ=255)
                                                             aaaTCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGc                                                                                <  2:190080/148‑1 (MQ=255)
                                                             aaaTCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGc                                                                                >  2:2261280/1‑148 (MQ=255)
                                                               aTCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAg                                                                              <  1:1052596/148‑1 (MQ=255)
                                                                  ccATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACa                                                                           <  2:3178036/148‑1 (MQ=255)
                                                                     tCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACaaaa                                                                        >  2:1681368/1‑148 (MQ=255)
                                                                     tCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACaaaa                                                                        <  2:1111827/148‑1 (MQ=255)
                                                                      cACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAg                                                                       <  2:761001/148‑1 (MQ=255)
                                                                      cACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAg                                                                       <  1:1090207/148‑1 (MQ=255)
                                                                         cAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCgg                                                                    >  2:3050315/1‑148 (MQ=255)
                                                                         cAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCgg                                                                    <  2:876997/148‑1 (MQ=255)
                                                                           gACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGat                                                                  <  2:1739488/148‑1 (MQ=255)
                                                                            aCCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGata                                                                 <  1:12339/148‑1 (MQ=255)
                                                                            aCCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGata                                                                 >  1:369861/1‑148 (MQ=255)
                                                                                tGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCg                                                             <  2:257137/148‑1 (MQ=255)
                                                                                   cAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTg                                                          <  2:793450/148‑1 (MQ=255)
                                                                                   cAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTg                                                          >  1:3476747/1‑148 (MQ=255)
                                                                                          cgTTAACGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCat                                                   <  2:1204360/148‑1 (MQ=255)
                                                                                            ttAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCagat                                                 >  2:2380829/1‑148 (MQ=255)
                                                                                             tAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATc                                                <  1:3512758/148‑1 (MQ=255)
                                                                                               aaGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCa                                              >  2:1727132/1‑148 (MQ=255)
                                                                                                aGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCCTATCCAc                                             >  2:301947/1‑148 (MQ=255)
                                                                                                aGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCAc                                             >  1:102345/1‑148 (MQ=255)
                                                                                                  cgcCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGc                                           <  1:1270970/148‑1 (MQ=255)
                                                                                                  cgcCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGc                                           <  1:676090/148‑1 (MQ=255)
                                                                                                     cTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCg                                        <  1:1484650/148‑1 (MQ=255)
                                                                                                     cTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCg                                        <  2:2917813/148‑1 (MQ=255)
                                                                                                     cTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCg                                        <  1:1797826/148‑1 (MQ=255)
                                                                                                       gAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCt                                      <  1:2002370/148‑1 (MQ=255)
                                                                                                        aGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTa                                     >  2:2693186/1‑148 (MQ=255)
                                                                                                        aGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTa                                     <  2:288329/148‑1 (MQ=255)
                                                                                                        aGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTa                                     >  1:2222943/1‑148 (MQ=255)
                                                                                                          ttGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTacg                                   >  2:850360/1‑148 (MQ=255)
                                                                                                             cTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTacgacg                                <  1:3261936/148‑1 (MQ=255)
                                                                                                             cTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTacgacg                                <  1:1574577/148‑1 (MQ=255)
                                                                                                               cGTCAATAATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGtt                              >  2:3124857/1‑148 (MQ=255)
                                                                                                                  cAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCCGACAAAAGCGGATATGCGTTGCCCCCCTATCCACGCCCGCTACGACGTTATg                           >  2:3126908/1‑148 (MQ=255)
                                                                                                                  cAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATg                           <  1:884329/148‑1 (MQ=255)
                                                                                                                   aaTCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGc                          >  2:487772/1‑148 (MQ=255)
                                                                                                                      cATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGt                       >  1:2301053/1‑148 (MQ=255)
                                                                                                                        tttCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGttt                     <  2:2665830/148‑1 (MQ=255)
                                                                                                                         ttCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGtttt                    >  2:926613/1‑148 (MQ=255)
                                                                                                                           cGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTg                  <  2:2948299/148‑1 (MQ=255)
                                                                                                                           cGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTg                  >  1:2798391/1‑148 (MQ=255)
                                                                                                                             ccGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCa                <  1:142705/148‑1 (MQ=255)
                                                                                                                               gATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCCTATCCACGCCCGCTACGACGGTATGCTGTTTTTGCATg              >  1:604304/1‑148 (MQ=255)
                                                                                                                               gATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGGTGTTTTTGCATg              >  1:604445/1‑148 (MQ=255)
                                                                                                                               gATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATg              <  1:3499239/148‑1 (MQ=255)
                                                                                                                                  aCCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATGAgg           <  1:1563380/148‑1 (MQ=255)
                                                                                                                                  aCCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCAGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATGAgg           <  1:528351/148‑1 (MQ=255)
                                                                                                                                   cccGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATGAGGc          <  1:2923703/148‑1 (MQ=255)
                                                                                                                                       acacCAGGCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATGAGGCTGTa      <  2:2383528/148‑1 (MQ=255)
                                                                                                                                       acacCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATGAGGCTGTa      >  1:930839/1‑148 (MQ=255)
                                                                                                                                         acCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATGAGGCTGTAcc    <  1:3350714/148‑1 (MQ=255)
                                                                                                                                           cAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATGAGGCTGTACCTc  >  2:3459480/1‑148 (MQ=255)
                                                                                                                                              |                                                                                                                                                
CTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATGAGGCTGTACCTC  >  CP009273/4297343‑4297629

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: