Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,764,932 A→G M593T (ATG→ACG)  ydiJ ← putative oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,764,9320AG83.1% 459.4 / 95.4 213M593T (ATG→ACG) ydiJputative oxidase
Reads supporting (aligned to +/- strand):  ref base A (15/21);  new base G (99/78);  total (114/99)
Fisher's exact test for biased strand distribution p-value = 1.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.54e-01

Only 100 of 217 total aligned reads displayed.

ATCCGCCAACAAACGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCGCAATGGGGATC  >  NC_000913/1764841‑1765021
                                                                                           |                                                                                         
aTCCGCCAACAAACGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACAt                                                                                          >  2:2719210/1‑93 (MQ=255)
 tCCGCCAACAAACGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACGTc                                                                                         >  1:4787232/1‑93 (MQ=255)
 tCCGCCAACAAACGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACGTc                                                                                         <  1:4646535/93‑1 (MQ=255)
    gCCAACAAACGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGa                                                                                      <  1:2238685/93‑1 (MQ=255)
     ccAACAAACGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGCCGGACACGTCGGAc                                                                                     <  2:6443789/93‑1 (MQ=255)
       aacaaACGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTa                                                                                   >  2:535594/1‑93 (MQ=255)
       aacaaACGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTa                                                                                   <  1:6122418/93‑1 (MQ=255)
       aacaaACGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTa                                                                                   >  1:2693166/1‑93 (MQ=255)
        acaaACGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTAc                                                                                  >  2:3837848/1‑93 (MQ=255)
         caaACGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTACg                                                                                 <  1:3599622/93‑1 (MQ=255)
          aaaCGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTACgg                                                                                <  1:1127410/93‑1 (MQ=255)
          aaaCGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTACgg                                                                                <  2:6810092/93‑1 (MQ=255)
          aaaCGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTACgg                                                                                <  2:5685997/93‑1 (MQ=255)
          aaaCGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTACgg                                                                                >  1:1986835/1‑93 (MQ=255)
           aaCGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTACggg                                                                               >  2:961921/1‑93 (MQ=255)
           aaCGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTACggg                                                                               <  1:928979/93‑1 (MQ=255)
            aCGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTACGGGc                                                                              <  1:5123459/93‑1 (MQ=255)
               cAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTACGGGCATc                                                                           <  2:5260315/93‑1 (MQ=255)
                aGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTACGGGCATCa                                                                          <  1:3755470/93‑1 (MQ=255)
                 gCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTACGGGCATCaa                                                                         >  2:243253/1‑93 (MQ=255)
                  ccATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTACGGGCATCaaa                                                                        <  2:7406810/93‑1 (MQ=255)
                    aTTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGt                                                                      >  2:1324464/1‑93 (MQ=255)
                     ttCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGtt                                                                     >  2:6605519/1‑93 (MQ=255)
                      tCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTg                                                                    <  1:7403992/93‑1 (MQ=255)
                        aCGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGAGCATCAAAGTTGaa                                                                  <  2:5491124/93‑1 (MQ=255)
                          gCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGc                                                                >  1:3285657/1‑93 (MQ=255)
                             ccAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATa                                                             >  2:448336/1‑93 (MQ=255)
                              cAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATaa                                                            >  1:3461698/1‑93 (MQ=255)
                               aGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATaaa                                                           >  1:2029832/1‑93 (MQ=255)
                                gCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAAc                                                          >  1:6316045/1‑93 (MQ=255)
                                gCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGAGCATCAAAGTTGAAGCATAAAc                                                          >  1:6543820/1‑93 (MQ=255)
                                gCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGGGATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAAc                                                          >  1:4236281/1‑93 (MQ=255)
                                   ttGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGt                                                       >  1:6249342/1‑93 (MQ=255)
                                   ttGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGt                                                       <  2:131256/93‑1 (MQ=255)
                                   ttGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGt                                                       <  1:289750/93‑1 (MQ=255)
                                    tGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGtt                                                      <  1:5051199/93‑1 (MQ=255)
                                    tGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGtt                                                      <  1:6074017/93‑1 (MQ=255)
                                    tGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGAGCATCAAAGTTGAAGCATAAACCGtt                                                      <  2:6519482/93‑1 (MQ=255)
                                      cgcgcCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGc                                                    <  2:4111702/93‑1 (MQ=255)
                                       gcgcCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGAGCATCAAAGTTGAAGCATAAACCGTTGcc                                                   >  2:339727/1‑93 (MQ=255)
                                        cgcCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCg                                                  <  2:1724706/93‑1 (MQ=255)
                                        cgcCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCg                                                  >  2:4022104/1‑93 (MQ=255)
                                        cgcCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCg                                                  >  2:4812483/1‑93 (MQ=255)
                                        cgcCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGAGCATCAAAGTTGAAGCATAAACCGTTGCCg                                                  <  1:1739166/93‑1 (MQ=255)
                                         gcCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGt                                                 >  1:2178696/1‑93 (MQ=255)
                                         gcCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGt                                                 <  1:4600696/93‑1 (MQ=255)
                                          cccTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGtt                                                <  1:1012594/93‑1 (MQ=255)
                                          cccTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGAGCATCAAAGTTGAAGCATAAACCGTTGCCGtt                                                <  2:7030419/93‑1 (MQ=255)
                                            cTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTac                                              >  2:11726/1‑93 (MQ=255)
                                             tttCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTaca                                             >  2:6881206/1‑93 (MQ=255)
                                               tCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACt                                           >  1:2661563/1‑93 (MQ=255)
                                                cGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTc                                          <  2:3029024/93‑1 (MQ=255)
                                                 ggTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTcc                                         <  2:6570866/93‑1 (MQ=255)
                                                   tgaatgaaTCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCAt                                       >  2:7305804/1‑93 (MQ=255)
                                                   tgaatgaaTCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCAt                                       >  2:2140705/1‑93 (MQ=255)
                                                   tgaatgaaTCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGAGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCAt                                       <  1:2556686/93‑1 (MQ=255)
                                                    gaatgaatCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATc                                      <  1:1399086/93‑1 (MQ=255)
                                                    gaatgaatCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGAGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATc                                      <  2:5830531/93‑1 (MQ=255)
                                                    gaatgaatCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGAGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATc                                      >  1:6010364/1‑93 (MQ=255)
                                                     aatgaatCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCg                                     >  1:7221495/1‑93 (MQ=255)
                                                        gaatCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCAc                                  >  1:7103625/1‑93 (MQ=255)
                                                        gaatCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCAc                                  <  1:5754667/93‑1 (MQ=255)
                                                        gaatCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCAc                                  <  1:826588/93‑1 (MQ=255)
                                                          atCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACcg                                >  2:6693401/1‑93 (MQ=255)
                                                           tCCGGTTCTGGGTGATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACcgc                               >  1:1603447/1‑93 (MQ=255)
                                                             cGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGAGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCgcgc                             <  2:1960380/93‑1 (MQ=255)
                                                               gTTCTGGGTGATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCa                           <  1:4710310/93‑1 (MQ=255)
                                                                ttCTGGGTGATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCAc                          <  1:1826118/93‑1 (MQ=255)
                                                                ttCTGGGTGATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCAc                          <  2:1083801/93‑1 (MQ=255)
                                                                ttCTGGGTGATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCAc                          <  2:1533845/93‑1 (MQ=255)
                                                                ttCTGGGTGATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCAc                          >  2:4397622/1‑93 (MQ=255)
                                                                   tGGGTGATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCActgc                       <  2:5161218/93‑1 (MQ=255)
                                                                   tGGGTGATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCActgc                       >  1:4137704/1‑93 (MQ=255)
                                                                   tGGGTGATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCActgc                       >  2:5864289/1‑93 (MQ=255)
                                                                   tGGGTGATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCActgc                       <  1:3143784/93‑1 (MQ=255)
                                                                   tGGGTGATCTTCATCGACGGACACATCGGACTACGAGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCActgc                       <  2:2610930/93‑1 (MQ=255)
                                                                     ggTGATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCActgctg                     <  1:2429887/93‑1 (MQ=255)
                                                                     ggTGATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCActgctg                     <  1:5636862/93‑1 (MQ=255)
                                                                     ggTGATCTTCATCGACGGACACATCGGACTACGAGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCActgctg                     <  2:1236247/93‑1 (MQ=255)
                                                                      gTGATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTgg                    >  2:6267272/1‑93 (MQ=255)
                                                                        gATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCg                  >  2:2854914/1‑93 (MQ=255)
                                                                        gATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCg                  >  1:5567437/1‑93 (MQ=255)
                                                                        gATCTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCg                  >  1:1240996/1‑93 (MQ=255)
                                                                           cTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTAc               >  2:7152370/1‑93 (MQ=255)
                                                                           cTTCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTAc               >  1:6378708/1‑93 (MQ=255)
                                                                            ttCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTAcc              >  1:5336688/1‑93 (MQ=255)
                                                                            ttCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTAcc              <  1:1343151/93‑1 (MQ=255)
                                                                            ttCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTAcc              >  1:2613495/1‑93 (MQ=255)
                                                                             tCATCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCg             >  2:3644240/1‑93 (MQ=255)
                                                                               aTCGACGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCGCa           >  1:6694252/1‑93 (MQ=255)
                                                                                   aCGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCGCAATgg       <  2:6097771/93‑1 (MQ=255)
                                                                                    cGGACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCGCAATggg      <  1:6791863/93‑1 (MQ=255)
                                                                                      gACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCGCAATGGGGa    >  2:2126255/1‑93 (MQ=255)
                                                                                      gACACGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACACCATCGCACCGCGCCACTGCTGGCGTACCGCAATGGGGa    <  2:2215978/93‑1 (MQ=255)
                                                                                      gACACATCGGACTACGAGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCGCAATGGGGa    >  2:4963565/1‑93 (MQ=255)
                                                                                      gACACATCGGACTACGAGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCGCAATGGGGa    <  2:168324/93‑1 (MQ=255)
                                                                                       acacGTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCGCAATGGGGAt   <  1:4462241/93‑1 (MQ=255)
                                                                                        cacgTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCGCAATGGGGATc  <  1:429427/93‑1 (MQ=255)
                                                                                        cacgTCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCGCAATGGGGATc  <  1:5226969/93‑1 (MQ=255)
                                                                                        cacaTCGGACTACGAGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCGCAATGGGGATc  <  1:5185293/93‑1 (MQ=255)
                                                                                           |                                                                                         
ATCCGCCAACAAACGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCGCAATGGGGATC  >  NC_000913/1764841‑1765021

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