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,3410TC99.2% 1274.7 / ‑3.5 370intergenic (‑342/‑343)ygdP/mutHputative invasion protein/methyl‑directed mismatch repair
Reads supporting (aligned to +/- strand):  ref base T (1/1);  major base C (214/151);  minor base A (1/2);  total (216/154)
Fisher's exact test for biased strand distribution p-value = 5.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

Only 100 of 374 total aligned reads displayed.

AAAATATTCATTATCAGTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGTTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCACTAAAAACTGGAGTTTCGTCGCAC  >  NC_000913/2967251‑2967430
                                                                                          |                                                                                         
aaaaTATTCATTATCAGTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTg                                                                                         <  1:768896/93‑1 (MQ=255)
 aaaTATTCATTATCAGTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTgg                                                                                        >  1:2919330/1‑93 (MQ=255)
   atatTCATTATCAGTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCg                                                                                      <  1:3025785/93‑1 (MQ=255)
   atatTCATTATCAGTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCg                                                                                      >  1:1720504/1‑93 (MQ=255)
   atatTCATTATCAGTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCg                                                                                      <  1:8055254/93‑1 (MQ=255)
   atatTCATTATCAGTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCg                                                                                      <  2:3796730/93‑1 (MQ=255)
    tatTCATTATCAGTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGa                                                                                     <  1:1261868/93‑1 (MQ=255)
    tatTCATTATCAGTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGa                                                                                     <  2:6197282/93‑1 (MQ=255)
       tCATTATCAGTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAg                                                                                  >  2:7609146/1‑93 (MQ=255)
        cATTATCAGTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGa                                                                                 <  1:613788/93‑1 (MQ=255)
        cATTATCAGTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGa                                                                                 <  1:5663200/93‑1 (MQ=255)
         aTTATCAGTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGAt                                                                                <  2:5537145/93‑1 (MQ=255)
         aTTATCAGTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGAt                                                                                >  2:6810445/1‑93 (MQ=255)
          ttATCAGTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATc                                                                               >  2:6039590/1‑93 (MQ=255)
           tATCAGTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCa                                                                              >  2:1883192/1‑93 (MQ=255)
            aTCAGTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCaa                                                                             >  1:8563544/1‑93 (MQ=255)
            aTCAGTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCaa                                                                             >  1:664072/1‑93 (MQ=255)
             tCAGTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAAc                                                                            <  2:2017747/93‑1 (MQ=255)
             tCAGTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAAc                                                                            >  1:3623880/1‑93 (MQ=255)
               aGTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCa                                                                          >  2:3401752/1‑93 (MQ=255)
               aGTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCa                                                                          >  1:3415483/1‑93 (MQ=255)
                gTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCaa                                                                         >  1:920517/1‑93 (MQ=255)
                 tCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAAt                                                                        <  2:8726342/93‑1 (MQ=255)
                 tCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAAt                                                                        <  2:1260057/93‑1 (MQ=255)
                  cAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATg                                                                       <  1:1625239/93‑1 (MQ=255)
                   aaTGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGc                                                                      >  2:4296627/1‑93 (MQ=255)
                   aaTGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGc                                                                      >  2:4537584/1‑93 (MQ=255)
                   aaTGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGc                                                                      <  1:5860909/93‑1 (MQ=255)
                    aTGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGcc                                                                     >  2:830231/1‑93 (MQ=255)
                    aTGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGcc                                                                     <  1:1062915/93‑1 (MQ=255)
                    aTGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGcc                                                                     >  1:2131237/1‑93 (MQ=255)
                     tGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCg                                                                    >  1:6915842/1‑93 (MQ=255)
                      gCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGt                                                                   <  1:2848566/93‑1 (MQ=255)
                        aTTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTAt                                                                 <  2:1518394/93‑1 (MQ=255)
                         ttACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTAtt                                                                <  1:8655072/93‑1 (MQ=255)
                         ttACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTAtt                                                                >  1:7159825/1‑93 (MQ=255)
                          tACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTAttt                                                               >  1:625362/1‑93 (MQ=255)
                            cATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTAt                                                             <  2:5065788/93‑1 (MQ=255)
                            cATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTAt                                                             >  2:1052611/1‑93 (MQ=255)
                             aTGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATc                                                            >  1:5321130/1‑93 (MQ=255)
                             aTGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATc                                                            >  2:4918422/1‑93 (MQ=255)
                             aTGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATc                                                            >  2:6644097/1‑93 (MQ=255)
                             aTGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATc                                                            <  2:4061311/93‑1 (MQ=255)
                              tGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATcc                                                           >  2:2193156/1‑93 (MQ=255)
                              tGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATcc                                                           >  1:539529/1‑93 (MQ=255)
                              tGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATcc                                                           >  1:2789866/1‑93 (MQ=255)
                                tAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCcac                                                         >  1:5493992/1‑93 (MQ=255)
                                 aaTTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCcaca                                                        <  2:8013898/93‑1 (MQ=255)
                                     ttAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAAt                                                    <  2:1659408/93‑1 (MQ=255)
                                     ttAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAAt                                                    <  2:7098774/93‑1 (MQ=255)
                                      tAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAAtg                                                   <  1:76329/93‑1 (MQ=255)
                                      tAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAAtg                                                   >  2:1401607/1‑93 (MQ=255)
                                            tGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCAc                                             >  2:3070233/1‑93 (MQ=255)
                                            tGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCAc                                             >  1:4426415/1‑93 (MQ=255)
                                            tGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCAc                                             >  2:7410745/1‑93 (MQ=255)
                                              aaaaaGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTa                                           <  2:8404871/93‑1 (MQ=255)
                                              aaaaaGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTa                                           <  1:1990030/93‑1 (MQ=255)
                                              aaaaaGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTa                                           >  2:2907791/1‑93 (MQ=255)
                                               aaaaGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTaa                                          >  1:5206814/1‑93 (MQ=255)
                                                aaaGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAg                                         >  1:4113480/1‑93 (MQ=255)
                                                 aaGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGt                                        <  1:2682286/93‑1 (MQ=255)
                                                 aaGAGTAATTCGTGACCCAGGTCACACCTCTCATGTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGt                                        <  2:3556295/93‑1 (MQ=255)
                                                  agagTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCAATAAGtt                                       >  2:2773529/1‑93 (MQ=255)
                                                   gagTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTa                                      >  2:155230/1‑93 (MQ=255)
                                                   gagTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTa                                      >  2:5646927/1‑93 (MQ=255)
                                                   gagTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTa                                      <  1:4613534/93‑1 (MQ=255)
                                                    agTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTaa                                     >  1:7383777/1‑93 (MQ=255)
                                                     gTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAg                                    >  2:8204342/1‑93 (MQ=255)
                                                      tAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGc                                   >  2:7515354/1‑93 (MQ=255)
                                                      tAATTCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGc                                   <  2:2306360/93‑1 (MQ=255)
                                                         ttCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACt                                <  1:4984757/93‑1 (MQ=255)
                                                         ttCGTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAAACAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACt                                >  1:4731171/1‑93 (MQ=255)
                                                            gTGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGaa                             >  1:3840974/1‑93 (MQ=255)
                                                             tGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAAc                            <  2:5356026/93‑1 (MQ=255)
                                                             tGACCCAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAAc                            >  2:2709757/1‑93 (MQ=255)
                                                                 ccAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCACt                        >  1:6496772/1‑93 (MQ=255)
                                                                  cAGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCACTa                       >  1:3306865/1‑93 (MQ=255)
                                                                   aGGTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCACTaa                      >  2:6890627/1‑93 (MQ=255)
                                                                    ggTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCACTaaa                     >  1:6381934/1‑93 (MQ=255)
                                                                    ggTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCACTaaa                     >  1:5933264/1‑93 (MQ=255)
                                                                    ggTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCACTaaa                     >  1:5080432/1‑93 (MQ=255)
                                                                     gTCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCACTaaaa                    >  1:8326794/1‑93 (MQ=255)
                                                                      tCACACCTTTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCACTaaaaa                   >  1:3756099/1‑93 (MQ=255)
                                                                      tCACACCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCACTaaaaa                   >  2:360727/1‑93 (MQ=255)
                                                                       cacaCCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCACTAAAAAc                  >  2:6221645/1‑93 (MQ=255)
                                                                        acacCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCACTAAAAACt                 >  2:1134188/1‑93 (MQ=255)
                                                                         cacCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCACTAAAAACTg                >  1:7880874/1‑93 (MQ=255)
                                                                         cacCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCACTAAAAACTg                <  1:5568620/93‑1 (MQ=255)
                                                                          acCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCACTAAAAACTgg               >  2:6475606/1‑93 (MQ=255)
                                                                          acCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCACTAAAAACTgg               >  1:5103982/1‑93 (MQ=255)
                                                                          acCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCACTAAAAACTgg               <  1:2480087/93‑1 (MQ=255)
                                                                          acCTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCACTAAAAACTgg               >  1:7233042/1‑93 (MQ=255)
                                                                           ccTCTCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCACTAAAAACTGGa              <  2:2488695/93‑1 (MQ=255)
                                                                            ctctCATTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCACTAAAAACTGGAg             <  1:3489047/93‑1 (MQ=255)
                                                                                 aTTTACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCACTAAAAACTGGAGTTtcg        >  2:8522463/1‑93 (MQ=255)
                                                                                  tttACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCACTAAAAACTGGAGTTtcgt       <  1:1574383/93‑1 (MQ=255)
                                                                                  tttACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCACTAAAAACTGGAGTTtcgt       <  1:3941680/93‑1 (MQ=255)
                                                                                  tttACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCACTAAAAACTGGAGTTtcgt       >  1:2263839/1‑93 (MQ=255)
                                                                                    tACGGGCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCACTAAAAACTGGAGTTtcgtcg     <  2:5877357/93‑1 (MQ=255)
                                                                                       gggCTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCACTAAAAACTGGAGTTTCGTCGCAc  >  2:4472041/1‑93 (MQ=255)
                                                                                          |                                                                                         
AAAATATTCATTATCAGTCAATGCATTACATGTAATTTTAGTAATGAAAAAGAGTAATTCGTGACCCAGGTCACACCTCTCATTTACGGGTTGGCGAAAGATCAACCAATGCCGTATTTATCCACAGAATGTGCCACTAAGTTAAGCACTGAACCACTAAAAACTGGAGTTTCGTCGCAC  >  NC_000913/2967251‑2967430

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