Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,167,940 G→A G33S (GGT→AGT)  qseB → quorum sensing regulator B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,167,9400GA99.6% 969.0 / ‑5.4 248G33S (GGT→AGT) qseBquorum sensing regulator B
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (130/117);  total (131/117)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

Only 100 of 250 total aligned reads displayed.

ATTTTACTGATAGAAGATGACATGCTGATTGGCGACGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAGGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCATGGATGGTCGCGATATT  >  NC_000913/3167850‑3168029
                                                                                          |                                                                                         
aTTTTACTGATAGAAGATGACATGCTGATTGGCGACGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAgt                                                                                         <  1:762294/93‑1 (MQ=255)
 ttttACTGATAGAAGATGACATGCTGATTGGCGACGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAgtc                                                                                        <  1:1173161/93‑1 (MQ=255)
 ttttACTGATAGAAGATGACATGCTGATTGGCGACGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAgtc                                                                                        <  2:4625152/93‑1 (MQ=255)
  tttACTGATAGAAGATGACATGCTGATTGGCGACGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAgtcg                                                                                       >  2:784366/1‑93 (MQ=255)
   ttACTGATAGAAGATGACATGCTGATTGGCGACGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAgtcgt                                                                                      >  1:5922859/1‑93 (MQ=255)
   ttACTGATAGAAGATGACATGCTGATTGGCGACGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAgtcgt                                                                                      >  2:233289/1‑93 (MQ=255)
   ttACTGATAGAAGATGACATGCTGATTGGCGACGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAgtcgt                                                                                      <  2:1328733/93‑1 (MQ=255)
    tACTGATAGAAGATGACATGCTGATTGGCGACGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAgtcgtc                                                                                     <  2:2776976/93‑1 (MQ=255)
    tACTGATAGAAGATGACATGCTGATTGGCGACGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAgtcgtc                                                                                     <  1:3986463/93‑1 (MQ=255)
     aCTGATAGAAGATGACATGCTGATTGGCGACGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCa                                                                                    <  1:1036583/93‑1 (MQ=255)
     aCTGATAGAAGATGACATGCTGATTGGCGACGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCa                                                                                    >  1:3767842/1‑93 (MQ=255)
     aCTGATAGAAGATGACATGCTGATTGGCGACGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCa                                                                                    >  2:3946299/1‑93 (MQ=255)
     aCTGATAGAAGATGACATGCTGATTGGCGACGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCa                                                                                    >  1:3329669/1‑93 (MQ=255)
      cTGATAGAAGATGACATGCTGATTGGCGACGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAg                                                                                   >  2:1164292/1‑93 (MQ=255)
       tGATAGAAGATGACATGCTGATTGGCGACGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACAAAAAGTCGTCAgg                                                                                  <  1:5098309/93‑1 (MQ=255)
        gATAGAAGATGACATGCTGATTGGCGACGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAggg                                                                                 >  2:3425144/1‑93 (MQ=255)
        gATAGAAGATGACATGCTGATTGGCGACGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAggg                                                                                 <  2:3910355/93‑1 (MQ=255)
         aTAGAAGATGACATGCTGATTGGCGACGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGa                                                                                <  2:3882326/93‑1 (MQ=255)
         aTAGAAGATGACATGCTGATTGGCGACGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGa                                                                                >  1:3221310/1‑93 (MQ=255)
         aTAGAAGATGACATGCTGATTGGCGACGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGa                                                                                <  1:3565345/93‑1 (MQ=255)
         aTAGAAGATGACAGGCTGATTGGCGACGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGTTTTACACAAAGTCGTCAGGGa                                                                                <  2:1636413/93‑1 (MQ=255)
          tAGAAGATGACATGCTGATTGGCGACGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGaa                                                                               >  1:1757101/1‑93 (MQ=255)
           agaagaTGACATGCTGATTGGCGACGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGaaa                                                                              >  1:5451241/1‑93 (MQ=255)
            gaagaTGACATGCTGATTGGCGACGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGaaaa                                                                             <  2:3964312/93‑1 (MQ=255)
            gaagaTGACATGCTGATTGGCGACGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGaaaa                                                                             <  1:5302156/93‑1 (MQ=255)
               gaTGACATGCTGATTGGCGACGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAagag                                                                          <  2:2685173/93‑1 (MQ=255)
               gaTGACATGCTGATTGGCGACGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAagag                                                                          >  2:5097537/1‑93 (MQ=255)
                 tGACATGCTGATTGGCGACGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGgc                                                                        >  2:299221/1‑93 (MQ=255)
                 tGACATGCTGATTGGCGACGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGgc                                                                        <  2:6128510/93‑1 (MQ=255)
                  gACATGCTGATTGGCGACGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGgcg                                                                       >  2:2837855/1‑93 (MQ=255)
                   aCATGCTGATTGGCGACGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGgcgc                                                                      >  1:21136/1‑93 (MQ=255)
                    cATGCTGATTGGCGACGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCt                                                                     >  1:4509534/1‑93 (MQ=255)
                     aTGCTGATTGGCGACGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAGGTCGTCAGGGAAAAGAGGCGCtt                                                                    >  2:3536146/1‑93 (MQ=255)
                     aTGCTGATTGGCGACGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCtt                                                                    <  1:1996449/93‑1 (MQ=255)
                      tGCTGATTGGCGACGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCttt                                                                   <  2:4475599/93‑1 (MQ=255)
                      tGCTGATTGGCGACGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCttt                                                                   >  1:3839627/1‑93 (MQ=255)
                       gCTGATTGGCGACGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTta                                                                  >  2:6093368/1‑93 (MQ=255)
                        cTGATTGGCGACGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTtat                                                                 <  1:4733164/93‑1 (MQ=255)
                          gATTGGCGACGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAg                                                               >  1:666602/1‑93 (MQ=255)
                            atGGCGACGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCATTATAgcg                                                             <  1:6288958/92‑1 (MQ=255)
                              ggCGACGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCa                                                           >  1:4680695/1‑93 (MQ=255)
                                cGACGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCAcc                                                         <  2:6308026/93‑1 (MQ=255)
                                 gACGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGATGCGCTTTATAGCGCACCt                                                        >  1:5557042/1‑93 (MQ=255)
                                 gACGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCt                                                        <  2:5054181/93‑1 (MQ=255)
                                 gACGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCt                                                        >  2:3367199/1‑93 (MQ=255)
                                   cGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTa                                                      <  1:4359729/93‑1 (MQ=255)
                                   cGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTa                                                      <  2:5398323/93‑1 (MQ=255)
                                   cGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTa                                                      <  2:5006773/93‑1 (MQ=255)
                                    gacATCAAAACGGGCCTTAGTAACATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTat                                                     <  2:3138330/91‑1 (MQ=255)
                                     gCATCAAAACGAGCCTTAGTAAAATGTGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTatg                                                    <  2:3450693/93‑1 (MQ=255)
                                      cATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCTACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTatga                                                   >  2:402182/1‑93 (MQ=255)
                                      cATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTatga                                                   >  2:3525628/1‑93 (MQ=255)
                                        tCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTatgatg                                                 >  1:2727847/1‑93 (MQ=255)
                                          aaaaCGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCg                                               <  1:2367677/93‑1 (MQ=255)
                                           aaaCGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCgg                                              >  2:4515869/1‑93 (MQ=255)
                                             aCGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTg                                            <  1:4455727/93‑1 (MQ=255)
                                             aCGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTg                                            >  2:1218647/1‑93 (MQ=255)
                                             aCGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTg                                            >  2:749747/1‑93 (MQ=255)
                                               gggCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGAt                                          >  2:5758063/1‑93 (MQ=255)
                                                ggCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATc                                         >  2:3597820/1‑93 (MQ=255)
                                                 gCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATcc                                        <  1:3732133/93‑1 (MQ=255)
                                                 gCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATcc                                        <  2:4765357/93‑1 (MQ=255)
                                                  ccTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCt                                       <  2:868401/93‑1 (MQ=255)
                                                  ccTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCt                                       <  1:2249245/93‑1 (MQ=255)
                                                    ttAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTgg                                     <  2:2210916/93‑1 (MQ=255)
                                                    ttAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTgg                                     >  2:2527758/1‑93 (MQ=255)
                                                    ttAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTgg                                     >  1:860180/1‑93 (MQ=255)
                                                       gTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGAtt                                  <  1:2165231/93‑1 (MQ=255)
                                                       gTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGAtt                                  >  1:3521491/1‑93 (MQ=255)
                                                       gTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGAtt                                  <  2:31518/93‑1 (MQ=255)
                                                        tAAAATGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGAttt                                 <  2:1727645/93‑1 (MQ=255)
                                                          aaaTGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTaa                               <  2:567469/93‑1 (MQ=255)
                                                          aaaTGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTaa                               >  1:5359244/1‑93 (MQ=255)
                                                            aTGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAAcc                             >  2:1817367/1‑93 (MQ=255)
                                                            aTGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAAcc                             >  1:3471197/1‑93 (MQ=255)
                                                             tGGGTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGTTTTAACCt                            >  2:5469184/1‑93 (MQ=255)
                                                             tGGCTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCCGGGTTTAACCt                            <  1:5617124/93‑1 (MQ=255)
                                                              gggTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCtt                           >  1:6054967/1‑93 (MQ=255)
                                                               ggTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTa                          <  1:1581122/93‑1 (MQ=255)
                                                               ggTTTTAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTa                          <  1:1946509/93‑1 (MQ=255)
                                                                 ttttAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTAcc                        >  1:6396972/1‑93 (MQ=255)
                                                                    tAGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAgg                     >  2:1468036/1‑93 (MQ=255)
                                                                     aGCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTAAAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGc                    <  1:363158/93‑1 (MQ=255)
                                                                      gCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCa                   >  2:3054278/1‑93 (MQ=255)
                                                                      gCGTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCa                   >  1:2473591/1‑93 (MQ=255)
                                                                        gTCGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCatg                 >  1:4969837/1‑93 (MQ=255)
                                                                          cGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCatgga               <  1:3146014/93‑1 (MQ=255)
                                                                          cGACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCatgga               <  1:3181875/93‑1 (MQ=255)
                                                                           gACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCatggat              >  2:1016380/1‑93 (MQ=255)
                                                                           gACTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCatggat              <  2:5194265/93‑1 (MQ=255)
                                                                            aCTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCatggatg             <  1:2855191/93‑1 (MQ=255)
                                                                             cTGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCatggatgg            >  1:5721794/1‑93 (MQ=255)
                                                                              tGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCatggatggc           >  1:5220520/1‑92 (MQ=255)
                                                                              tGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCATGGATGGt           <  2:3193683/93‑1 (MQ=255)
                                                                              tGGTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCATGGATGGt           <  1:1086820/93‑1 (MQ=255)
                                                                               ggTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCATGGATGGTc          <  2:5518750/93‑1 (MQ=255)
                                                                               ggTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCATGGATGGTc          >  1:2560544/1‑93 (MQ=255)
                                                                               ggTTTACACAAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCATGGATGGTc          <  1:1711587/93‑1 (MQ=255)
                                                                                     cacaAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCATGGATGGTCGCGata    >  2:4896307/1‑93 (MQ=255)
                                                                                       caAAGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCATGGATGGTCGCGATAtt  <  1:3163480/93‑1 (MQ=255)
                                                                                          |                                                                                         
ATTTTACTGATAGAAGATGACATGCTGATTGGCGACGGCATCAAAACGGGCCTTAGTAAAATGGGTTTTAGCGTCGACTGGTTTACACAAGGTCGTCAGGGAAAAGAGGCGCTTTATAGCGCACCTTATGATGCGGTGATCCTGGATTTAACCTTACCAGGCATGGATGGTCGCGATATT  >  NC_000913/3167850‑3168029

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