Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,768,689 C→A 100% S215* (TCA→TAA)  ydiB → quinate/shikimate 5‑dehydrogenase, NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,768,6890CA98.2% 880.2 / 5.6 277S215* (TCA→TAA) ydiBquinate/shikimate 5‑dehydrogenase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base C (4/1);  new base A (115/157);  total (119/158)
Fisher's exact test for biased strand distribution p-value = 1.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-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 282 total aligned reads displayed.

CTTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGGG  >  CP009273/1768543‑1768834
                                                                                                                                                  |                                                                                                                                                 
cTTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGaataat                                                                                                                                                 >  1:114830/1‑149 (MQ=255)
cTTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGaataat                                                                                                                                                 >  1:113308/1‑149 (MQ=255)
 ttCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAAtt                                                                                                                                                >  1:3228674/1‑149 (MQ=255)
    gcgcAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGt                                                                                                                                             >  1:2497402/1‑149 (MQ=255)
      gcAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTa                                                                                                                                           >  2:3118672/1‑149 (MQ=255)
       cAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTaa                                                                                                                                          >  2:3673508/1‑149 (MQ=255)
       cAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTaa                                                                                                                                          <  2:3423584/149‑1 (MQ=255)
        aGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAAt                                                                                                                                         <  1:3627944/149‑1 (MQ=255)
          cGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCTATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGa                                                                                                                                       <  2:2912778/149‑1 (MQ=255)
          cGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGa                                                                                                                                       <  1:2770999/149‑1 (MQ=255)
           gggTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGat                                                                                                                                      >  2:2616982/1‑149 (MQ=255)
            ggTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGata                                                                                                                                     <  2:1180967/149‑1 (MQ=255)
              ttAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATc                                                                                                                                   <  2:3434697/149‑1 (MQ=255)
                aaTGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAg                                                                                                                                 >  1:3705666/1‑149 (MQ=255)
                 aTGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGt                                                                                                                                >  2:1697543/1‑149 (MQ=255)
                  tGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTc                                                                                                                               >  2:2519154/1‑149 (MQ=255)
                  tGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTc                                                                                                                               >  1:3722195/1‑149 (MQ=255)
                    aaaaCACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTg                                                                                                                             >  2:1518091/1‑149 (MQ=255)
                     aaaCACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGt                                                                                                                            <  2:1038253/149‑1 (MQ=255)
                      aaCTCCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCATTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGACTAATTGGTTAATGATATCAGTCTGtt                                                                                                                           <  1:3463415/149‑1 (MQ=255)
                      aaCACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGtt                                                                                                                           >  1:1777661/1‑149 (MQ=255)
                          ccGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACAt                                                                                                                       <  2:2651534/149‑1 (MQ=255)
                            gATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATcc                                                                                                                     >  1:4091918/1‑149 (MQ=255)
                                gtgtCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGa                                                                                                                 >  1:1332894/1‑149 (MQ=255)
                                gtgtCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGa                                                                                                                 >  1:980339/1‑149 (MQ=255)
                                   tcgtcACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACtt                                                                                                              >  1:3855221/1‑149 (MQ=255)
                                     gtcACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCt                                                                                                            <  1:2534618/149‑1 (MQ=255)
                                     gtcACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCt                                                                                                            <  1:2974494/149‑1 (MQ=255)
                                       cACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTgg                                                                                                          <  1:3011428/149‑1 (MQ=255)
                                        aCGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGt                                                                                                         <  1:3066489/149‑1 (MQ=255)
                                         cGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTc                                                                                                        >  1:1914090/1‑149 (MQ=255)
                                           gTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCAc                                                                                                      >  2:3959218/1‑149 (MQ=255)
                                            tCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACt                                                                                                     >  1:903078/1‑149 (MQ=255)
                                            tCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACt                                                                                                     <  1:793181/149‑1 (MQ=255)
                                            tCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACt                                                                                                     <  2:2005740/149‑1 (MQ=255)
                                               ccGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGaa                                                                                                  >  1:3764061/1‑149 (MQ=255)
                                               ccGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGaa                                                                                                  <  2:2358098/149‑1 (MQ=255)
                                                cGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAAt                                                                                                 <  2:1827712/149‑1 (MQ=255)
                                                  aTCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGc                                                                                               >  2:3208371/1‑149 (MQ=255)
                                                  aTCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGc                                                                                               >  1:1562383/1‑149 (MQ=255)
                                                       gCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGta                                                                                          <  2:1258268/149‑1 (MQ=255)
                                                          gATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATaa                                                                                       >  2:2765867/1‑149 (MQ=255)
                                                          gATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATaa                                                                                       >  1:2378384/1‑149 (MQ=255)
                                                          gATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATaa                                                                                       <  2:153982/149‑1 (MQ=255)
                                                               gCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGc                                                                                  <  2:2840016/149‑1 (MQ=255)
                                                                 aaGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCat                                                                                >  1:2309428/1‑149 (MQ=255)
                                                                  aGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCata                                                                               >  2:84071/1‑149 (MQ=255)
                                                                  aGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCata                                                                               <  2:3183049/149‑1 (MQ=255)
                                                                   gCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCatat                                                                              <  2:805765/149‑1 (MQ=255)
                                                                    ccTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATg                                                                             <  1:1841230/149‑1 (MQ=255)
                                                                      tttGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGAc                                                                           <  1:63610/149‑1 (MQ=255)
                                                                      tttGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGAc                                                                           <  1:182813/149‑1 (MQ=255)
                                                                         gCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGaa                                                                        <  1:1306281/149‑1 (MQ=255)
                                                                            gAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGtt                                                                     >  2:1341113/1‑149 (MQ=255)
                                                                            gAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGtt                                                                     <  2:2919593/149‑1 (MQ=255)
                                                                            gAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGtt                                                                     <  2:2446783/149‑1 (MQ=255)
                                                                             aaGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTa                                                                    <  1:765379/149‑1 (MQ=255)
                                                                              aGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTAt                                                                   >  2:3252626/1‑149 (MQ=255)
                                                                              aGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTAt                                                                   >  2:2225134/1‑149 (MQ=255)
                                                                              aGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTAt                                                                   <  1:326832/149‑1 (MQ=255)
                                                                                  cTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTgca                                                               <  2:3250758/149‑1 (MQ=255)
                                                                                  cTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTgca                                                               <  1:4048191/149‑1 (MQ=255)
                                                                                  cTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTgca                                                               <  2:1198463/149‑1 (MQ=255)
                                                                                       ttCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAgg                                                          >  1:1968369/1‑149 (MQ=255)
                                                                                        tCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGgc                                                         <  1:2168644/149‑1 (MQ=255)
                                                                                            ccgACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGcaa                                                     <  2:2782073/149‑1 (MQ=255)
                                                                                             cgACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGcaac                                                    <  2:482799/149‑1 (MQ=255)
                                                                                                cATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAg                                                 <  2:3792154/149‑1 (MQ=255)
                                                                                                   tttAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTg                                              <  1:1499138/149‑1 (MQ=255)
                                                                                                      aaCCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGtt                                           <  1:1400251/149‑1 (MQ=255)
                                                                                                        ccAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGc                                         >  2:2043424/1‑149 (MQ=255)
                                                                                                             ggCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAAc                                    >  1:2027763/1‑149 (MQ=255)
                                                                                                               cacaAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGa                                  <  2:2393231/149‑1 (MQ=255)
                                                                                                               cacaAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGa                                  <  1:1388942/149‑1 (MQ=255)
                                                                                                                acaacaGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGAt                                 <  1:3311359/149‑1 (MQ=255)
                                                                                                                 caAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGAtt                                >  2:756128/1‑149 (MQ=255)
                                                                                                                  aaaaGTGTGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTg                               <  2:2676689/149‑1 (MQ=255)
                                                                                                                   aaaGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGa                              <  1:1697543/149‑1 (MQ=255)
                                                                                                                   aaaGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGa                              <  1:1184572/149‑1 (MQ=255)
                                                                                                                   aaaGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGa                              <  2:3586659/149‑1 (MQ=255)
                                                                                                                          gTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATAc                       <  2:1587992/149‑1 (MQ=255)
                                                                                                                           tATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACg                      >  2:169379/1‑149 (MQ=255)
                                                                                                                           tATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACg                      <  2:242770/149‑1 (MQ=255)
                                                                                                                           tATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACg                      <  2:3492410/149‑1 (MQ=255)
                                                                                                                             tGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACTCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGc                    >  1:3923962/1‑149 (MQ=255)
                                                                                                                             tGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGc                    <  2:2943206/149‑1 (MQ=255)
                                                                                                                              gAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCa                   >  2:3925021/1‑149 (MQ=255)
                                                                                                                                  ccccTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCAtgtt               <  2:557057/149‑1 (MQ=255)
                                                                                                                                  ccccTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCAtgtt               <  2:209925/149‑1 (MQ=255)
                                                                                                                                   cccTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCAtgttg              >  2:415371/1‑149 (MQ=255)
                                                                                                                                    ccTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCAtgttgt             >  1:527479/1‑149 (MQ=255)
                                                                                                                                     cTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCAtgttgtt            <  2:2736799/149‑1 (MQ=255)
                                                                                                                                       tGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCNGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCAtgttgttgt          <  1:2084330/149‑1 (MQ=255)
                                                                                                                                         agaATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTgg        >  2:3341038/1‑149 (MQ=255)
                                                                                                                                          gaatgaatAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGc       >  2:2471976/1‑149 (MQ=255)
                                                                                                                                           aatgaatAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCa      >  1:834898/1‑149 (MQ=255)
                                                                                                                                            atgaatAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCaa     >  1:1673029/1‑149 (MQ=255)
                                                                                                                                            atgaatAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCaa     <  2:2572200/149‑1 (MQ=255)
                                                                                                                                            atgaatAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCAGCATATGACGAAGTTATTGAAGAAGGCGCAACAAGCTGTTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCaa     >  1:1429131/1‑149 (MQ=255)
                                                                                                                                              gaatAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAgg   <  2:665990/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CTTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGGG  >  CP009273/1768543‑1768834

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: