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,6890CA95.2% 635.0 / 21.0 218S215* (TCA→TAA) ydiBquinate/shikimate 5‑dehydrogenase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base C (6/5);  new base A (123/84);  total (129/89)
Fisher's exact test for biased strand distribution p-value = 7.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 221 total aligned reads displayed.

TTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGGGG  >  CP009273/1768544‑1768835
                                                                                                                                                 |                                                                                                                                                  
ttCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGaataat                                                                                                                                                  >  2:2118413/1‑148 (MQ=255)
    cgcAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGACTGAATAATTGGt                                                                                                                                              >  2:3128600/1‑148 (MQ=255)
    cgcAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGt                                                                                                                                              <  1:2759784/148‑1 (MQ=255)
      cAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTa                                                                                                                                            <  1:887906/148‑1 (MQ=255)
        gCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAAt                                                                                                                                          >  1:36507/1‑148 (MQ=255)
          gggTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGa                                                                                                                                        <  1:1851171/148‑1 (MQ=255)
           ggTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGat                                                                                                                                       <  1:2137438/148‑1 (MQ=255)
           ggTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGat                                                                                                                                       <  1:2065475/148‑1 (MQ=255)
            gTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGata                                                                                                                                      >  2:673843/1‑148 (MQ=255)
            gTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAAGAATTGGTTAATGata                                                                                                                                      >  2:3037634/1‑148 (MQ=255)
             ttAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGatat                                                                                                                                     <  2:1538879/148‑1 (MQ=255)
             ttAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGatat                                                                                                                                     <  1:1010561/148‑1 (MQ=255)
             ttAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGatat                                                                                                                                     <  2:2275059/148‑1 (MQ=255)
              tAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATc                                                                                                                                    >  2:1547680/1‑148 (MQ=255)
              tAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATc                                                                                                                                    >  1:2414634/1‑148 (MQ=255)
               aaTGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTTCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCa                                                                                                                                   <  2:3314879/148‑1 (MQ=255)
               aaTGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCa                                                                                                                                   <  2:2212077/148‑1 (MQ=255)
               aaTGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCa                                                                                                                                   <  1:2466086/148‑1 (MQ=255)
                aTGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAg                                                                                                                                  <  2:2485855/148‑1 (MQ=255)
                 tGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGt                                                                                                                                 >  1:746500/1‑148 (MQ=255)
                 tGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGt                                                                                                                                 <  2:1687085/148‑1 (MQ=255)
                 tGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGt                                                                                                                                 >  2:1490486/1‑148 (MQ=255)
                   aaaaCACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCt                                                                                                                               <  2:1349040/148‑1 (MQ=255)
                      acacCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGtt                                                                                                                            >  1:1440993/1‑148 (MQ=255)
                       cacCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTa                                                                                                                           >  2:1130014/1‑148 (MQ=255)
                       cacCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTa                                                                                                                           >  2:1130027/1‑148 (MQ=255)
                       cacCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTa                                                                                                                           >  2:121202/1‑148 (MQ=255)
                         ccGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACa                                                                                                                         <  1:686257/148‑1 (MQ=255)
                          cGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACAt                                                                                                                        <  2:1902103/148‑1 (MQ=255)
                              tgtgTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCgg                                                                                                                    >  2:278392/1‑148 (MQ=255)
                                tgtCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGa                                                                                                                  <  2:1095283/148‑1 (MQ=255)
                                 gtcgtcACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGAc                                                                                                                 <  1:2999936/148‑1 (MQ=255)
                                   cgtcACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACtt                                                                                                               >  2:1821082/1‑148 (MQ=255)
                                    gtcACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTc                                                                                                              <  2:3211115/148‑1 (MQ=255)
                                    gtcACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTc                                                                                                              <  2:365761/148‑1 (MQ=255)
                                     tcACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCt                                                                                                             >  1:3063459/1‑148 (MQ=255)
                                       aCGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTgg                                                                                                           >  1:642977/1‑148 (MQ=255)
                                        cGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGt                                                                                                          <  2:2722479/148‑1 (MQ=255)
                                         ggTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTc                                                                                                         <  2:2564214/148‑1 (MQ=255)
                                           tCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCAc                                                                                                       <  1:373674/148‑1 (MQ=255)
                                             aCCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTg                                                                                                     <  1:3218539/148‑1 (MQ=255)
                                              ccGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGa                                                                                                    <  1:1550726/148‑1 (MQ=255)
                                               cGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGaa                                                                                                   >  1:496288/1‑148 (MQ=255)
                                                gATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAAt                                                                                                  >  1:3527191/1‑148 (MQ=255)
                                                    tcGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCgt                                                                                              <  1:2157645/148‑1 (MQ=255)
                                                        cGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGtat                                                                                          >  1:1145355/1‑148 (MQ=255)
                                                         gATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGtata                                                                                         <  2:1954621/148‑1 (MQ=255)
                                                          aTCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATaa                                                                                        <  2:2056414/148‑1 (MQ=255)
                                                           tCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAAc                                                                                       >  2:921978/1‑148 (MQ=255)
                                                            cAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAAcc                                                                                      >  2:132242/1‑148 (MQ=255)
                                                             aGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAAccc                                                                                     >  2:62650/1‑148 (MQ=255)
                                                                aaGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCa                                                                                  >  1:149357/1‑148 (MQ=255)
                                                                   ccTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCatat                                                                               >  1:2340984/1‑148 (MQ=255)
                                                                    cTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATg                                                                              >  1:1205114/1‑148 (MQ=255)
                                                                     tttGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGa                                                                             <  2:497417/148‑1 (MQ=255)
                                                                           gAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGt                                                                       >  2:3040509/1‑148 (MQ=255)
                                                                           gAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGt                                                                       >  2:332967/1‑148 (MQ=255)
                                                                           gAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGt                                                                       >  2:3522364/1‑148 (MQ=255)
                                                                             aGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGtga                                                                     >  2:3419968/1‑146 (MQ=255)
                                                                               cccTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTAtt                                                                   >  2:3032787/1‑148 (MQ=255)
                                                                               cccTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTAtt                                                                   <  1:2699559/148‑1 (MQ=255)
                                                                                  tGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTgca                                                                <  1:1288673/148‑1 (MQ=255)
                                                                                       tCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAgg                                                           <  2:1598306/148‑1 (MQ=255)
                                                                                        ccgccgACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGACTAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGgc                                                          >  1:2387348/1‑148 (MQ=255)
                                                                                          gccgACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGgcgc                                                        >  1:2878411/1‑148 (MQ=255)
                                                                                            cgACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGcaa                                                      <  2:3010488/148‑1 (MQ=255)
                                                                                                 ttttAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGc                                                 >  2:1053293/1‑148 (MQ=255)
                                                                                                     aaCCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGt                                             >  1:1424560/1‑148 (MQ=255)
                                                                                                     aaCCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGt                                             <  1:3473837/148‑1 (MQ=255)
                                                                                                       ccAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTg                                           >  2:3146415/1‑148 (MQ=255)
                                                                                                       ccAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTg                                           >  1:1118762/1‑148 (MQ=255)
                                                                                                        cAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGc                                          >  1:3088323/1‑148 (MQ=255)
                                                                                                           tGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCaaa                                       <  1:2856767/148‑1 (MQ=255)
                                                                                                            ggCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCaaaa                                      <  1:1999726/148‑1 (MQ=255)
                                                                                                             gCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAAc                                     >  2:2258637/1‑148 (MQ=255)
                                                                                                              cacaAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACg                                    >  1:1472019/1‑148 (MQ=255)
                                                                                                               acaAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGa                                   <  2:2000936/148‑1 (MQ=255)
                                                                                                               acaAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGa                                   >  2:2617569/1‑148 (MQ=255)
                                                                                                                caAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGAt                                  <  1:227153/148‑1 (MQ=255)
                                                                                                                 aaaaGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGAtt                                 <  2:395352/148‑1 (MQ=255)
                                                                                                                   aaGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGa                               >  2:1294657/1‑148 (MQ=255)
                                                                                                                   aaGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGa                               >  1:2632954/1‑148 (MQ=255)
                                                                                                                    aGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGAt                              >  2:2423209/1‑148 (MQ=255)
                                                                                                                       gggTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGa                           <  2:1269146/148‑1 (MQ=255)
                                                                                                                        ggTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGAt                          <  2:183558/148‑1 (MQ=255)
                                                                                                                              aaaCCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCa                    >  2:3545502/1‑148 (MQ=255)
                                                                                                                                aCCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCAtg                  <  1:1945330/148‑1 (MQ=255)
                                                                                                                                  cccTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCAtgtt                <  1:1957717/148‑1 (MQ=255)
                                                                                                                                   ccTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCAtgttg               >  2:888588/1‑148 (MQ=255)
                                                                                                                                      tGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCAtgttgttg            <  2:1554626/148‑1 (MQ=255)
                                                                                                                                       gagaATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCAtgttgttgt           <  1:2441193/148‑1 (MQ=255)
                                                                                                                                        agaATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTtgtg          >  2:1741707/1‑148 (MQ=255)
                                                                                                                                         gaatgaatAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGAGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTgg         >  2:864459/1‑148 (MQ=255)
                                                                                                                                          aatgaatAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGc        >  2:1040781/1‑148 (MQ=255)
                                                                                                                                          aatgaatAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGc        >  1:1678466/1‑148 (MQ=255)
                                                                                                                                           atgaatAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCa       <  1:3128600/148‑1 (MQ=255)
                                                                                                                                           atgaatAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCa       >  2:766539/1‑148 (MQ=255)
                                                                                                                                           atgaatAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCa       <  2:2853306/148‑1 (MQ=255)
                                                                                                                                              aataatTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGCTGGATACGGCATGTTGTTGTGGCAAgg    <  2:436096/148‑1 (MQ=255)
                                                                                                                                                tCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAgggg  >  2:297546/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
TTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGGGG  >  CP009273/1768544‑1768835

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: