Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,483,906 T→A 15.6% I51N (ATT→AAT)  evgA → response regulator in two‑component regulatory system with EvgS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,483,9060TA15.6% 354.8 / 52.7 135I51N (ATT→AAT) evgAresponse regulator in two‑component regulatory system with EvgS
Reads supporting (aligned to +/- strand):  ref base T (56/58);  new base A (8/13);  total (64/71)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.05e-01

Only 100 of 139 total aligned reads displayed.

TATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATG  >  NC_000913/2483810‑2484004
                                                                                                |                                                                                                  
tATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGa                                                                                                 <  2:512201/100‑1 (MQ=255)
 aTTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGAt                                                                                                <  1:2703639/100‑1 (MQ=255)
   tGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCTTGATATCGTCATCATTGATGt                                                                                              >  2:3009717/1‑100 (MQ=255)
    gATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTc                                                                                             <  2:157370/100‑1 (MQ=255)
    gATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTc                                                                                             >  2:3256198/1‑100 (MQ=255)
    gATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTc                                                                                             >  2:737188/1‑100 (MQ=255)
    gATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTc                                                                                             >  1:1297384/1‑100 (MQ=255)
     aTCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCg                                                                                            >  2:2108472/1‑100 (MQ=255)
       cAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGat                                                                                          >  2:188023/1‑100 (MQ=255)
       cAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGat                                                                                          >  1:1426679/1‑100 (MQ=255)
       cAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGat                                                                                          >  1:2917218/1‑100 (MQ=255)
       cAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGat                                                                                          <  2:3761505/100‑1 (MQ=255)
       cAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGACATCAATGATGACGat                                                                                          <  1:3373125/100‑1 (MQ=255)
        aaaaaaGATATTGAAATCTTAGCAGAGATGCCTGAAGGCGGAAGTCCCGTTCAGCGGGTGGAAACACTTAAGCCTGATCCCGTCACCATTGATGTCGata                                                                                         <  1:176695/100‑1 (MQ=255)
        aaaaaCGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGata                                                                                         <  2:1638454/100‑1 (MQ=255)
        aaaaaCGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGata                                                                                         >  2:2914925/1‑100 (MQ=255)
         aaaaCGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGGGTCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGACATCAATGATGACGa                                                                                           <  2:1724899/97‑1 (MQ=255)
          aaaCGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGGCGatgt                                                                                        >  1:3512988/1‑97 (MQ=255)
           aaCGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATcc                                                                                      >  2:1853647/1‑100 (MQ=255)
            aCGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATccc                                                                                     >  2:3577159/1‑100 (MQ=255)
            aCGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGACATCAATGATGACGATATcag                                                                                     <  2:2702669/100‑3 (MQ=38)
            aCGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGACATCAATGATGACGATATcag                                                                                     <  2:1982802/100‑3 (MQ=38)
             cGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATcccc                                                                                    >  2:2923033/1‑100 (MQ=255)
             cGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGACATCAATGATGACGATATcagg                                                                                    <  1:419637/100‑4 (MQ=38)
              gATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGACATCAATGATGACGATATcaggc                                                                                   <  1:1252784/100‑5 (MQ=38)
              gATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGACATCAATGATGACGATATcaggc                                                                                   <  1:3239348/100‑5 (MQ=38)
              gATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGACATCAATGATGACGATATcaggc                                                                                   <  1:1947695/100‑5 (MQ=38)
                 atTGAAATCTTAGCAGATTTGACTGAAGGCGTAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGACATCAATGATTACGATATcaggc                                                                                   <  1:803613/97‑5 (MQ=37)
                  ttGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGt                                                                               <  2:2186966/100‑1 (MQ=255)
                   tGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAg                                                                                >  2:1726280/1‑98 (MQ=255)
                   tGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGtt                                                                              <  1:1317088/100‑1 (MQ=255)
                   tGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGACATCAATGATGACGATATcaggcttaag                                                                              <  2:3845937/100‑10 (MQ=25)
                    gAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTa                                                                             <  1:250960/100‑1 (MQ=255)
                      aaTCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAAc                                                                           >  2:5971/1‑100 (MQ=255)
                       aTCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACg                                                                          >  2:2380765/1‑100 (MQ=255)
                        tCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACgn                                                                         <  2:3342257/100‑2 (MQ=255)
                        tCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACgg                                                                         <  2:1297392/100‑1 (MQ=255)
                        tCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACgg                                                                         <  2:1628456/100‑1 (MQ=255)
                         cTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGt                                                                        <  2:622776/100‑1 (MQ=255)
                            aGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAAcagatcg                                                                     >  1:2625717/1‑94 (MQ=255)
                             gCAGAGTTGAATGAAGGCGGTAGTGCCGTTGAGCGGGTGGAAGCCATTAAGGCTGATATCGTCATCATTGATGTCGATACCCCCGGAGTTAACGGTATc                                                                     <  1:3233817/99‑1 (MQ=255)
                              cAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCa                                                                   <  2:48133/100‑1 (MQ=255)
                               agagTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCa                                                                   <  1:2526772/99‑1 (MQ=255)
                               agagTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAg                                                                  >  1:354674/1‑100 (MQ=255)
                               agagTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCACCGGAGTTAACGGTATCCAg                                                                  >  1:2355668/1‑100 (MQ=255)
                                gagTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAgg                                                                 >  2:2095844/1‑100 (MQ=255)
                                 agTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGCAAACACTTAAGCCTGATATCGTCATCATTGCTGTCGATATCCCCGGAGTTAACGGTATCCAGgt                                                                >  1:3012808/1‑100 (MQ=255)
                                  gTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGgtg                                                               >  2:1811152/1‑100 (MQ=255)
                                    tGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGtt                                                             <  1:122322/100‑1 (MQ=255)
                                      aCTGAAGGCGGAAGTGCCGTTCAGGGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAg                                                           <  1:1718745/100‑1 (MQ=255)
                                      aCTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAg                                                           <  2:1450407/100‑1 (MQ=255)
                                       cTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGa                                                          >  1:2570754/1‑100 (MQ=255)
                                       cTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGa                                                          >  2:2322150/1‑100 (MQ=255)
                                       cTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGa                                                          >  2:118022/1‑100 (MQ=255)
                                       cTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGa                                                          <  1:3780911/100‑1 (MQ=255)
                                        ngAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGaa                                                         >  2:3057343/2‑100 (MQ=255)
                                         gAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGaaa                                                        >  1:3126991/1‑100 (MQ=255)
                                              cGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTg                                                   <  2:2103142/100‑1 (MQ=255)
                                                gAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAg                                                 >  2:2693040/1‑100 (MQ=255)
                                                gAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAg                                                 <  2:546904/100‑1 (MQ=255)
                                                    tGCCGTGCAGCGGGTGGAAACACTTAAGCCTGATATCGTCGTCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAg                                             <  1:3706119/100‑1 (MQ=255)
                                                     gCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAgc                                            <  1:737182/100‑1 (MQ=255)
                                                     gCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAgc                                            >  1:3230870/1‑100 (MQ=255)
                                                       cGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAgcgc                                          >  2:1696892/1‑100 (MQ=255)
                                                        gTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGcc                                         <  1:1568528/100‑1 (MQ=255)
                                                         ttCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCa                                        <  2:1987731/100‑1 (MQ=255)
                                                          tCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCaa                                       <  2:2065496/100‑1 (MQ=255)
                                                           cAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAat                                      <  2:1182797/100‑1 (MQ=255)
                                                            aGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAata                                     <  2:3230895/100‑1 (MQ=255)
                                                             gCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAatat                                    <  2:3440783/100‑1 (MQ=255)
                                                              cGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAAtata                                   <  2:61525/100‑1 (MQ=255)
                                                                ggTGGAAACACTTAAGCCTGATATCGCAATTATTGATGTCGATATCCCCGTAGTTAACGGTAGCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGc                                 <  1:538535/100‑1 (MQ=255)
                                                                   ggAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGa                              >  2:2582301/1‑100 (MQ=255)
                                                                   ggAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGa                              <  2:1328095/100‑1 (MQ=255)
                                                                    gAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGaa                             >  2:217986/1‑100 (MQ=255)
                                                                     aaaCACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGa                              <  2:1124544/98‑1 (MQ=255)
                                                                      aaCACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAatt                           <  1:1845429/100‑1 (MQ=255)
                                                                      aaCACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAat                            >  1:1968048/1‑99 (MQ=255)
                                                                         acTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAattatt                        >  2:8619/1‑100 (MQ=255)
                                                                         acTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAattatt                        <  2:1301190/100‑1 (MQ=255)
                                                                         acTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAATTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAatgatt                        >  2:2583415/1‑100 (MQ=255)
                                                                          cTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAattatta                       >  1:1759216/1‑100 (MQ=255)
                                                                           ttAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAattattat                      <  2:1052480/100‑1 (MQ=255)
                                                                              aGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATtattattat                   <  2:1913862/100‑1 (MQ=255)
                                                                              aGCCTGATATCGTAATCATTTATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGATGAAGCGCCAATATAGCGGAATtattattat                   >  1:221236/1‑100 (MQ=255)
                                                                               gCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATc                  >  2:638031/1‑100 (MQ=255)
                                                                               gCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATc                  >  2:1943953/1‑100 (MQ=255)
                                                                                ccTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCg                 <  1:3110885/100‑1 (MQ=255)
                                                                                ccTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCg                 <  2:311495/100‑1 (MQ=255)
                                                                                   gATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGCATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGtct              >  2:270120/1‑100 (MQ=255)
                                                                                     tatCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTcc            <  1:2108457/100‑1 (MQ=255)
                                                                                     tatCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTcc            <  1:3041045/100‑1 (MQ=255)
                                                                                     tatCGACATCAATGATGACGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTcc            >  2:2540219/1‑100 (MQ=255)
                                                                                      atCGACATCAATGATGACGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCg           >  1:1122118/1‑100 (MQ=255)
                                                                                       tCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGc          <  2:112461/100‑1 (MQ=255)
                                                                                       tCGACATCAATGATGACGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGc          >  2:3515316/1‑100 (MQ=255)
                                                                                              catTGATGTCGATATCCCCGGAGTTAACGGTATCCCGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTaaaact   >  1:345696/1‑98 (MQ=255)
                                                                                              catTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAAt   >  2:3708288/1‑100 (MQ=255)
                                                                                              catTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAAt   >  2:3346561/1‑100 (MQ=255)
                                                                                               atTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATg  <  2:3260287/100‑1 (MQ=255)
                                                                                                |                                                                                                  
TATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATG  >  NC_000913/2483810‑2484004

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: