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.5% 1017.7 / 3.3 322S215* (TCA→TAA) ydiBquinate/shikimate 5‑dehydrogenase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base C (0/5);  new base A (124/190);  total (124/198)
Fisher's exact test for biased strand distribution p-value = 1.61e-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 329 total aligned reads displayed.

CCTTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGGGG  >  CP009273/1768542‑1768835
                                                                                                                                                   |                                                                                                                                                  
 cTTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAAtt                                                                                                                                                 >  1:457538/1‑150 (MQ=255)
    cgcgCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGt                                                                                                                                              <  2:4016416/150‑1 (MQ=255)
    cgcgCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGt                                                                                                                                              <  2:4168424/150‑1 (MQ=255)
       gcAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTaa                                                                                                                                           <  1:3477266/150‑1 (MQ=255)
       gcAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTaa                                                                                                                                           <  2:4584342/150‑1 (MQ=255)
        cAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAAt                                                                                                                                          <  2:2612779/150‑1 (MQ=255)
         aGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATg                                                                                                                                         >  2:628901/1‑150 (MQ=255)
             ggTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGatat                                                                                                                                     <  1:2218699/150‑1 (MQ=255)
              gTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATc                                                                                                                                    <  1:3665938/150‑1 (MQ=255)
              gTTAAAGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATc                                                                                                                                    <  2:3530316/150‑1 (MQ=255)
                tAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAg                                                                                                                                  >  2:3849126/1‑150 (MQ=255)
                tAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAg                                                                                                                                  <  2:1853819/150‑1 (MQ=255)
                   tGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCt                                                                                                                               <  2:3159936/150‑1 (MQ=255)
                      aaaCACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGtt                                                                                                                            <  1:5204531/150‑1 (MQ=255)
                            cGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGTCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATcc                                                                                                                      >  1:1573187/1‑150 (MQ=255)
                            cGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATcc                                                                                                                      <  2:5043908/150‑1 (MQ=255)
                            cGAGTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATcc                                                                                                                      <  1:1847463/150‑1 (MQ=255)
                               ttGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCggg                                                                                                                   <  1:2952927/150‑1 (MQ=255)
                                 gtgtCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGAc                                                                                                                 <  2:2380172/150‑1 (MQ=255)
                                   gtcgtcACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACtt                                                                                                               <  1:3513324/150‑1 (MQ=255)
                                     cgtcACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCt                                                                                                             <  1:379388/150‑1 (MQ=255)
                                     cgtcACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCt                                                                                                             <  2:2148518/150‑1 (MQ=255)
                                      gtcACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTg                                                                                                            <  1:3627197/150‑1 (MQ=255)
                                      gtcACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTg                                                                                                            >  1:1249496/1‑150 (MQ=255)
                                        cACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGt                                                                                                          <  1:4018986/150‑1 (MQ=255)
                                            gTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACt                                                                                                      <  2:3407041/150‑1 (MQ=255)
                                            gTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACt                                                                                                      <  2:1443076/150‑1 (MQ=255)
                                             tCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTGCATCCGGGACTTCTGGTCACTg                                                                                                     <  2:7772/150‑1 (MQ=255)
                                               aCCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGaa                                                                                                   >  2:3025993/1‑150 (MQ=255)
                                               aCCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACACCTAGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGaa                                                                                                   >  1:151383/1‑150 (MQ=255)
                                                 cGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATg                                                                                                 <  2:1085229/150‑1 (MQ=255)
                                                  gATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGc                                                                                                >  2:1272105/1‑150 (MQ=255)
                                                  gATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGc                                                                                                >  2:4551232/1‑150 (MQ=255)
                                                  gATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTGGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGc                                                                                                >  2:309480/1‑150 (MQ=255)
                                                        gCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGtat                                                                                          <  1:3597824/150‑1 (MQ=255)
                                                         ccGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGtata                                                                                         <  2:1798921/150‑1 (MQ=255)
                                                           gATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAAc                                                                                       >  1:1806448/1‑150 (MQ=255)
                                                             tCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAAccc                                                                                     <  2:1636767/150‑1 (MQ=255)
                                                                 cAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTCACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCat                                                                                 <  2:4460617/150‑1 (MQ=255)
                                                                    gCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATg                                                                              <  2:2707292/150‑1 (MQ=255)
                                                                     ccTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGa                                                                             <  2:313645/150‑1 (MQ=255)
                                                                      cTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGAc                                                                            <  1:2305189/150‑1 (MQ=255)
                                                                      cTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGAc                                                                            <  2:1366276/150‑1 (MQ=255)
                                                                         tGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGaa                                                                         <  1:4133707/150‑1 (MQ=255)
                                                                          gCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAg                                                                        <  1:2787840/150‑1 (MQ=255)
                                                                            tGAAGCCCTGGCTTCCGCCGACATTTTAACCAATTGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATACTTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGtt                                                                      <  2:1824970/150‑1 (MQ=255)
                                                                            tGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTTGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGAGTGTATAACCCGCATATGACGAAGtt                                                                      >  1:1892944/1‑150 (MQ=255)
                                                                             gAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTa                                                                     >  1:2148710/1‑150 (MQ=255)
                                                                             gAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTa                                                                     <  2:4249178/150‑1 (MQ=255)
                                                                               aGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTAtt                                                                   >  2:1194345/1‑150 (MQ=255)
                                                                                 cccTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTgc                                                                 <  2:3619627/150‑1 (MQ=255)
                                                                                  ccTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTgca                                                                <  1:957568/150‑1 (MQ=255)
                                                                                    tGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTgcagc                                                              >  1:5278694/1‑150 (MQ=255)
                                                                                    tGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTgcagc                                                              >  1:3915788/1‑150 (MQ=255)
                                                                                    tGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTgcagc                                                              <  1:2441836/150‑1 (MQ=255)
                                                                                      gCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGcagcag                                                            <  2:3888030/150‑1 (MQ=255)
                                                                                             ccgACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGcaac                                                     <  2:1035367/150‑1 (MQ=255)
                                                                                                aCATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAg                                                  <  1:1115826/150‑1 (MQ=255)
                                                                                                aCATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAg                                                  <  1:201932/150‑1 (MQ=255)
                                                                                                  aTTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCt                                                >  2:1008189/1‑150 (MQ=255)
                                                                                                  aTTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCt                                                >  1:2634814/1‑150 (MQ=255)
                                                                                                  aTTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCt                                                <  1:1182479/150‑1 (MQ=255)
                                                                                                   ttttAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTg                                               <  1:1410959/150‑1 (MQ=255)
                                                                                                    tttAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTgg                                              >  2:2768352/1‑150 (MQ=255)
                                                                                                       aaCCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTg                                           >  1:2324856/1‑150 (MQ=255)
                                                                                                         ccAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCa                                         >  1:2540755/1‑150 (MQ=255)
                                                                                                         ccAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCa                                         <  1:3488604/150‑1 (MQ=255)
                                                                                                           aaTGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCaaa                                       >  1:3063560/1‑150 (MQ=255)
                                                                                                             tGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAAc                                     <  1:3807708/150‑1 (MQ=255)
                                                                                                              ggCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACg                                    <  1:4331956/150‑1 (MQ=255)
                                                                                                                cacaAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGAt                                  >  2:5202309/1‑150 (MQ=255)
                                                                                                                 acaAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGAtt                                 <  2:4861632/150‑1 (MQ=255)
                                                                                                                  caAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTg                                <  2:3296308/150‑1 (MQ=255)
                                                                                                                  caAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTg                                >  1:3377255/1‑150 (MQ=255)
                                                                                                                   aaaaGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGa                               <  1:75467/150‑1 (MQ=255)
                                                                                                                   aaaaGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGa                               >  2:5317732/1‑150 (MQ=255)
                                                                                                                    aaaGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGAt                              <  2:22293/150‑1 (MQ=255)
                                                                                                                     aaGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATg                             <  2:4298451/150‑1 (MQ=255)
                                                                                                                        tGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGAt                          <  1:3697405/150‑1 (MQ=255)
                                                                                                                             aTGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGc                     <  1:4909884/150‑1 (MQ=255)
                                                                                                                               gAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCAt                   >  2:1383952/1‑150 (MQ=255)
                                                                                                                                aaaCCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCAtg                  <  1:169510/150‑1 (MQ=255)
                                                                                                                                aaaCCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCAtg                  >  2:659794/1‑150 (MQ=255)
                                                                                                                                  aCCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCAtgtt                <  2:3777679/150‑1 (MQ=255)
                                                                                                                                   ccccTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCAtgttg               <  2:2419368/150‑1 (MQ=255)
                                                                                                                                   ccccTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCAtgttg               <  1:1058307/150‑1 (MQ=255)
                                                                                                                                   ccccTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCAtgttg               <  2:1936287/150‑1 (MQ=255)
                                                                                                                                    cccTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCAtgttgt              <  2:264689/150‑1 (MQ=255)
                                                                                                                                     ccTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCAtgttgtt             >  2:3675133/1‑150 (MQ=255)
                                                                                                                                      cTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCAtgttgttg            <  2:4742070/150‑1 (MQ=255)
                                                                                                                                      cTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCAtgttgttg            >  1:3197708/1‑150 (MQ=255)
                                                                                                                                       ttGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCAtgttgttgt           >  1:4602210/1‑150 (MQ=255)
                                                                                                                                        tGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTtgtg          >  1:4726804/1‑150 (MQ=255)
                                                                                                                                        tGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTtgtg          <  2:2884695/150‑1 (MQ=255)
                                                                                                                                          agaATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGc        >  2:841703/1‑150 (MQ=255)
                                                                                                                                              tgaatAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAgg    >  2:4746343/1‑150 (MQ=255)
                                                                                                                                                aataatTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAgggg  <  1:2715930/150‑1 (MQ=255)
                                                                                                                                                aataatTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAgggg  >  1:906679/1‑150 (MQ=255)
                                                                                                                                                aataatTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAgggg  >  2:2205365/1‑150 (MQ=255)
                                                                                                                                                aataatTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAgggg  <  2:2120385/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CCTTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGGGG  >  CP009273/1768542‑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: