Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,483,912 T→A 16.2% V53D (GTC→GAC)  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,9120TA16.2% 369.9 / 59.3 153V53D (GTC→GAC) evgAresponse regulator in two‑component regulatory system with EvgS
Reads supporting (aligned to +/- strand):  ref base T (62/66);  new base A (12/13);  total (74/79)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.46e-01

Only 100 of 154 total aligned reads displayed.

TGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCAT  >  NC_000913/2483813‑2484009
                                                                                                   |                                                                                                 
tGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGt                                                                                                   <  1:2605832/100‑1 (MQ=255)
 gATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTc                                                                                                  >  1:3158629/1‑100 (MQ=255)
 gATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATAGCGTCATCATTGATGTc                                                                                                  >  1:2564570/1‑100 (MQ=255)
  aTCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCg                                                                                                 >  1:904260/1‑100 (MQ=255)
   tCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGa                                                                                                >  1:4300958/1‑100 (MQ=255)
   tCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGa                                                                                                >  1:1753164/1‑100 (MQ=255)
    cAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGat                                                                                               <  2:392587/100‑1 (MQ=255)
    cAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGat                                                                                               <  2:3905380/100‑1 (MQ=255)
    cAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGat                                                                                               >  2:3706776/1‑100 (MQ=255)
     aaaaaCGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGata                                                                                              >  2:322174/1‑100 (MQ=255)
     aaaaaCGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGACATCAATGATGACGata                                                                                              <  2:1794053/100‑1 (MQ=255)
     aaaaaCGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGACATCAATGATGACGata                                                                                              <  2:755761/100‑1 (MQ=255)
       aaaCGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATc                                                                                            >  2:3190550/1‑100 (MQ=255)
        aaCGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGACATCAATGATGACGATATca                                                                                           <  2:4168572/100‑2 (MQ=38)
          cGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATcccc                                                                                         >  1:1081009/1‑100 (MQ=255)
          cGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGACATCAATGATGACGATATcagg                                                                                         <  2:2558536/100‑4 (MQ=38)
          cGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGACATCAATGATGACGATATcagg                                                                                         <  2:4047716/100‑4 (MQ=38)
          cGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGACATCAATGATGACGATATcagg                                                                                         <  2:4148028/100‑4 (MQ=38)
          cGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGACATCAATGATGACGATATcagg                                                                                         <  1:2319583/100‑4 (MQ=38)
           gATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGACATCAATGATGACGATATcaggc                                                                                        <  2:2040633/100‑5 (MQ=38)
             tatTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGa                                                                                      <  2:904264/100‑1 (MQ=255)
              atTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATGCCCGGAg                                                                                     <  2:2303216/100‑1 (MQ=255)
              atTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAg                                                                                     >  1:1700536/1‑100 (MQ=255)
              atTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGACATCAATGATGACGATATcaggcttn                                                                                     <  1:2890267/100‑8 (MQ=37)
               ttGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGt                                                                                    >  1:519792/1‑100 (MQ=255)
                tGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGtt                                                                                   >  2:2760159/1‑100 (MQ=255)
                 gAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTa                                                                                  >  1:149981/1‑100 (MQ=255)
                   aaTCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAAc                                                                                >  2:1563715/1‑100 (MQ=255)
                   aaTCTTAGCAGAGTTGACAGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAAc                                                                                <  1:2412606/100‑1 (MQ=255)
                    aTCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACg                                                                               <  2:292636/100‑1 (MQ=255)
                    aTCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACg                                                                               >  2:696628/1‑100 (MQ=255)
                    aTCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACg                                                                               >  1:104453/1‑100 (MQ=255)
                    aTCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCCGAGTTAACg                                                                               >  1:803499/1‑100 (MQ=255)
                       ttAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACggca                                                                            <  2:4073087/100‑3 (MQ=255)
                        tAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTAt                                                                           >  2:2635313/1‑100 (MQ=255)
                         aGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATc                                                                          >  2:3861514/1‑100 (MQ=255)
                          gCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATcc                                                                         >  2:1385404/1‑100 (MQ=255)
                           cAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATcc                                                                         <  1:641282/99‑1 (MQ=255)
                           cAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCa                                                                        <  1:2651783/100‑1 (MQ=255)
                            agagTTGACTGAAGGCGGAAGGGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAg                                                                       <  2:1760037/100‑1 (MQ=255)
                             gagTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAgg                                                                      <  2:409314/100‑1 (MQ=255)
                               gTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGgtg                                                                    >  2:2236920/1‑100 (MQ=255)
                               gTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGgtg                                                                    >  1:3424583/1‑100 (MQ=255)
                                  gACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTa                                                                 >  2:1616198/1‑100 (MQ=255)
                                   aCTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAg                                                                <  1:2156343/100‑1 (MQ=255)
                                    cTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGa                                                               <  1:829494/100‑1 (MQ=255)
                                      gAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGaaa                                                             >  1:2909347/1‑100 (MQ=255)
                                        aGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAACCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACg                                                           >  2:3104706/1‑100 (MQ=255)
                                         ggCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGc                                                          >  1:439195/1‑100 (MQ=255)
                                         ggCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATACAGGTGTTAGAAACGc                                                          >  1:1768275/1‑100 (MQ=255)
                                           cGGAAGTGCCGTTCAGCGGGTGGAAACACTTAATCCTGATATCGTCATCATTGATGTTGATATCCCCGGAGTTAACGGTATCCAGTTGTTAGAAACGCTg                                                        <  1:1620631/100‑1 (MQ=255)
                                           cGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTg                                                        <  2:480916/100‑1 (MQ=255)
                                           cGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTg                                                        <  1:720492/100‑1 (MQ=255)
                                             gAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAg                                                      <  2:1186482/100‑1 (MQ=255)
                                             gAAGTGCAGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAg                                                      <  2:149982/100‑1 (MQ=255)
                                                gTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGaa                                                   <  2:1567984/100‑1 (MQ=255)
                                                gTGCCGTTCAGCGGGTGGAAACACGTAAGCCTGAAATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCtn                                                        <  1:3775940/95‑2 (MQ=255)
                                                  gCCGTTCAGCGGTTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGaacggca                                                 >  2:2891245/1‑94 (MQ=255)
                                                  gCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAgc                                                 >  2:708766/1‑100 (MQ=255)
                                                   ccGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAgcg                                                >  2:1042338/1‑100 (MQ=255)
                                                   ccGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAgcg                                                >  2:2822489/1‑100 (MQ=255)
                                                     gTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGcc                                              >  1:1289209/1‑100 (MQ=255)
                                                      ttCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCa                                             <  1:3197545/100‑1 (MQ=255)
                                                         aGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAata                                          <  2:3438934/100‑1 (MQ=255)
                                                            gggTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAg                                       <  1:777444/100‑1 (MQ=255)
                                                              gTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCg                                     <  1:2314961/100‑1 (MQ=255)
                                                               tGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCgg                                    <  1:2035112/100‑1 (MQ=255)
                                                                  aaaCACGTAAGCCTGATATCGGTAGCGGTGATGTCGATGTCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAat                                 >  1:3020822/1‑100 (MQ=255)
                                                                   aaCACTTAAGCCTGATATCGTCATCATTGATGTCGATATCACCGGAGTTAACGGGATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAat                                 >  1:1895150/1‑99 (MQ=255)
                                                                    acacTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAatta                               <  1:4278153/100‑1 (MQ=255)
                                                                    acacTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAatta                               <  2:2309655/100‑1 (MQ=255)
                                                                     cacTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAattat                              >  1:2229859/1‑100 (MQ=255)
                                                                       cTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAattatta                            <  1:848507/100‑1 (MQ=255)
                                                                       cTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAattatta                            <  2:1362757/100‑1 (MQ=255)
                                                                       cTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAattatta                            <  2:3424601/100‑1 (MQ=255)
                                                                       cTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAattatta                            <  1:8976/100‑1 (MQ=255)
                                                                       cTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAattatta                            <  1:848570/100‑1 (MQ=255)
                                                                         tAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAattattatt                          <  2:1745209/100‑1 (MQ=255)
                                                                          tccggggATATCGACATCAATGATGACGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGCAttattatta                         >  1:1712964/7‑100 (MQ=37)
                                                                            gCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATc                       <  2:442163/100‑1 (MQ=255)
                                                                            gCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATc                       >  2:1084869/1‑100 (MQ=255)
                                                                            cggggATATCGACATCAATGATGACGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATc                       >  2:4369966/5‑100 (MQ=38)
                                                                            cggggATATCGACATCAATGATGACGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATc                       >  2:3654703/5‑100 (MQ=38)
                                                                              gggATATCGACATCAATGATGACGAAATCCACGGAGTTAACGGTATCCAGATGTTAGAAACACTCAGGAAGCGCCAATATAGCGGAATtattattat                        >  1:316338/3‑97 (MQ=37)
                                                                              cTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCg                      <  1:2259019/99‑1 (MQ=255)
                                                                              cTGATATCGTCATAATTGATGTCGATATCCCCGGAGTTAACGGTATTCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCg                      <  1:1370710/99‑1 (MQ=255)
                                                                               tGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGtc                    >  2:4363840/1‑100 (MQ=255)
                                                                               ggATATCGACATCAATGATGACGATATCCCCGGAGTTAACGGCATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATCGCGGAAATATTATTATc                       >  1:2381082/2‑97 (MQ=37)
                                                                                gATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGtct                   >  1:750910/1‑100 (MQ=255)
                                                                                gATATCGACATCAATGATGACGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGtct                   >  2:2171288/1‑100 (MQ=39)
                                                                                  tatCGACATCAATGATGACGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTcc                 >  2:774417/1‑100 (MQ=255)
                                                                                   atCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCg                <  2:2417804/100‑1 (MQ=255)
                                                                                    tCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGc               <  2:1309102/100‑1 (MQ=255)
                                                                                    tCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGc               <  1:3892577/100‑1 (MQ=255)
                                                                                      gTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTa             >  1:2079793/1‑100 (MQ=255)
                                                                                           catTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAAt        >  2:3025306/1‑100 (MQ=255)
                                                                                            atTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATg       <  2:3004513/100‑1 (MQ=255)
                                                                                               gATGTCGATATCCCCGGGGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGAcc    >  1:656647/1‑100 (MQ=255)
                                                                                               gATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGAcc    >  2:2432547/1‑100 (MQ=255)
                                                                                                aTGTCAACATCCCCGGACTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCa   <  2:3943782/100‑1 (MQ=255)
                                                                                                   |                                                                                                 
TGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCAT  >  NC_000913/2483813‑2484009

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