Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 638,137 A→G intergenic (‑281/‑31) dsbG ← / → ahpC thiol:disulfide interchange protein/alkyl hydroperoxide reductase C22 subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913638,1370AG98.2% 1132.4 / 4.6 334intergenic (‑281/‑31)dsbG/ahpCthiol:disulfide interchange protein/alkyl hydroperoxide reductase C22 subunit
Reads supporting (aligned to +/- strand):  ref base A (5/1);  new base G (153/175);  total (158/176)
Fisher's exact test for biased strand distribution p-value = 1.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 339 total aligned reads displayed.

CACAAAGTTGTCTTACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTACATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAGGCATTCAAAAACGGCGA  >  NC_000913/638047‑638226
                                                                                          |                                                                                         
  caAAGTTGTCTTACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCAtc                                                                                       >  1:8697804/1‑93 (MQ=255)
   aaaGTTGTCTTACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGTATTAGCCGAATCGGTAAAAATTGGTTACCTTGCAtct                                                                                      >  1:5012274/1‑93 (MQ=255)
   aaaGTTGTCTTACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCAtct                                                                                      >  1:538977/1‑93 (MQ=255)
    aaGTTGTCTTACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCAtctc                                                                                     <  1:8077054/93‑1 (MQ=255)
      gTTGTCTTACGAAGGTTGTAAGGTAAAACTTATCCATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCAt                                                                                   <  2:833180/93‑1 (MQ=255)
       ttGTCTTACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATc                                                                                  <  1:4685814/93‑1 (MQ=255)
         gTCTTACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGa                                                                                >  2:6027530/1‑93 (MQ=255)
          tCTTACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGaa                                                                               >  2:766180/1‑93 (MQ=255)
            ttACGAAGGTTGTAAGGTAAACCTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGaaaa                                                                             <  2:5473074/93‑1 (MQ=255)
            ttACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTTGTTACCTTGCATCTCATCGaaaa                                                                             >  1:3657965/1‑93 (MQ=255)
             tACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAac                                                                            >  2:2836539/1‑93 (MQ=255)
              aCGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAaca                                                                           <  2:1337347/93‑1 (MQ=255)
              aCGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAaca                                                                           <  2:8072938/93‑1 (MQ=255)
               cGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAacac                                                                          <  1:336699/93‑1 (MQ=255)
               cGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAacac                                                                          >  2:3639219/1‑93 (MQ=255)
               cGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAacac                                                                          <  1:1907735/93‑1 (MQ=255)
                gAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACg                                                                         >  2:2353820/1‑93 (MQ=255)
                  aGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACgga                                                                       <  2:1222648/93‑1 (MQ=255)
                  aGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACgga                                                                       <  1:8510890/93‑1 (MQ=255)
                  aGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACgga                                                                       >  1:4597491/1‑93 (MQ=255)
                  aGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACgga                                                                       >  1:392257/1‑93 (MQ=255)
                   ggTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACggag                                                                      <  2:5048022/93‑1 (MQ=255)
                   ggTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACggag                                                                      <  1:7338263/93‑1 (MQ=255)
                   ggTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACggag                                                                      <  1:6284168/93‑1 (MQ=255)
                   ggTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACggag                                                                      >  1:7116123/1‑93 (MQ=255)
                    gTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTTGTTACCTTGCATCTCATCGAAAACACggagg                                                                     >  1:3565923/1‑93 (MQ=255)
                       gTAAGGTAACACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAACACACggaggacg                                                                  >  1:4137471/1‑91 (MQ=255)
                       gTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAg                                                                  >  1:5722759/1‑93 (MQ=255)
                       gTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAg                                                                  <  2:6672561/93‑1 (MQ=255)
                       gTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAg                                                                  >  1:2750149/1‑93 (MQ=255)
                        tAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGt                                                                 >  2:5915146/1‑93 (MQ=255)
                          aGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGtat                                                               <  2:5170363/93‑1 (MQ=255)
                          aGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGtat                                                               >  1:8158201/1‑93 (MQ=255)
                          aGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGtat                                                               >  2:3515227/1‑93 (MQ=255)
                           ggTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGtata                                                              >  2:2231501/1‑93 (MQ=255)
                           ggTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGtata                                                              <  1:2819580/93‑1 (MQ=255)
                             tAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGa                                                            <  2:682079/93‑1 (MQ=255)
                             tAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGa                                                            <  2:2992698/93‑1 (MQ=255)
                             tAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGa                                                            <  1:5341275/93‑1 (MQ=255)
                               aaaCTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGATGTATAGATg                                                          >  2:7856708/1‑93 (MQ=255)
                                aaCTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGt                                                         >  2:624863/1‑93 (MQ=255)
                                 aCTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTc                                                        >  1:3158399/1‑93 (MQ=255)
                                  cTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTcc                                                       >  1:6717486/1‑93 (MQ=255)
                                   ttATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCt                                                      >  2:2144948/1‑93 (MQ=255)
                                    tATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAAGACGGAGGAAGTATAGATGTCCtt                                                     >  2:4756197/1‑93 (MQ=255)
                                    tATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCtt                                                     >  2:6722596/1‑93 (MQ=255)
                                    tATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCtt                                                     <  2:4372072/93‑1 (MQ=255)
                                    tATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTACATCTCATCGAAAACACGGAGGAAGTATAGATGTCCtt                                                     <  1:1385706/93‑1 (MQ=255)
                                     aTCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTg                                                    <  2:1109395/93‑1 (MQ=255)
                                      tCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGa                                                   <  1:4765942/93‑1 (MQ=255)
                                      tCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGa                                                   <  2:3902690/93‑1 (MQ=255)
                                      tCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGa                                                   <  2:5284323/93‑1 (MQ=255)
                                      tCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGa                                                   <  2:1775975/93‑1 (MQ=255)
                                       cGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGAt                                                  <  1:3463531/93‑1 (MQ=255)
                                          tttGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTaa                                               >  2:419157/1‑93 (MQ=255)
                                           ttGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAac                                              <  1:1655806/93‑1 (MQ=255)
                                              aTAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACAcc                                           <  2:4024827/93‑1 (MQ=255)
                                               tAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCa                                          <  1:2412221/93‑1 (MQ=255)
                                               tAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTACATCTCATCGAAAACACGGAGGAAGAATAGATGTCCTTGATTAACACCa                                          >  1:7076554/1‑93 (MQ=255)
                                                aaTGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCaa                                         <  1:6979515/93‑1 (MQ=255)
                                                aaTGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCaa                                         <  1:4268071/93‑1 (MQ=255)
                                                 aTGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCaaa                                        >  1:216171/1‑93 (MQ=255)
                                                    gAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAAtt                                     >  2:5773303/1‑93 (MQ=255)
                                                    gAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAAtt                                     >  1:91398/1‑93 (MQ=255)
                                                     aaaCGCATTAGCCGAATCGGCAAAAATTGGTTACCTTACATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTa                                    >  1:5995198/1‑93 (MQ=255)
                                                        cGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAAc                                 <  2:253804/93‑1 (MQ=255)
                                                        cGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAAc                                 <  1:55894/93‑1 (MQ=255)
                                                         gCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAAcc                                <  2:7532510/93‑1 (MQ=255)
                                                          cATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCt                               >  2:3243228/1‑93 (MQ=255)
                                                          cATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCt                               >  2:7582643/1‑93 (MQ=255)
                                                          cATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCt                               <  2:2031945/93‑1 (MQ=255)
                                                            ttAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCttt                             <  1:7980189/93‑1 (MQ=255)
                                                             tAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACTAAAATTAAACCtttt                            >  1:5923401/1‑93 (MQ=255)
                                                             tAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCtttt                            <  1:3067720/93‑1 (MQ=255)
                                                                  gAATCTGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTaaaaa                       <  2:4673461/93‑1 (MQ=255)
                                                                  gAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTaaaaa                       >  2:470557/1‑93 (MQ=255)
                                                                  gAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTaaaaa                       >  2:6306331/1‑93 (MQ=255)
                                                                  gAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTaaaaa                       >  1:1289924/1‑93 (MQ=255)
                                                                  gAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTaaaaa                       >  2:8497444/1‑93 (MQ=255)
                                                                   aaTCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAAc                      >  1:6047084/1‑93 (MQ=255)
                                                                     tCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCa                    <  1:3209325/93‑1 (MQ=255)
                                                                      cGGCTAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATAAAACCTTTTAAAAACCAg                   >  1:6955295/1‑93 (MQ=255)
                                                                      cGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAg                   <  1:8595370/93‑1 (MQ=255)
                                                                      cGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAg                   <  2:5517876/93‑1 (MQ=255)
                                                                       ggCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAgg                  <  1:4973150/93‑1 (MQ=255)
                                                                        gCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAGGc                 <  2:8209241/93‑1 (MQ=255)
                                                                          aaaaaTTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAGGCAt               <  1:6434708/93‑1 (MQ=255)
                                                                          aaaaaTTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAGGCAt               >  2:379036/1‑93 (MQ=255)
                                                                            aaaTTGGTTACCTTGCATGTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTGACACCAAAATTAAACCTTTTAAAAACCAGGCATTc             <  1:6196346/93‑1 (MQ=255)
                                                                             aaTTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAGGCATTCa            <  2:8287544/93‑1 (MQ=255)
                                                                              aTTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAGGCATTCaa           <  1:6851917/93‑1 (MQ=255)
                                                                               ttGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAGGCATTCaaa          <  1:6538244/93‑1 (MQ=255)
                                                                               ttGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAGGCATTCaaa          >  1:7532129/1‑93 (MQ=255)
                                                                                  gTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAGGCATTCAAAAAc       <  2:7236736/93‑1 (MQ=255)
                                                                                  gTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAGGCATTCAAAAAc       >  1:5795367/1‑93 (MQ=255)
                                                                                    tACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAGGCATTCAAAAACgg     >  2:8734377/1‑93 (MQ=255)
                                                                                     aCCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAGGCATTCAAAAACGGc    >  1:4878539/1‑93 (MQ=255)
                                                                                       cTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAGGCATTCAAAAACGGCGa  >  2:6985967/1‑93 (MQ=255)
                                                                                       cTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAGGCATTCAAAAACGGCGa  >  2:153750/1‑93 (MQ=255)
                                                                                       cTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAGGCATTCAAAAACGGCGa  >  1:8353040/1‑93 (MQ=255)
                                                                                          |                                                                                         
CACAAAGTTGTCTTACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTACATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAGGCATTCAAAAACGGCGA  >  NC_000913/638047‑638226

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