Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,893,446 T→A 100% K63* (AAA→TAA)  ygcW ← putative dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,893,4460TA98.7% 1253.3 / 2.3 387K63* (AAA→TAA) ygcWputative dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (1/4);  new base A (168/214);  total (169/218)
Fisher's exact test for biased strand distribution p-value = 3.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-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 392 total aligned reads displayed.

CCAGCACCTTATTCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGA  >  CP009273/2893299‑2893592
                                                                                                                                                   |                                                                                                                                                  
ccAGCACCTTATTCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGTAACGCTCACTGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGttatt                                                                                                                                                  <  1:847371/150‑1 (MQ=255)
ccAGCACCTTATTCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGttatt                                                                                                                                                  <  2:491824/150‑1 (MQ=255)
ccAGCACCTTATTCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGttatt                                                                                                                                                  >  1:2617457/1‑150 (MQ=255)
     aCCTTATTCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATc                                                                                                                                             >  1:4624600/1‑150 (MQ=255)
     aCCTTATTCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATc                                                                                                                                             >  2:467266/1‑150 (MQ=255)
       cTTATTCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCAGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCAt                                                                                                                                           >  2:2737388/1‑150 (MQ=255)
        ttATTCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCAtt                                                                                                                                          <  2:4830429/150‑1 (MQ=255)
         tATTCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCAttt                                                                                                                                         >  2:867546/1‑150 (MQ=255)
          aTTCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTc                                                                                                                                        <  2:1738368/150‑1 (MQ=255)
          aTTCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTc                                                                                                                                        <  2:3444846/150‑1 (MQ=255)
             cAGCTTACAAATACCGTCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCtt                                                                                                                                     <  2:3329610/150‑1 (MQ=255)
             cAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCtt                                                                                                                                     <  2:4357524/150‑1 (MQ=255)
               gCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTg                                                                                                                                   <  2:4186619/150‑1 (MQ=255)
                 ttACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGtt                                                                                                                                 <  1:4750933/150‑1 (MQ=255)
                 ttACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGtt                                                                                                                                 <  1:3675011/150‑1 (MQ=255)
                 ttACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGtt                                                                                                                                 <  2:1654413/150‑1 (MQ=255)
                  tACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGttt                                                                                                                                <  1:4825586/150‑1 (MQ=255)
                    cAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCACGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCg                                                                                                                              >  1:5130726/1‑150 (MQ=255)
                    cAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCg                                                                                                                              <  1:2220906/150‑1 (MQ=255)
                    cAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCg                                                                                                                              >  1:1574865/1‑150 (MQ=255)
                    cAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCg                                                                                                                              >  2:1142408/1‑150 (MQ=255)
                    cAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCAGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCg                                                                                                                              >  1:567708/1‑150 (MQ=255)
                           cGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTAt                                                                                                                       <  2:3524324/150‑1 (MQ=255)
                            ggCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATc                                                                                                                      >  1:1000795/1‑150 (MQ=255)
                             gCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCt                                                                                                                     >  1:3384945/1‑150 (MQ=255)
                             gCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCt                                                                                                                     <  1:3413828/150‑1 (MQ=255)
                               aTTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCttt                                                                                                                   <  1:4960745/150‑1 (MQ=255)
                               aTTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCttt                                                                                                                   <  2:2813900/150‑1 (MQ=255)
                                  gTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGAc                                                                                                                >  2:1582731/1‑150 (MQ=255)
                                  gTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGAc                                                                                                                >  1:4805830/1‑150 (MQ=255)
                                   ttAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACg                                                                                                               <  2:5000826/150‑1 (MQ=255)
                                   ttAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACg                                                                                                               <  1:3201810/150‑1 (MQ=255)
                                    tAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGa                                                                                                              <  1:3612782/150‑1 (MQ=255)
                                     aaCCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGaa                                                                                                             <  2:1934062/150‑1 (MQ=255)
                                       ccAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAAc                                                                                                           <  1:2402331/150‑1 (MQ=255)
                                        cAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACt                                                                                                          <  2:1300515/150‑1 (MQ=255)
                                          gAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAg                                                                                                        >  2:4452096/1‑150 (MQ=255)
                                            atatCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGa                                                                                                      >  1:4316684/1‑150 (MQ=255)
                                              atCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAAt                                                                                                    >  1:144352/1‑150 (MQ=255)
                                              atCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATAACCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACTAAAATAGTAAt                                                                                                    >  1:1987367/1‑150 (MQ=255)
                                                  aCTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAg                                                                                                <  1:235758/150‑1 (MQ=255)
                                                   cTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGa                                                                                               >  1:689170/1‑150 (MQ=255)
                                                      tACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGatat                                                                                            <  2:2433723/150‑1 (MQ=255)
                                                       aCCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATAtt                                                                                           <  2:3122644/150‑1 (MQ=255)
                                                        ccGACACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATAttt                                                                                          <  1:3843634/150‑1 (MQ=255)
                                                         cGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTg                                                                                         <  2:1438409/150‑1 (MQ=255)
                                                          gAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTgc                                                                                        <  2:3796184/150‑1 (MQ=255)
                                                              cGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCa                                                                                    <  1:381841/150‑1 (MQ=255)
                                                               gCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAg                                                                                   >  1:1647478/1‑150 (MQ=255)
                                                                 tCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCt                                                                                 <  2:2214586/150‑1 (MQ=255)
                                                                 tCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCAGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCt                                                                                 <  1:3105296/150‑1 (MQ=255)
                                                                   acaGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCttt                                                                               <  1:736382/150‑1 (MQ=255)
                                                                   acaGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCttt                                                                               <  2:646539/150‑1 (MQ=255)
                                                                     agcaagcaGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTgg                                                                             >  1:2801665/1‑150 (MQ=255)
                                                                           cagcGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACg                                                                       >  2:5125086/1‑150 (MQ=255)
                                                                               gATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAGCTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCAt                                                                   >  2:3962724/1‑150 (MQ=255)
                                                                               gATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCAt                                                                   >  2:374427/1‑150 (MQ=255)
                                                                               gATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCAt                                                                   <  2:2617457/150‑1 (MQ=255)
                                                                               gATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCAt                                                                   >  1:2950925/1‑150 (MQ=255)
                                                                               gATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCAt                                                                   <  2:2344614/150‑1 (MQ=255)
                                                                               gATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCAt                                                                   >  1:1108990/1‑150 (MQ=255)
                                                                                aTAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATg                                                                  >  1:897379/1‑150 (MQ=255)
                                                                                aTAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCGTTGGCCAACGCCATg                                                                  >  1:1692943/1‑150 (MQ=255)
                                                                                  aaTCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGc                                                                >  2:524964/1‑150 (MQ=255)
                                                                                   aTCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCa                                                               >  2:4149735/1‑150 (MQ=255)
                                                                                   aTCTTCTGCGGCGCGCCTTCAGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCa                                                               >  1:1869559/1‑150 (MQ=255)
                                                                                      ttctGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAAt                                                            <  2:4038588/150‑1 (MQ=255)
                                                                                       tctGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATg                                                           <  1:3920695/150‑1 (MQ=255)
                                                                                       tctGCGGCGCGCCTTCAGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATg                                                           <  2:2189592/150‑1 (MQ=255)
                                                                                        ctGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGc                                                          <  2:4624600/150‑1 (MQ=255)
                                                                                          gcggcgCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCAGTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCt                                                        <  2:2981300/150‑1 (MQ=255)
                                                                                              cgcgcCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGcc                                                    <  2:3605644/150‑1 (MQ=255)
                                                                                                         cGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGAAGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTAt                                         <  1:366965/150‑1 (MQ=255)
                                                                                                          cgTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGTCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTAta                                        >  1:4428167/2‑149 (MQ=255)
                                                                                                           gTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTc                                       >  1:2478507/1‑150 (MQ=255)
                                                                                                            tGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTcc                                      <  1:2115899/150‑1 (MQ=255)
                                                                                                              aTACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCa                                    <  1:5021490/150‑1 (MQ=255)
                                                                                                                aCCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCAcc                                  <  1:412690/150‑1 (MQ=255)
                                                                                                                aCCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCAcc                                  >  1:92670/1‑150 (MQ=255)
                                                                                                                 cccACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCg                                 <  1:755449/150‑1 (MQ=255)
                                                                                                                  ccACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCgg                                >  1:1753707/1‑150 (MQ=255)
                                                                                                                      cTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCACATGCCTGGCCTAAACCGCTCTTCCCACCGGTAAc                            >  1:3029095/1‑150 (MQ=255)
                                                                                                                        gCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACaa                          <  1:1246334/150‑1 (MQ=255)
                                                                                                                          aTGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAAtt                        <  2:3433226/150‑1 (MQ=255)
                                                                                                                           tGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTg                       >  1:3532763/1‑150 (MQ=255)
                                                                                                                           tGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTg                       >  2:496893/1‑150 (MQ=255)
                                                                                                                             aaGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTCTTCCCACCGGTAACAATTGCg                     >  2:718904/1‑150 (MQ=255)
                                                                                                                               gTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGt                   <  2:3387260/150‑1 (MQ=255)
                                                                                                                                   aCCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTa               <  1:538584/150‑1 (MQ=255)
                                                                                                                                    ccTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTAc              <  2:1831787/150‑1 (MQ=255)
                                                                                                                                    ccTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCAGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTACCACCGGTAACAATTGCGGTTTTAc              <  2:4100886/150‑1 (MQ=255)
                                                                                                                                     cTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTAcc             >  2:149846/1‑150 (MQ=255)
                                                                                                                                        aaCACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCttt          >  2:4130024/1‑150 (MQ=255)
                                                                                                                                          cacCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCa        <  1:487316/150‑1 (MQ=255)
                                                                                                                                           acCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAg       <  1:4481892/150‑1 (MQ=255)
                                                                                                                                           acCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAg       <  2:935449/150‑1 (MQ=255)
                                                                                                                                            cccTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAgg      <  1:4919575/150‑1 (MQ=255)
                                                                                                                                             ccTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAggg     <  1:3477348/150‑1 (MQ=255)
                                                                                                                                             ccTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAggg     <  2:1331940/150‑1 (MQ=255)
                                                                                                                                              cTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGga    <  2:897387/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CCAGCACCTTATTCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGA  >  CP009273/2893299‑2893592

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: