Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,967,341 T→C intergenic (‑342/‑343) ygdP ← / → mutH putative invasion protein/methyl‑directed mismatch repair

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,967,3410TC84.7% 666.2 / 103.9 287intergenic (‑342/‑343)ygdP/mutHputative invasion protein/methyl‑directed mismatch repair
Reads supporting (aligned to +/- strand):  ref base T (27/17);  new base C (140/103);  total (167/120)
Fisher's exact test for biased strand distribution p-value = 7.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 287 total aligned reads displayed.

AAAATATTCATTATCAGTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGTTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCACTAAAAACTGGAGTTTCGTCGCAC  >  NC_000913/2967251‑2967430
                                                                                          |                                                                                         
aaaaTATTCATTATCAGTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGTTg                                                                                         >  1:910573/1‑93 (MQ=255)
aaaaTATTCATTATCAGTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTg                                                                                         <  2:1210452/93‑1 (MQ=255)
aaaaTATTCATTATCAGTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTg                                                                                         <  2:6120325/93‑1 (MQ=255)
 aaaTATTCATTATCAGTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTgg                                                                                        >  1:5122257/1‑93 (MQ=255)
 aaaTATTCATTATCAGTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTgg                                                                                        <  2:6796305/93‑1 (MQ=255)
 aaaTATTCATTATCAGTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTgg                                                                                        >  2:5223811/1‑93 (MQ=255)
  aaTATTCATTATCAGTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGc                                                                                       <  1:2124499/93‑1 (MQ=255)
  aaTATTCATTATCAGTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGc                                                                                       >  1:3476475/1‑93 (MQ=255)
   atatTCATTATCAGTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGTTGGCg                                                                                      >  2:5992507/1‑93 (MQ=255)
      ttCATTATCAGTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGTTGGCGaaa                                                                                   <  1:6985693/93‑1 (MQ=255)
       tCATTATCAGTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAg                                                                                  <  1:5457336/93‑1 (MQ=255)
       tCATTATCAGTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAg                                                                                  >  2:6403606/1‑93 (MQ=255)
         aTTATCAGTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGAt                                                                                <  2:2893337/93‑1 (MQ=255)
          ttATCAGTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGTTGGCGAAAGATc                                                                               >  1:1933211/1‑93 (MQ=255)
          ttATCAGTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATc                                                                               <  2:4404578/93‑1 (MQ=255)
          ttATCAGTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATc                                                                               >  2:4026425/1‑93 (MQ=255)
          ttATCAGTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATc                                                                               <  2:828681/93‑1 (MQ=255)
          ttATCAGTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATc                                                                               >  1:2732443/1‑93 (MQ=255)
            aTCAGTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGCCCCAGGTCACCCCTCTCATTTACGGGCTGGCGAAAGATCaa                                                                             <  1:3660439/93‑1 (MQ=255)
            aTCAGTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCaa                                                                             <  2:211643/93‑1 (MQ=255)
             tCAGTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAAc                                                                            <  2:5085852/93‑1 (MQ=255)
              cAGTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGTTGGCGAAAGATCAAcc                                                                           <  1:5075218/93‑1 (MQ=255)
              cAGTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAAcc                                                                           <  1:1042187/93‑1 (MQ=255)
                gTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGTTGGCGAAAGATCAACCaa                                                                         >  1:2802396/1‑93 (MQ=255)
                gTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCaa                                                                         >  2:7182166/1‑93 (MQ=255)
                gTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCaa                                                                         <  1:3720775/93‑1 (MQ=255)
                gTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCaa                                                                         >  2:652619/1‑93 (MQ=255)
                 tCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAAt                                                                        >  2:7314332/1‑93 (MQ=255)
                  cAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATg                                                                       <  1:5680051/93‑1 (MQ=255)
                   aaTGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGTTGGCGAAAGATCAACCAATGc                                                                      >  2:1705965/1‑93 (MQ=255)
                   aaTGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGTTGGCGAAAGATCAACCAATGc                                                                      >  2:2665008/1‑93 (MQ=255)
                   aaTGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGc                                                                      <  2:2400928/93‑1 (MQ=255)
                     tGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCg                                                                    >  2:5583433/1‑93 (MQ=255)
                     tGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCg                                                                    >  1:2913137/1‑93 (MQ=255)
                        aTTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTAt                                                                 >  1:42120/1‑93 (MQ=255)
                        aTTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTAt                                                                 >  1:456213/1‑93 (MQ=255)
                         ttACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGTTGGCGAAAGATCAACCAATGCCGTAtt                                                                <  1:6815490/93‑1 (MQ=255)
                         ttACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTAtt                                                                <  2:1566935/93‑1 (MQ=255)
                          tACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTAttt                                                               <  1:440038/93‑1 (MQ=255)
                          tACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTAttt                                                               >  1:1604400/1‑93 (MQ=255)
                           aCATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGTTGGCGAAAGATCAACCAATGCCGTATTTa                                                              <  1:7440656/93‑1 (MQ=255)
                              tGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCGATGCCGTATTTATcc                                                           <  1:7121208/93‑1 (MQ=255)
                              tGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATcc                                                           >  2:6866777/1‑93 (MQ=255)
                              tGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATcc                                                           >  2:4253640/1‑93 (MQ=255)
                                tAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCcac                                                         >  2:925561/1‑93 (MQ=255)
                                 aaTTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGTTGGCGAAAGATCAACCAATGCCGTATTTATCcaca                                                        <  2:2548941/93‑1 (MQ=255)
                                   ttttAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGa                                                      >  2:3428045/1‑93 (MQ=255)
                                    ttttGTAATGAAAAAGAGTAATTCGTGACCCAGGGCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGGATTTTTCCACAGaa                                                     >  1:3003860/1‑93 (MQ=255)
                                    tttAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGaa                                                     <  1:410419/93‑1 (MQ=255)
                                     ttAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGTTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAAt                                                    >  1:1066174/1‑93 (MQ=255)
                                     ttAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAAt                                                    >  1:6598337/1‑93 (MQ=255)
                                     ttAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAAt                                                    <  1:6340696/93‑1 (MQ=255)
                                      tAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAAtg                                                   <  2:1049688/93‑1 (MQ=255)
                                       aGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGTTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAAtgt                                                  <  2:5336868/93‑1 (MQ=255)
                                       aGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAAtgt                                                  <  2:4664319/93‑1 (MQ=255)
                                         tAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGc                                                >  2:3174750/1‑93 (MQ=255)
                                          aaTGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGTTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGcc                                               >  1:5446475/1‑93 (MQ=255)
                                          aaTGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGcc                                               >  1:2321240/1‑93 (MQ=255)
                                            tGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCAc                                             >  1:3968155/1‑93 (MQ=255)
                                            tGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCAc                                             >  1:5880210/1‑93 (MQ=255)
                                             gAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCTGTATTTATCCACAGAATGTGCCANt                                            >  2:7007325/1‑93 (MQ=255)
                                               aaaaGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTaa                                          >  1:3408410/1‑93 (MQ=255)
                                                aaaGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAg                                         >  2:453847/1‑93 (MQ=255)
                                                 aaGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGTTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGt                                        >  2:5897466/1‑93 (MQ=255)
                                                 aaGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGt                                        >  1:2987964/1‑93 (MQ=255)
                                                 aaGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGt                                        <  2:1971899/93‑1 (MQ=255)
                                                   gagTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTa                                      <  1:1696402/93‑1 (MQ=255)
                                                   gagTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTa                                      <  1:3582408/93‑1 (MQ=255)
                                                    agTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTaa                                     <  1:2600426/93‑1 (MQ=255)
                                                    agTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTaa                                     <  1:3298634/93‑1 (MQ=255)
                                                     gTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGTTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAg                                    >  2:7427191/1‑93 (MQ=255)
                                                      tAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGc                                   <  2:232032/93‑1 (MQ=255)
                                                      tAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGc                                   <  1:4924945/93‑1 (MQ=255)
                                                          tCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTg                               <  1:298341/93‑1 (MQ=255)
                                                           cGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGa                              <  1:1294757/93‑1 (MQ=255)
                                                             tGACCCAGGTCACACCTCTCATTTACGGGTTGGCGAGAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAAc                            <  2:2503174/93‑1 (MQ=255)
                                                             tGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAAc                            >  2:3247258/1‑93 (MQ=255)
                                                              gACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAAcc                           <  1:7531490/93‑1 (MQ=255)
                                                               aCCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCa                          >  2:3614333/1‑93 (MQ=255)
                                                               aCCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCa                          <  2:2936022/93‑1 (MQ=255)
                                                                cccAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCAc                         >  1:4781119/1‑93 (MQ=255)
                                                                cccAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCAc                         >  1:1428212/1‑93 (MQ=255)
                                                                   aGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCACTaa                      >  2:2795719/1‑93 (MQ=255)
                                                                   aGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCACTaa                      >  1:6053120/1‑93 (MQ=255)
                                                                     gTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCACTaaaa                    >  1:5268472/1‑93 (MQ=255)
                                                                      tCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCACTaaaaa                   <  2:5377331/93‑1 (MQ=255)
                                                                       cactCCTCTCATTTACGGGTTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCACTAAAAAc                  <  2:2856290/93‑1 (MQ=255)
                                                                        acacCTCTCATTTACGGGTTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCACTAAAAACt                 <  2:6209995/93‑1 (MQ=255)
                                                                        acacCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCACTAAAAACt                 >  2:3914834/1‑93 (MQ=255)
                                                                         cacCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCACTAAAAACTg                >  2:3593650/1‑93 (MQ=255)
                                                                         cacCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCACTAAAAACTg                <  1:3223584/93‑1 (MQ=255)
                                                                          acCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCACTAAAAACTga               >  2:1650435/1‑92 (MQ=255)
                                                                          acCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATATATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCACTAAAAACTgg               >  2:5687908/1‑93 (MQ=255)
                                                                             tctcATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCACTAAAAACTGGAGt            <  2:5009792/93‑1 (MQ=255)
                                                                              ctcATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCACTAAAAACTGGAGtt           <  2:4388586/93‑1 (MQ=255)
                                                                               tcATTTACGGGTTGGCGAGAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCACTAAAAACTGGAGttt          <  2:6703718/93‑1 (MQ=255)
                                                                                cATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCACTAAAAACTGGAGTTtc         <  1:1823241/93‑1 (MQ=255)
                                                                                cATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCACTAAAAACTGGAGTTtc         <  1:2460612/93‑1 (MQ=255)
                                                                                 aTTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCACTAAAAACTGGAGTTtcg        >  1:6684958/1‑93 (MQ=255)
                                                                                     aCGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCACTAAAAACTGGAGTTTCGTCGc    <  1:7416315/93‑1 (MQ=255)
                                                                                          |                                                                                         
AAAATATTCATTATCAGTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGTTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCACTAAAAACTGGAGTTTCGTCGCAC  >  NC_000913/2967251‑2967430

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