Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 332,019 T→C V142A (GTG→GCG)  yahA → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913332,0190TC99.6% 822.0 / NA 255V142A (GTG→GCG) yahAhypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (114/140);  minor base A (1/0);  total (115/140)
Fisher's exact test for biased strand distribution p-value = 4.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 263 total aligned reads displayed.

CCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAA  >  NC_000913/331929‑332108
                                                                                          |                                                                                         
ccATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGgcgc                                                                                         <  1:2864945/93‑1 (MQ=255)
 cATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCt                                                                                        <  1:1796829/93‑1 (MQ=255)
  aTCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCtt                                                                                       >  1:5450846/1‑93 (MQ=255)
   tCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTg                                                                                      <  1:1982505/93‑1 (MQ=255)
    cAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGt                                                                                     <  1:2904462/93‑1 (MQ=255)
     aGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTc                                                                                    >  1:787771/1‑93 (MQ=255)
     aGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTc                                                                                    <  2:3669770/93‑1 (MQ=255)
     aGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTc                                                                                    <  1:3171060/93‑1 (MQ=255)
      gTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTcc                                                                                   <  1:6299881/93‑1 (MQ=255)
        cTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGc                                                                                 <  2:4631603/93‑1 (MQ=255)
        cTGGCGTTGGAAAACCACGAATCCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGc                                                                                 >  2:6214220/1‑93 (MQ=255)
         tGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCt                                                                                >  2:524211/1‑93 (MQ=255)
             gTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGa                                                                            >  1:1707494/1‑93 (MQ=255)
              ttGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGaa                                                                           >  1:2160717/1‑93 (MQ=255)
                ggAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACa                                                                         >  1:161747/1‑93 (MQ=255)
                  aaaaCCACGAATTCAAACCGCGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGGACATc                                                                       >  2:2469690/1‑93 (MQ=255)
                  aaaaCCACGAATGCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATc                                                                       <  1:2543453/93‑1 (MQ=255)
                    aaCCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATGca                                                                     >  2:2545991/1‑93 (MQ=255)
                    aaCCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCca                                                                     <  1:4890690/93‑1 (MQ=255)
                    aaCCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCca                                                                     <  1:4919954/93‑1 (MQ=255)
                     aCCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCcac                                                                    >  1:3576258/1‑93 (MQ=255)
                       cACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACaa                                                                  >  1:467140/1‑93 (MQ=255)
                       cACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACaa                                                                  >  1:6483086/1‑93 (MQ=255)
                        aCGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACaaa                                                                 <  2:5677994/93‑1 (MQ=255)
                        aCGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACaaa                                                                 <  1:5027039/93‑1 (MQ=255)
                         cGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAAc                                                                <  2:4442235/93‑1 (MQ=255)
                         cGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGAGCTTGTCCGCTGGGAACATCCACAAAc                                                                >  1:2610545/1‑93 (MQ=255)
                          gAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAACCg                                                               <  1:213836/93‑1 (MQ=255)
                            aTTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACggg                                                             <  2:2132931/93‑1 (MQ=255)
                              tCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGaa                                                           <  1:1399879/93‑1 (MQ=255)
                              tCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGaa                                                           >  2:4564901/1‑93 (MQ=255)
                                aaaCCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAAtt                                                         <  1:2803751/93‑1 (MQ=255)
                                aaaCCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAAtt                                                         >  2:2865975/1‑93 (MQ=255)
                                aaaCCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAAtt                                                         <  1:4069755/93‑1 (MQ=255)
                                aaaCCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAAtt                                                         >  2:6105895/1‑93 (MQ=255)
                                aaaCCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAAtt                                                         >  1:5775088/1‑93 (MQ=255)
                                 aaCCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTa                                                        >  2:1725860/1‑93 (MQ=255)
                                  aCCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGGACATCCACAAACGGGAATTAt                                                       >  1:3825540/1‑93 (MQ=255)
                                  aCCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTAt                                                       <  1:3083928/93‑1 (MQ=255)
                                   ccGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATc                                                      <  2:4915344/93‑1 (MQ=255)
                                    cGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATcc                                                     <  2:1006236/93‑1 (MQ=255)
                                     gTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATccc                                                    <  1:274556/93‑1 (MQ=255)
                                       ggATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCAc                                                  >  2:1173978/1‑93 (MQ=255)
                                       ggATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCAc                                                  <  1:2231782/93‑1 (MQ=255)
                                       ggATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCAc                                                  >  2:1436143/1‑93 (MQ=255)
                                        gATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCAcc                                                 >  1:1939246/1‑93 (MQ=255)
                                        gATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCAcc                                                 >  1:4280524/1‑93 (MQ=255)
                                         aTCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATGCACAAACGGGAATTATCCCTCCg                                                >  2:5650599/1‑93 (MQ=255)
                                         aTCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCg                                                >  2:3871183/1‑93 (MQ=255)
                                           ccAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGa                                              <  2:3358845/93‑1 (MQ=255)
                                            cAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGac                                             >  1:3755843/1‑92 (MQ=255)
                                            cAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGAt                                             <  2:5368591/93‑1 (MQ=255)
                                             aaCCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATc                                            >  1:6017730/1‑93 (MQ=255)
                                               ccGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAg                                          <  2:6622596/93‑1 (MQ=255)
                                                cGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGGATTATCCCACCGGATCAGt                                         >  1:3021480/1‑93 (MQ=255)
                                                 ggTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGtt                                        <  2:1058016/93‑1 (MQ=255)
                                                 ggTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGtt                                        <  1:1263745/93‑1 (MQ=255)
                                                  gTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCCGttt                                       >  2:3231273/1‑93 (MQ=255)
                                                   ttttCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTa                                      <  1:1163425/93‑1 (MQ=255)
                                                     ttCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTAtt                                    <  2:744673/93‑1 (MQ=255)
                                                     ttCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTAtt                                    >  1:100879/1‑93 (MQ=255)
                                                     ttCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTAtt                                    <  1:3448923/93‑1 (MQ=255)
                                                     ttCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTAtt                                    <  1:6634670/93‑1 (MQ=255)
                                                       cTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTcc                                  <  2:180306/93‑1 (MQ=255)
                                                       cTGCGCGCAGACTGGCGTACTGACGGGCTGGGAGGCGCTTGTCCGCTGGTAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTcc                                  <  2:2139317/93‑1 (MQ=255)
                                                        tGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCct                                 <  1:4438651/93‑1 (MQ=255)
                                                         gcgcgcAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCctc                                <  2:3000814/93‑1 (MQ=255)
                                                          cgcgcAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTGTCCCACCGGATCAGTTTATTCctct                               <  2:2910092/93‑1 (MQ=255)
                                                          cgcgcAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCctct                               <  1:1911784/93‑1 (MQ=255)
                                                          cgcgcAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCctct                               <  1:5092788/93‑1 (MQ=255)
                                                           gcgcAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTg                              <  1:5717059/93‑1 (MQ=255)
                                                            cgcAGCCTGCCGTACTGACGGGCTGTGAGGCGCCTGCCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTgg                             <  2:1271648/93‑1 (MQ=255)
                                                            cgcAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTgg                             >  1:6155176/1‑93 (MQ=255)
                                                             gcAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGc                            <  1:5282081/93‑1 (MQ=255)
                                                              cAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCg                           >  2:5508759/1‑93 (MQ=255)
                                                               aGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCgg                          >  1:1568153/1‑93 (MQ=255)
                                                                gACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGa                         >  1:2083225/1‑93 (MQ=255)
                                                                   tGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTc                      <  2:4379731/93‑1 (MQ=255)
                                                                    ggCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCa                     >  2:4959102/1‑93 (MQ=255)
                                                                      cGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATc                   <  2:4232611/93‑1 (MQ=255)
                                                                       gTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATcc                  <  2:6703163/93‑1 (MQ=255)
                                                                       gTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATcc                  <  1:986216/93‑1 (MQ=255)
                                                                       gTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATcc                  >  2:5309330/1‑93 (MQ=255)
                                                                        tACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCg                 >  2:4056251/1‑93 (MQ=255)
                                                                        tACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCg                 <  2:6266833/93‑1 (MQ=255)
                                                                          cTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGt               >  1:5411533/1‑93 (MQ=255)
                                                                           tGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTc              <  1:3334510/93‑1 (MQ=255)
                                                                            gACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCt             <  1:6335862/93‑1 (MQ=255)
                                                                               gggCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTAt          >  2:264160/1‑93 (MQ=255)
                                                                               gggCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTAt          >  2:5816952/1‑93 (MQ=255)
                                                                               gggCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTAt          <  1:581327/93‑1 (MQ=255)
                                                                                ggCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTAtt         <  2:4058893/93‑1 (MQ=255)
                                                                                 gCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTg        <  1:1570996/93‑1 (MQ=255)
                                                                                 gCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTg        >  2:1885379/1‑93 (MQ=255)
                                                                                  cTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGt       >  2:5926745/1‑93 (MQ=255)
                                                                                    gtgAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCa     >  2:5145834/1‑93 (MQ=255)
                                                                                    gtgAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCa     <  1:5981317/93‑1 (MQ=255)
                                                                                    gtgAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCa     >  2:4100228/1‑93 (MQ=255)
                                                                                    gtgAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCa     >  1:1531645/1‑93 (MQ=255)
                                                                                    gtgAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCa     >  2:3687503/1‑93 (MQ=255)
                                                                                          |                                                                                         
CCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAA  >  NC_000913/331929‑332108

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