Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,661,110 C→A 5.5% D123Y (GAT→TAT)  ynfK ← putative dethiobiotin synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,661,1100CA 5.5% 790.2 / 16.7 201D123Y (GAT→TAT) ynfKputative dethiobiotin synthetase
Reads supporting (aligned to +/- strand):  ref base C (91/99);  new base A (4/7);  total (95/106)
Fisher's exact test for biased strand distribution p-value = 5.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 217 total aligned reads displayed.

GTTAGCCACCCAGCCAATGAGCGGCAGCCCGTCGTTGGCGATCGCCTGAGCTGTTAGCAGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGCCAGGCCGTTTGAAATGAGGGTGTAATTGATTGGGCAACTGTGCGCCACGCTACTTTCTTCTTCGCTTAACGCGATAGGAT  >  CP009273/1660961‑1661259
                                                                                                                                                     |                                                                                                                                                     
gTTAGCCACCCAGCCAATGAGCGGCAGCCCGTCGTTGGCGATCGCCTGAGCTGTTAGCAGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAAtc                                                                                                                                                       >  2:1472029/1‑150 (MQ=255)
gTTAGCCACCCAGCCAATGAGCGGCAGCCCGTCGTTGGCGATCGCCTGAGCTGTTAGCAGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAAtc                                                                                                                                                       >  2:1056934/1‑150 (MQ=255)
 ttAGCCACCCAGCCAATGAGCGGCAGCCCGTCGTTGGCGATCGCCAGAGCTGTTAGCAGTTCATTGTTAATGCAACCATCCTGAAAACCGAAAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAAtca                                                                                                                                                      >  2:722408/1‑150 (MQ=255)
  tAGCCACCCAGCCAATGAGCGGCAGCCCGTCGTTGGCGATCGCCTGAGCTGTTAGCAGTGCATGGTTAATGCAACCTTCCTGAATACCGACACCCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAAtcat                                                                                                                                                     <  2:1722173/150‑1 (MQ=255)
   aGCCACCCAGCCAATGAGCGGCAGCCCGTCGTTGGCGATCGCCTGAGCTGTTAGCAGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAAtcatt                                                                                                                                                    <  2:3230508/150‑1 (MQ=255)
   aGCCACCCAGCCAATGAGCGGCAGCCCGTCGTTGGCGATCGCCTGAGCTGTTAGCAGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAAtcatt                                                                                                                                                    <  2:2297916/150‑1 (MQ=255)
    gCCACCCAGCCAATGAGCGGCAGCCCGTCGTTGGCGATCGCCTGAGCTGTTAGCAGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAAtaattc                                                                                                                                                   <  1:1963619/150‑1 (MQ=255)
     ccACCCAGCCAATGAGCGGCAGCCCGTCGTTGGCGATCGCCTGAGCTGTTAGCAGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAAtcattca                                                                                                                                                  <  1:1599946/150‑1 (MQ=255)
       aCCCAGCCAATGAGCGGCAGCCCGTCGTTGGCGATCGCCTGAGCTGTTAGCAGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATtcatc                                                                                                                                                >  1:1037998/1‑150 (MQ=255)
        cccAGCCAATGAGCGGCAGCCCGTCGTTGGCGATCGCCTGAGCTGTTAGCAGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATtcatca                                                                                                                                               <  1:1618168/150‑1 (MQ=255)
           aGCCAATGAGCGGCAGCCCGTCGTTGGCGATCGCCTGAGCTGTTAGCAGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGAc                                                                                                                                            <  2:3795060/150‑1 (MQ=255)
              cAATGAGCGGCAGCCCGTCGTTGGCGATCGCCTGAGCTGTTAGCAGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTgc                                                                                                                                         <  1:4080028/150‑1 (MQ=255)
              cAATGAGCGGCAGCCCGTCGTTGGCGATCGCCTGAGCTGTTAGCAGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTgc                                                                                                                                         <  1:4022942/150‑1 (MQ=255)
               aaTGAGCGGCAGCCCGTCGTTGGCGATCGCCTGAGCTGTTAGCAGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTgcg                                                                                                                                        >  1:1474576/1‑150 (MQ=255)
               aaTGAGCGGCAGCCCGTCGTTGGCGATCGCCTGAGCTGTTAGCAGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTgcg                                                                                                                                        >  2:4077873/1‑150 (MQ=255)
                             cgtcgtTGGCGATCGCCTGAGCTGTTAGCAGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTcc                                                                                                                          <  2:638155/150‑1 (MQ=255)
                               tcgtTGGCGATCGCCTGAGCTGTTAGCAGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCt                                                                                                                        <  1:814470/150‑1 (MQ=255)
                               tcgtTGGCGATCGCCTGAGCTGTTAGCAGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCt                                                                                                                        <  1:2903692/150‑1 (MQ=255)
                                cgtTGGCGATCGCCTGAGCTGTTAGCAGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCtt                                                                                                                       <  1:1215075/150‑1 (MQ=255)
                                cgtTGGCGATCGCCTGAGCTGTTAGCAGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCtt                                                                                                                       >  1:3682945/1‑150 (MQ=255)
                                  ttGGCGATCGCCTGAGCTGTTAGCAGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCt                                                                                                                     <  1:3945049/150‑1 (MQ=255)
                                  ttGGCGATCGCCTGAGCTGTTAGCAGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCt                                                                                                                     <  1:3727853/150‑1 (MQ=255)
                                    ggCGATCGCCTGAGCTGTTAGCAGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTAc                                                                                                                   >  2:3068062/1‑150 (MQ=255)
                                        aTCGCCTGAGCTGTTAGCAGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACg                                                                                                               <  2:3471481/150‑1 (MQ=255)
                                          cGCCTGAGCTGTTAGCAGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGAc                                                                                                             >  1:214365/1‑150 (MQ=255)
                                           gCCTGAGCTGTTAGCAGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACa                                                                                                            >  2:1774584/1‑150 (MQ=255)
                                           gCCTGAGCTGTTAGCAGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACa                                                                                                            >  1:1252680/1‑150 (MQ=255)
                                            ccTGAGCTGTTAGCAGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACAt                                                                                                           >  2:1800373/1‑150 (MQ=255)
                                             cTGAGCTGTTAGCAGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATg                                                                                                          >  2:1694040/1‑150 (MQ=255)
                                             cTGAGCTGTTAGCAGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATg                                                                                                          <  1:2270627/150‑1 (MQ=255)
                                             cTGAGCTGTTAGCAGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATg                                                                                                          >  1:1865963/1‑150 (MQ=255)
                                              tGAGCTGTTAGCAGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGa                                                                                                         >  2:2272067/1‑150 (MQ=255)
                                                aGCTGTTAGCAGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATc                                                                                                       <  2:2454690/150‑1 (MQ=255)
                                                  cTGTTAGCAGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGa                                                                                                     <  2:3935367/150‑1 (MQ=255)
                                                  cTGTTAGCAGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATAATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGa                                                                                                     >  2:2752136/1‑150 (MQ=255)
                                                   tGTTAGCAGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATAATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGAc                                                                                                    <  2:4025262/150‑1 (MQ=255)
                                                      tAGCAGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACttt                                                                                                 >  1:2813057/1‑150 (MQ=255)
                                                       aGCAGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATAATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACtttt                                                                                                <  1:2801921/150‑1 (MQ=255)
                                                          aGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCg                                                                                             >  2:3577116/1‑150 (MQ=255)
                                                          aGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCg                                                                                             <  2:435261/150‑1 (MQ=255)
                                                           gTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCgg                                                                                            <  1:193803/150‑1 (MQ=255)
                                                            tGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGt                                                                                           <  1:1363871/150‑1 (MQ=255)
                                                            tGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGt                                                                                           <  1:1363963/150‑1 (MQ=255)
                                                            tGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGt                                                                                           >  2:3374961/1‑150 (MQ=255)
                                                               aTGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAg                                                                                        >  2:3446627/1‑150 (MQ=255)
                                                                tGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTAACCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAgg                                                                                       >  2:3657548/1‑150 (MQ=255)
                                                                 ggTTAATGCAACCTTCCTGAAGACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGt                                                                                      <  2:3780198/150‑1 (MQ=255)
                                                                     aaTGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGc                                                                                  >  1:1926058/1‑150 (MQ=255)
                                                                         cAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGCCAgg                                                                              <  1:1138255/150‑1 (MQ=255)
                                                                          aaCCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGCCCGGc                                                                             >  2:1560842/1‑150 (MQ=255)
                                                                           aCCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGCCAGGcc                                                                            >  2:1045155/1‑150 (MQ=255)
                                                                           aCCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGCCAGGcc                                                                            >  2:2319834/1‑150 (MQ=255)
                                                                              ttCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGCCTGGCCGtt                                                                         >  2:3739141/1‑150 (MQ=255)
                                                                               tCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGCCAGGCCGttt                                                                        <  2:362379/150‑1 (MQ=255)
                                                                               tCCTGAAGACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGCCAGGCCGttt                                                                        <  2:2566506/150‑1 (MQ=255)
                                                                                 cTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGCCAGGCCGTTTGa                                                                      >  2:2191834/1‑150 (MQ=255)
                                                                                  tGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGCCAGGCCGTTTGaa                                                                     >  2:3750186/1‑150 (MQ=255)
                                                                                   gAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGCCAGGCCGTTTGaaa                                                                    <  1:25341/150‑1 (MQ=255)
                                                                                    aagaCCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGCCAGGCCGTTTGAAAt                                                                   <  1:2624209/147‑1 (MQ=255)
                                                                                          gACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCCCGACATGATCGACTTTTTCGGTCAGGTTTGCCAGGCCGTTTGAAATGAGGgt                                                             <  2:849763/150‑1 (MQ=255)
                                                                                           aCAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGCCAGGCCGTTTGAAATGAGGgtg                                                            <  1:579984/150‑1 (MQ=255)
                                                                                           aCAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGCCAGGCCGTTTGAAATGAGGgtg                                                            <  1:3345822/150‑1 (MQ=255)
                                                                                                cATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTATTTCGGTCAGGTTTGCCAGGCCGTTTGAAATGAGGGTGTAatt                                                       >  2:4024072/1‑150 (MQ=255)
                                                                                                   cAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGCCCGGCCGTTTGAAATGAGGGTGTAattgat                                                    >  2:3799608/1‑150 (MQ=255)
                                                                                                    aaCACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGCCAGGCCGTTTGAAATGAGGGTGTAattgatt                                                   >  1:547175/1‑150 (MQ=255)
                                                                                                       acCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGCCAGGCCGTTTGAAATGAGGGTGTAATTGATTggg                                                >  2:2851339/1‑150 (MQ=255)
                                                                                                       acCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGCCAGGCCGTTTGAAATGAGGGTGTAATTGATTggg                                                >  2:2851322/1‑150 (MQ=255)
                                                                                                        ccGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGCCAGGCCGTTTGAAATGAGGGTGTAATTGATTGGGc                                               <  2:2707484/150‑1 (MQ=255)
                                                                                                            cAGTTGTTCCTGCACTACCCATTCAGAGAGTGGTCGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGCCAGGCCGTTTGAAATGAGGGTGTAATTGATTGGGCAACt                                           <  2:1588334/150‑1 (MQ=255)
                                                                                                            cAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGCCAGGCCGTTTGAAATGAGGGTGTAATTGATTGGGCAACt                                           <  2:562532/150‑1 (MQ=255)
                                                                                                                 gttCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGCCAGGCCGTTTGAACTGAGGGTGTAATTGATTGGGCAACTGTgcg                                      >  2:767924/1‑150 (MQ=255)
                                                                                                                  ttCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGCCAGGCCGTTTGAAATGAGGGTGTAATTGATTGGGCAACTGTgcgc                                     >  2:101105/1‑150 (MQ=255)
                                                                                                                  ttCCTGCACTACCCATTCAGAGAGTGGACGCAAATAATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGCCAGGCCGTTTGAAATGAGGGTGTAATTGATTGGGCAACTGTgcgc                                     <  1:4024019/150‑1 (MQ=255)
                                                                                                                    ccTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGCCAGGCCGTTTGAAATGAGGGTGTAATTGATTGGGCAACTGTGCGCCa                                   >  2:3953939/1‑150 (MQ=255)
                                                                                                                        cACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGCCAGGCCGTTTGAAATGAGGGTGTAATTGATTGGGCAACTGTGCGCCACGCt                               >  2:3488460/1‑150 (MQ=255)
                                                                                                                        cACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGCCAGGCCGTTTGAAATGAGGGTGTAATTGATTGGGCAACTGTGCGCCACGCt                               <  1:507754/150‑1 (MQ=255)
                                                                                                                          cTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGCCAGGCCGTTTGAAATGAGGGTGTAATTGATTGGGCAACTGTGCGCCACGCTAc                             <  2:1153833/150‑1 (MQ=255)
                                                                                                                             cccATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGCCAGGCCGTTTGAAATGAGGGTGTAATTGATTGGGCAACTGTGCGCCACGCTACttt                          >  2:3997763/1‑150 (MQ=255)
                                                                                                                                aTTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGCCAGGCCGTTTGAAATGAGGGTGTAATTGATTGGGCAACTGTGCGCCACGCTACTttctt                       <  1:1056934/150‑1 (MQ=255)
                                                                                                                                aTTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGCCAGGCCGTTTGAAATGAGGGTGTAATTGATTGGGCAACTGTGCGCCACGCTACTttctt                       >  1:2789058/1‑150 (MQ=255)
                                                                                                                                aTTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGCCAGGCCGTTTGAAATGAGGGTGTAATTGATTGGGCAACTGTGCGCCACGCTACTttctt                       >  2:1190445/1‑150 (MQ=255)
                                                                                                                                 ttCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGCCAGGCCGTTTGAAATGAGGGTGTAATTGATTGGGCAACTGTGCGCCACGCTACTttcttc                      >  2:2873041/1‑150 (MQ=255)
                                                                                                                                   cAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGTACAGGTTTGCCAGGCCGTTTGAAATTAGGGTGTAATATATTGGGCAACTGTGCGCCACGCTACTttcttctt                    >  1:160424/1‑150 (MQ=255)
                                                                                                                                    agagagTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGCCAGGCCGTTTGAAATGAGGGTGTAATTGATTGGGCAACTGTGCGCCACGCTACTttcttcttc                   <  2:3030421/150‑1 (MQ=255)
                                                                                                                                    agagagTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGCCAGGCCGTTTGAAATGAGGGTGTAATTGATTGGGCAACTGTGCGCCACGCTACTttcttcttc                   >  1:2945669/1‑150 (MQ=255)
                                                                                                                                     gagagTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGCCAGGCCGTTTGAAATGAGGGTGTAATTGATTGGGCAACTGTGCGCCACGCTACTTTCTTCTTCg                  >  1:12417/1‑150 (MQ=255)
                                                                                                                                       gagTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGCCAGGCCGTTTGAAATGAGGGTGTAATTGATTGGGCAACTGTGCGCCACGCTACTTTCTTCTTCGCt                <  1:455531/150‑1 (MQ=255)
                                                                                                                                        agTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGCCAGGCCGTTTGAAATGAGGGTGTAATTGATTGGGCAACTGTGCGCCACGCTACTTTCTTCTTCGCtt               >  2:3300244/1‑150 (MQ=255)
                                                                                                                                         gTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCCCGACATGATCGACTTTTTCGGTCAGGTTTGCCAGGCCGTTTGAAATGAGGGTGTAATTGATTGGGCAACTGTGCGCCACGCTACTTTCTTCTTCGCTTa              <  1:1892045/150‑1 (MQ=255)
                                                                                                                                          tGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGCCAGGCCGTTTGAAATGAGGGTGTAATTGATTGGGCAACTGTGCGCCACGCTACTTTCTTCTTCGCTTaa             <  1:3472686/150‑1 (MQ=255)
                                                                                                                                          tGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGCCAGGCCGTTTGAAATGAGGGTGTAATTGATTGGGCAACTGTGCGCCACGCTACTTTCTTCTTCGCTTaa             <  1:2242246/150‑1 (MQ=255)
                                                                                                                                           ggACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGCCAGGCCGTTTGAAATGAGGGTGTAATTGATTGGGCAACTGTGCGCCACGCGACTTTCTTCTTCGCGTAAc            >  1:2637364/1‑150 (MQ=255)
                                                                                                                                           ggACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCGACCACGACATGATCGACTTGTTCGGTCAGGTTTGCCAGGCCGTTTGAAATGAGGGTGTAATTGATTGGGCAACTGTGCGCCACGCTACTTTCTTCTTCGCTTAAc            <  2:1455099/150‑1 (MQ=255)
                                                                                                                                            gACGCAACCCAGACAACAGACTGCGCCAGCCGCCCGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGCCAGGCCGTTTGAAATGAGGGTGTAATTGATTGGGCAACTGTGCGCCACGCTACTTTCTTCTTCGCTTAAcg           <  2:2407448/150‑1 (MQ=255)
                                                                                                                                                cAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGCCAGGCCGTTTGAAATGAGGGTGTAATTGATTGGGCAACTGTGCGCCACGCTACTTTCTTCTTCGCTTAACGCGAt       <  2:2834090/150‑1 (MQ=255)
                                                                                                                                                 aaaTCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGCCAGGCCGTTTGAAATGAGGGTGTAATTGATTGGGCAACTGTGCGCCACGCTACTTTCTTCTTCGCTTAACGCGATa      >  2:1479216/1‑150 (MQ=255)
                                                                                                                                                 aaaTCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGCCAGGCCGTTTGAAATGAGGGTGTAATTGATTGGGCAACTGTGCGCCACGCTACTTTCTTCTTCGCTTAACGCGATa      <  1:193524/150‑1 (MQ=255)
                                                                                                                                                  aaTCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGCCAGGCCGTTTGAAATGAGGGTGTAATTGATTGGGCAACTGTGCGCCACGCTACTTTCTTCTTCGCTTAACGCGATAg     <  1:3073813/150‑1 (MQ=255)
                                                                                                                                                   aTCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGCCAGGCCGTTTGAAATGAGGGTGTAATTGATTGGGCAACTGTGCGCCACGCTACTTTCTTCTTCGCTTAACGCGATAgg    <  1:2713891/150‑1 (MQ=255)
                                                                                                                                                     cattcatCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGCCAGGCCGTTTGAAATGAGGGTGTAATTGATTGGGCAACTGTGCGCCACGCTACTTTCTTCTTCGCTTAACGCGATAGGAt  <  1:3612281/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTTAGCCACCCAGCCAATGAGCGGCAGCCCGTCGTTGGCGATCGCCTGAGCTGTTAGCAGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTTTCGGTCAGGTTTGCCAGGCCGTTTGAAATGAGGGTGTAATTGATTGGGCAACTGTGCGCCACGCTACTTTCTTCTTCGCTTAACGCGATAGGAT  >  CP009273/1660961‑1661259

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: