Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,568,135 G→A P127S (CCC→TCC)  xasA ← acid sensitivity protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,568,1350GA89.5% 685.3 / 41.0 220P127S (CCC→TCC) xasAacid sensitivity protein
Reads supporting (aligned to +/- strand):  ref base G (8/15);  new base A (86/111);  total (94/126)
Fisher's exact test for biased strand distribution p-value = 5.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.03e-01

Only 100 of 221 total aligned reads displayed.

TTTAGCAATTCGCGCCGTGTATTTCGTGCCACCAAACTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGGGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCATCGGAATAAAACCAATGGCGATTTGCA  >  NC_000913/1568043‑1568224
                                                                                            |                                                                                         
tttAGCAATTCGCGCCGTGTATTTCGTGCCACCAAACTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATggg                                                                                           >  2:7134653/1‑93 (MQ=255)
   aGCAATTCGCGCCGTGTATTTCGTGCCACCAAACTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTc                                                                                        <  1:1787111/93‑1 (MQ=255)
    nnaaTTCGCGCCGTGTATTTCGTGCCACCAAACTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCt                                                                                       <  1:2420742/91‑1 (MQ=255)
    gCAATTCGCGCCGTGTATTTCGTGCCACCAAACTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCt                                                                                       <  1:2236492/93‑1 (MQ=255)
     cAATTCGCGCCGTGTATTTCGTGCCACCAAACTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCtt                                                                                      <  2:2498328/93‑1 (MQ=255)
     cAATTCGCGCCGTGTATTTCGTGCCACCAAACTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCtt                                                                                      >  1:3869449/1‑93 (MQ=255)
      aaTTCGCGCCGTGTATTTCGTGCCACCAAACTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCttc                                                                                     <  1:4619369/93‑1 (MQ=255)
       aTTCGCGCCGTGTATTTCGTGCCACCAAACTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCttca                                                                                    <  2:4918689/93‑1 (MQ=255)
        ttCGCGCCGTGTATTTCGTGCCACCAAACTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCttcat                                                                                   <  1:4203588/93‑1 (MQ=255)
        ttCGCGCCGTGTATTTCGTGCCACCAAACTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCttcat                                                                                   >  1:3474392/1‑93 (MQ=255)
          cgcgCCGTGTATTTCGTGCCACCAAACTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGGGTCttcattc                                                                                 <  2:6725540/93‑1 (MQ=255)
           gcgcCGTGTATTTCGTGCCACCAAACTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCttcattca                                                                                <  2:1406931/93‑1 (MQ=255)
             gcCGTGTATTTCGTGCCACCAAACTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAgc                                                                              >  1:5646976/1‑93 (MQ=255)
             gcCGTGTATTTCGTGCCACCAAACTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAgc                                                                              >  1:5536190/1‑93 (MQ=255)
               cGTGTATTTCGTGCCACCAAACTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAgcgc                                                                            <  1:5222405/93‑1 (MQ=255)
               cGTGTATTTCGTGCCACCAAACTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAgcgc                                                                            >  1:4100364/1‑93 (MQ=255)
                 tgtATTTCGTGCCACCAAACTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTg                                                                          <  1:5266789/93‑1 (MQ=255)
                 tgtATTTCGTGCCACCAAACTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTg                                                                          >  1:1603211/1‑93 (MQ=255)
                  gtATTTCGTGCCACCAAACTGCGTTAATGCCAGCGCTCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTgg                                                                         <  2:6487530/93‑1 (MQ=255)
                     tttCGTGCCACCAAACTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCa                                                                      >  2:75062/1‑93 (MQ=255)
                      ttCGTGCCACCAAACTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCAt                                                                     >  1:3962010/1‑93 (MQ=255)
                         gTGCCACCAAACTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTc                                                                  <  2:5929415/93‑1 (MQ=255)
                          tGCCACCAAACTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCa                                                                 >  2:5800380/1‑93 (MQ=255)
                            ccaccaAACTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAgg                                                               >  2:3735674/1‑93 (MQ=255)
                            ccaccaAACTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAgg                                                               >  1:5001611/1‑93 (MQ=255)
                            ccaccaAACTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAgg                                                               >  1:5577823/1‑93 (MQ=255)
                             caccaAACTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGGGTCTTCATTCAGCGCTGGCCATTTCAGGa                                                              >  1:6702251/1‑93 (MQ=255)
                              accaAACTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGAt                                                             >  2:4423463/1‑93 (MQ=255)
                              accaAACTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGAt                                                             <  2:6143495/93‑1 (MQ=255)
                               ccaAACTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATg                                                            >  2:1885590/1‑93 (MQ=255)
                                caAACTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGt                                                           >  2:557692/1‑93 (MQ=255)
                                caAACTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGt                                                           >  2:3411405/1‑93 (MQ=255)
                                 aaaCTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTa                                                          >  1:1299949/1‑93 (MQ=255)
                                 aaaCTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTa                                                          <  2:2202981/93‑1 (MQ=255)
                                 aaaCTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTa                                                          <  1:4313517/93‑1 (MQ=255)
                                  aaCTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGGGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAg                                                         <  2:6229509/93‑1 (MQ=255)
                                  aaCTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAg                                                         <  1:7234364/93‑1 (MQ=255)
                                    cTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGga                                                       <  2:5020383/93‑1 (MQ=255)
                                    cTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGga                                                       <  1:6158474/93‑1 (MQ=255)
                                     tGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGgag                                                      >  2:4096757/1‑93 (MQ=255)
                                      gCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGgaga                                                     <  2:5561324/93‑1 (MQ=255)
                                      gCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGgaga                                                     <  1:7078052/93‑1 (MQ=255)
                                      gCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGgaga                                                     <  1:333303/93‑1 (MQ=255)
                                       cGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGagag                                                    <  1:4000031/93‑1 (MQ=255)
                                       cGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGagag                                                    <  1:692200/93‑1 (MQ=255)
                                         ttAATGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTg                                                  >  1:3914311/1‑93 (MQ=255)
                                         ttAATGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTg                                                  >  2:3542580/1‑93 (MQ=255)
                                         ttAATGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTg                                                  <  2:298205/93‑1 (MQ=255)
                                           aaTGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGcc                                                >  1:4483031/1‑93 (MQ=255)
                                             tGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGcccc                                              <  2:6342447/93‑1 (MQ=255)
                                              gCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCt                                             <  2:4002336/93‑1 (MQ=255)
                                              gCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCt                                             >  1:6421166/1‑93 (MQ=255)
                                              gCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCt                                             <  1:2040884/93‑1 (MQ=255)
                                                cAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTaa                                           <  2:2360813/93‑1 (MQ=255)
                                                 aGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAac                                          >  1:2766479/1‑93 (MQ=255)
                                                  gcgcCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGGGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAaca                                         <  1:4283258/93‑1 (MQ=255)
                                                  gcgcCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAaca                                         >  1:3582947/1‑93 (MQ=255)
                                                   cgcCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAacac                                        <  1:3001217/93‑1 (MQ=255)
                                                   cgcCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAacac                                        <  1:5084405/93‑1 (MQ=255)
                                                     cccAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGa                                      >  1:3177934/1‑93 (MQ=255)
                                                     cccAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGa                                      <  1:7020688/93‑1 (MQ=255)
                                                       cAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGaaa                                    >  1:1696464/1‑93 (MQ=255)
                                                        aaaGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAAt                                   <  2:5030652/93‑1 (MQ=255)
                                                        aaaGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAAt                                   >  2:6900149/1‑93 (MQ=255)
                                                         aaGAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATa                                  <  1:5856617/93‑1 (MQ=255)
                                                           gAATGATGAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATaga                                <  2:2639273/93‑1 (MQ=255)
                                                             atgatgAGTGCTGCAATAGTTTTGGTAATGGGGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGc                              >  2:6282859/1‑93 (MQ=255)
                                                             atgatgAGTGCTGCAATAGTTTTGGTAATGGGGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGc                              <  2:657189/93‑1 (MQ=255)
                                                             atgatgAGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGc                              <  2:5404448/93‑1 (MQ=255)
                                                               gatgaGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCAt                            <  1:4632224/93‑1 (MQ=255)
                                                               gatgaGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCAt                            <  2:1092179/93‑1 (MQ=255)
                                                               gatgaGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCAt                            <  1:2706920/93‑1 (MQ=255)
                                                               gatgaGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCAt                            <  2:4500452/93‑1 (MQ=255)
                                                                atgaGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCATc                           <  2:7529282/93‑1 (MQ=255)
                                                                atgaGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCATc                           <  2:2559572/93‑1 (MQ=255)
                                                                 tgaGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCATCg                          <  2:3202065/93‑1 (MQ=255)
                                                                  gaGTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTACCACGAAATAGAGCATCgg                         <  1:2830784/93‑1 (MQ=255)
                                                                    gTGCTGCAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCATCGGaa                       <  1:7209549/93‑1 (MQ=255)
                                                                     tgctgcAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCATCGGAAt                      <  2:3669430/93‑1 (MQ=255)
                                                                      gctgcAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCATCGGAATa                     >  2:179035/1‑93 (MQ=255)
                                                                       ctgcAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCATCGGAATaa                    >  2:6391475/1‑93 (MQ=255)
                                                                        tgcAATAGTTTTGGTAATGGGGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCATCGGAATaaa                   >  2:749743/1‑93 (MQ=255)
                                                                        tgcAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCATCGGAATaaa                   >  1:4382559/1‑93 (MQ=255)
                                                                          cAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCATCGGAATAAAAc                 >  1:4250534/1‑93 (MQ=255)
                                                                          cAATAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCATCGGAATAAAAc                 <  1:5447985/93‑1 (MQ=255)
                                                                            aTAGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCATCGGAATAAAACCa               <  1:823440/93‑1 (MQ=255)
                                                                              aGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGGGTGCCCCTAACACGAAATAGAGCATCGGAATAAAACCAAt             <  2:4878963/93‑1 (MQ=255)
                                                                              aGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCATCGGAATAAAACCAAt             <  1:3108727/93‑1 (MQ=255)
                                                                              aGTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCATCGGAATAAAACCAAt             <  1:6373043/93‑1 (MQ=255)
                                                                               gTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCATCGGAATAAAACCAATg            <  2:6090401/93‑1 (MQ=255)
                                                                               gTTTTGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCATCGGAATAAAACCAATg            <  1:6628586/93‑1 (MQ=255)
                                                                                ttttGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCATCGGAATAAAACCAATgg           >  1:284528/1‑93 (MQ=255)
                                                                                ttttGGTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCATCGGAATAAAACCAATgg           <  2:2700206/93‑1 (MQ=255)
                                                                                    ggTAATGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCATCGGAATAAAACCAATGGCGAt       <  1:2406308/93‑1 (MQ=255)
                                                                                        aTGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCATCGGAATAAAACCAATGGCGATTTGc   >  1:795488/1‑93 (MQ=255)
                                                                                        aTGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCATCGGAATAAAACCAATGGCGATTTGc   <  1:4679388/93‑1 (MQ=255)
                                                                                        aTGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCATCGGAATAAAACCAATGGCGATTTGc   <  1:2305190/93‑1 (MQ=255)
                                                                                         tGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCATCGGAATAAAACCAATGGCGATTTGCa  <  2:1246786/93‑1 (MQ=255)
                                                                                         tGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCATCGGAATAAAACCAATGGCGATTTGCa  <  1:102877/93‑1 (MQ=255)
                                                                                         tGGAGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCATCGGAATAAAACCAATGGCGATTTGCa  >  2:4646148/1‑93 (MQ=255)
                                                                                            |                                                                                         
TTTAGCAATTCGCGCCGTGTATTTCGTGCCACCAAACTGCGTTAATGCCAGCGCCCAAAGAATGATGAGTGCTGCAATAGTTTTGGTAATGGGGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCATCGGAATAAAACCAATGGCGATTTGCA  >  NC_000913/1568043‑1568224

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