Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,181,284 A→T 5.6% L63F (TTA→TTT) ‡ yegV → putative kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,181,2840AT 5.6% 262.4 / 4.2 89L63F (TTA→TTT) ‡yegVputative kinase
Reads supporting (aligned to +/- strand):  ref base A (43/41);  new base T (2/3);  total (45/44)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.48e-01

ACCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTG  >  NC_000913/2181135‑2181432
                                                                                                                                                     |                                                                                                                                                    
aCCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCAGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTa                                                                                                                                                      >  1:489471/1‑150 (MQ=255)
    cggttCAGTCGGTGAGGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCACTCCCCTGGCTTGTGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAgc                                                                                                                                                  <  2:79843/146‑1 (MQ=255)
    cgcgTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAgc                                                                                                                                                  <  2:40087/150‑1 (MQ=255)
    cgcgTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAgc                                                                                                                                                  <  2:1160078/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:1213600/1‑150 (MQ=255)
                   tGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAg                                                                                                                                   <  1:1630567/150‑1 (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                                                                                                                                          >  1:887796/1‑127 (MQ=255)
                                   ggTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGc                                                                                                                                          <  2:887796/127‑1 (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                                                                                          >  2:1598779/1‑150 (MQ=255)
                                                            tATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGgag                                                                                          >  1:1036891/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:1407340/150‑1 (MQ=255)
                                                               gCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAgatg                                                                                       <  1:584043/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                                                                                  >  2:961591/1‑150 (MQ=255)
                                                                    ccccTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCg                                                                                  <  1:128725/150‑1 (MQ=255)
                                                                    ccccTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCg                                                                                  >  2:1704208/1‑150 (MQ=255)
                                                                    ccccTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCg                                                                                  >  1:720540/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:1428386/150‑1 (MQ=255)
                                                                                        tCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAgg                                                              >  2:365203/1‑150 (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:1692795/150‑1 (MQ=255)
                                                                                                                              gcACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAAc                        <  1:1598779/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                              >  1:956681/1‑131 (MQ=255)
                                                                                                                                           ccGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGt                              <  2:956681/131‑1 (MQ=255)
                                                                                                                                           ccGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTgg           >  2:1331540/1‑150 (MQ=255)
                                                                                                                                            cGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGgc          <  2:13570/150‑1 (MQ=255)
                                                                                                                                               ggCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCt       <  1:958196/150‑1 (MQ=255)
                                                                                                                                                gCGTTAATGCGCCTCGGCATCGATGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTg      <  1:598434/150‑1 (MQ=255)
                                                                                                                                                gCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTg      <  2:799740/150‑1 (MQ=255)
                                                                                                                                                  gTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGt    <  2:33628/150‑1 (MQ=255)
                                                                                                                                                   ttAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGtt   >  1:134532/1‑150 (MQ=255)
                                                                                                                                                    tAAAGCGCCTCGGCATCGAATCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTg  <  1:570911/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ACCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTG  >  NC_000913/2181135‑2181432

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: