Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,181,277 T→A 5.7% V61E (GTG→GAG)  yegV → putative kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,181,2770TA 5.7% 251.6 / 4.3 88V61E (GTG→GAG) yegVputative kinase
Reads supporting (aligned to +/- strand):  ref base T (43/39);  new base A (2/3);  total (46/42)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.80e-01

TGAATTAACCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGC  >  NC_000913/2181128‑2181424
                                                                                                                                                     |                                                                                                                                                   
tGAATTAACCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGt                                                                                                                                                     >  2:1407340/1‑150 (MQ=255)
 gAATTAACCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTg                                                                                                                                                    >  2:1507657/1‑150 (MQ=255)
  aaTTTACAACGCGTCCGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGAGTTAACGTTGGCGGCTTCGCACTGAATATTGCCGTgg                                                                                                                                                   <  2:1278248/150‑1 (MQ=255)
   aTTAACCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGc                                                                                                                                                  <  2:1114010/150‑1 (MQ=255)
   aTTAACCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGc                                                                                                                                                  <  2:199950/150‑1 (MQ=255)
       aCCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCAGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTa                                                                                                                                              >  1:489471/1‑150 (MQ=255)
           cggttCAGTCGGTGAGGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCACTCCCCTGGCTTGTGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAgc                                                                                                                                          <  2:79843/146‑1 (MQ=255)
           cgcgTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAgc                                                                                                                                          <  2:1160078/150‑1 (MQ=255)
           cgcgTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAgc                                                                                                                                          <  2:40087/150‑1 (MQ=255)
             cgTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAgcgc                                                                                                                                        >  1:727163/1‑150 (MQ=255)
              gTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGcc                                                                                                                                       <  2:902533/150‑1 (MQ=255)
               tCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCt                                                                                                                                      <  2:12817/150‑1 (MQ=255)
                  gTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCgg                                                                                                                                   >  1:37464/1‑150 (MQ=255)
                    cGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGAGGCGCTTTAACGCCACGGCa                                                                                                                                 <  2:57599/150‑1 (MQ=255)
                     ggTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCAt                                                                                                                                >  1:1330613/1‑150 (MQ=255)
                     ggTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGAGGCGCTTTAACGCCACGGCaa                                                                                                                                <  1:882459/150‑2 (MQ=255)
                     ggTGATGGGCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGAGGCGCTTTAACGCCACGGCaa                                                                                                                                <  1:271651/150‑2 (MQ=255)
                      gTGATGGTCGTCGGCGCGGCGGTCATTGCCGTGATGGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGACGTGGCGTTAAAGCGCCTCGGCATc                                                                                                                               >  1:118011/1‑150 (MQ=255)
                      gTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCGTc                                                                                                                               >  1:1222286/1‑150 (MQ=255)
                      gTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATc                                                                                                                               >  2:1481467/1‑150 (MQ=255)
                         aTGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCACCGaa                                                                                                                            <  2:1336229/150‑1 (MQ=255)
                          tGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAg                                                                                                                           <  1:1630567/150‑1 (MQ=255)
                          tGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAg                                                                                                                           >  1:1213600/1‑150 (MQ=255)
                               gtcgGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGt                                                                                                                      >  2:1458916/1‑150 (MQ=255)
                                tcgGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCATCCCCTGGCGTGGGTGCGATATCGAACTGAATCAGCAGAGCGGTAACGTTGGCGGCTGCGCACTGAATATTGCTGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTa                                                                                                                     >  2:1077573/1‑150 (MQ=255)
                                 cggcgcggcgGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTaa                                                                                                                    <  2:15038/150‑1 (MQ=255)
                                   gcgcggcgGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATg                                                                                                                  <  1:307705/150‑1 (MQ=255)
                                      cggcggTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCt                                                                                                               >  2:658112/1‑150 (MQ=255)
                                        gcggTCATTGACGTGATCGCCGACGCTTATGACCTCCCCTGGCGTGGGTGCGATATCGATCTGAAACAGCAGAGCGTTAACGTTGCCGCCTGCGCACTCAATATTGCCGTGCCGTTAAAGCGCCTCGGCATCGAACCGGGAAAAGCCTTg                                                                                                             >  2:1555411/1‑150 (MQ=255)
                                        gcggGCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTg                                                                                                             <  2:117528/150‑1 (MQ=255)
                                          ggTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGc                                                                                                                                  <  2:887796/127‑1 (MQ=255)
                                          ggTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGc                                                                                                                                  >  1:887796/1‑127 (MQ=255)
                                          ggTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCAGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGcc                                                                                                           >  2:245009/1‑150 (MQ=255)
                                           gTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCg                                                                                                          >  1:452812/1‑150 (MQ=255)
                                                tGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCgg                                                                                                     <  2:1330613/150‑1 (MQ=255)
                                                tGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCGCGGCATCGAAGCGGGTAATGCCTTGCCGCTCgg                                                                                                     <  2:138183/150‑1 (MQ=255)
                                                     tGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAgg                                                                                                >  1:13570/1‑150 (MQ=255)
                                                      gATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGACACAGCAGAGCGTTAACGTTGGCGGCGGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAggg                                                                                               <  1:213047/150‑1 (MQ=255)
                                                            cGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGtgtg                                                                                         >  2:996872/1‑150 (MQ=255)
                                                              aCGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTggg                                                                                       <  2:727163/150‑1 (MQ=255)
                                                                 cTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCgg                                                                                    >  1:546229/1‑150 (MQ=255)
                                                                   tATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGgag                                                                                  >  1:1036891/1‑150 (MQ=255)
                                                                   tATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGgag                                                                                  >  2:1598779/1‑150 (MQ=255)
                                                                    aTGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGgaga                                                                                 >  2:1591244/1‑150 (MQ=255)
                                                                    aTGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGgaga                                                                                 >  1:1181209/1‑150 (MQ=255)
                                                                     tGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAgat                                                                                >  1:1101847/1‑150 (MQ=255)
                                                                      gCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAgatg                                                                               <  1:584043/150‑1 (MQ=255)
                                                                      gCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAgatg                                                                               <  1:1407340/150‑1 (MQ=255)
                                                                        ccTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAgatgat                                                                             <  1:677604/150‑1 (MQ=255)
                                                                         cTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGAtt                                                                            <  1:1507657/150‑1 (MQ=255)
                                                                           ccccTGGCGTGGGTGCGATATCGAAGTGAAACAGCAGAGCGTTAAGGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCg                                                                          <  1:111364/150‑1 (MQ=255)
                                                                           ccccTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCg                                                                          <  1:128725/150‑1 (MQ=255)
                                                                           ccccTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCg                                                                          >  2:1704208/1‑150 (MQ=255)
                                                                           ccccTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCg                                                                          >  1:720540/1‑150 (MQ=255)
                                                                           ccccTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCg                                                                          >  2:961591/1‑150 (MQ=255)
                                                                            cccTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGc                                                                         <  1:658112/150‑1 (MQ=255)
                                                                             ccTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCa                                                                        <  2:452812/150‑1 (MQ=255)
                                                                               tGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAAc                                                                      <  1:1403871/150‑1 (MQ=255)
                                                                                ggCGGGGGGGCGATATCGAACTGAAACAGCAGAGCGTTACCGTTGGCGGCTGCGCACTGAATATTGCCGGGGCGTTAAAGCGCCTCGGCATCGACGCGGGTAGTGCATTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAAcc                                                                     >  1:1611729/1‑150 (MQ=255)
                                                                                  cGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCgg                                                                   >  2:1325766/1‑150 (MQ=255)
                                                                                   gTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGa                                                                  <  2:1181209/150‑1 (MQ=255)
                                                                                       gTGCGATATCGAACTGAAACAGCAGAGCGTTACCGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCCACCGGATGGc                                                              >  2:1059593/1‑150 (MQ=255)
                                                                                             tatCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAaga                                                        <  2:720540/150‑1 (MQ=255)
                                                                                             tatCGAACTGAAACAGCAGAGCGTTAACGTTGGAGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCCGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTCGGCGGAGATGATTCGCAACCGGATGGCAATaga                                                        >  1:275651/1‑150 (MQ=255)
                                                                                               tCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAgg                                                      >  2:365203/1‑150 (MQ=255)
                                                                                               tCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAgg                                                      <  2:1428386/150‑1 (MQ=255)
                                                                                                cGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAggg                                                     <  2:37464/150‑1 (MQ=255)
                                                                                                 gAACTGTATCAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAATTCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGTATGGCAAAAGAGGGc                                                    <  2:682722/150‑1 (MQ=255)
                                                                                                  aaCTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTTCCGTGTCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATTCCTAGCCGCTCGGTCAGGGCGTGTGTGCGGACATGATTCGCAACCGGATGGCAAAAGAGGCCt                                                   >  2:598434/1‑150 (MQ=255)
                                                                                                  aaCTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCt                                                   <  2:1055787/150‑1 (MQ=255)
                                                                                                     ggcaaCAGCAGAGCTGTAACGTTGGCGGCGGCGCACTGAATTTTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTaa                                                <  1:339798/147‑1 (MQ=255)
                                                                                                      gAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGACGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAAt                                               >  1:33628/1‑150 (MQ=255)
                                                                                                        aaCAGCAGAGCGTTAACGTTTGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCa                                             >  2:147320/1‑150 (MQ=255)
                                                                                                                agCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATc                                     <  1:1201097/150‑1 (MQ=255)
                                                                                                                   gTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGAt                                  >  1:554725/1‑150 (MQ=255)
                                                                                                                      aaCGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTAGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAAc                               >  1:1480453/1‑150 (MQ=255)
                                                                                                                             gcggcTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAg                        >  2:571112/1‑150 (MQ=255)
                                                                                                                                gcTGCGCACTGAATATTGCCGTGGCGTTAACGCGCCTCGGCATAGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCTAACGCGGGCTTACTCAGCCTGATCGATTGCGCGGTAGGTg                     >  1:131051/1‑150 (MQ=255)
                                                                                                                                   gcgcACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATa                  >  2:86139/1‑150 (MQ=255)
                                                                                                                                     gcACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAAc                <  1:1598779/150‑1 (MQ=255)
                                                                                                                                     gcACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAAc                <  1:1692795/150‑1 (MQ=255)
                                                                                                                                          gAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTg           >  1:1582281/1‑150 (MQ=255)
                                                                                                                                             tatTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGgtg        >  2:869926/1‑150 (MQ=255)
                                                                                                                                              gatGCCGAGGCGCTTTAACGCCACGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGgtgt       >  1:1427418/3‑150 (MQ=255)
                                                                                                                                               atGCCGAGGCGCTTTAACGCCACGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTc      >  2:696749/2‑150 (MQ=255)
                                                                                                                                                  ccGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGt                      <  2:956681/131‑1 (MQ=255)
                                                                                                                                                  ccGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGt                      >  1:956681/1‑131 (MQ=255)
                                                                                                                                                  ccGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTgg   >  2:1331540/1‑150 (MQ=255)
                                                                                                                                                   cGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGgc  <  2:13570/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TGAATTAACCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGC  >  NC_000913/2181128‑2181424

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑