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,1370AG88.1% 781.9 / 90.0 304intergenic (‑281/‑31)dsbG/ahpCthiol:disulfide interchange protein/alkyl hydroperoxide reductase C22 subunit
Reads supporting (aligned to +/- strand):  ref base A (16/20);  new base G (106/161);  total (123/181)
Fisher's exact test for biased strand distribution p-value = 5.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.04e-01

Only 100 of 307 total aligned reads displayed.

TTCACAAAGTTGTCTTACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTACATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAGGCATTCAAAAACGGCGAATT  >  NC_000913/638045‑638229
                                                                                            |                                                                                            
  cacaAAGTTGTCTTACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCa                                                                                            >  2:5013899/1‑93 (MQ=255)
  cacaAAGTTGTCTTACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCa                                                                                            <  2:4126875/93‑1 (MQ=255)
  cacaAAGTTGTCTTACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCa                                                                                            >  2:3454350/1‑93 (MQ=255)
  cacaAAGTTGTCTTACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTACa                                                                                            <  1:6669449/93‑1 (MQ=255)
   acaAGGTTGTCTTACGAAGGTTGTAAGTTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTACAt                                                                                           >  2:6421200/1‑93 (MQ=255)
    caAGGTTGTCTTACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTACAtc                                                                                          <  2:5484948/93‑1 (MQ=255)
      aaGTTGTCTTACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCAtctc                                                                                        <  2:4951116/93‑1 (MQ=255)
      aaGTTGTCTTACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCAtctc                                                                                        <  2:1151491/93‑1 (MQ=255)
       aGTTGTCTTACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCa                                                                                       >  1:5882849/1‑93 (MQ=255)
       aGTTGTCTTACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCa                                                                                       <  2:1057230/93‑1 (MQ=255)
       aGTTGTCTTACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCa                                                                                       >  2:3755500/1‑93 (MQ=255)
       aGTTGTCTTACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCa                                                                                       >  2:1311434/1‑93 (MQ=255)
        gTTGTCTTACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCAt                                                                                      <  1:4474962/93‑1 (MQ=255)
        gTTGTCTTACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCAt                                                                                      <  1:228654/93‑1 (MQ=255)
        gTTGTCTTACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTACATCTCAt                                                                                      >  1:2617167/1‑93 (MQ=255)
         ttGTCTTACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATc                                                                                     >  1:211875/1‑93 (MQ=255)
            tCTTACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGaa                                                                                  >  2:1687346/1‑93 (MQ=255)
              ttACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGaaaa                                                                                <  1:854525/93‑1 (MQ=255)
              ttACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTACATCTCATCGaaaa                                                                                >  2:5266391/1‑93 (MQ=255)
              ttACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTACATCTCATCGaaaa                                                                                <  1:3965144/93‑1 (MQ=255)
               tACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTACATCTCATCGAAAac                                                                               >  1:114422/1‑93 (MQ=255)
                aCGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAaca                                                                              <  2:1875245/93‑1 (MQ=255)
                aCGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAaca                                                                              <  2:4071054/93‑1 (MQ=255)
                 cGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAacac                                                                             <  2:3083814/93‑1 (MQ=255)
                  gAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACg                                                                            <  1:6903756/93‑1 (MQ=255)
                  gAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTACATCTCATCGAAAACACg                                                                            <  1:3850659/93‑1 (MQ=255)
                   aaGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACgg                                                                           >  2:5358099/1‑93 (MQ=255)
                    aGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACgga                                                                          <  1:2542505/93‑1 (MQ=255)
                    aGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACgga                                                                          >  2:2829250/1‑93 (MQ=255)
                        tGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGaa                                                                      >  2:5074588/1‑93 (MQ=255)
                         gTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAg                                                                     <  2:155425/93‑1 (MQ=255)
                         gTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAg                                                                     <  1:6263968/93‑1 (MQ=255)
                          tAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGt                                                                    <  1:4881922/93‑1 (MQ=255)
                          tAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGt                                                                    >  1:6398861/1‑93 (MQ=255)
                           aaGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGta                                                                   >  2:6776613/1‑93 (MQ=255)
                           aaGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGta                                                                   <  1:4276785/93‑1 (MQ=255)
                            aGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGtat                                                                  >  2:1959323/1‑93 (MQ=255)
                             ggTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGtata                                                                 <  2:5163727/93‑1 (MQ=255)
                               tAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGa                                                               >  2:4840142/1‑93 (MQ=255)
                               tAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGa                                                               <  2:7070215/93‑1 (MQ=255)
                               tAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGa                                                               <  2:7492120/93‑1 (MQ=255)
                                aaaaCTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGAt                                                              <  1:6802160/93‑1 (MQ=255)
                                   aCTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTc                                                           <  2:1471747/93‑1 (MQ=255)
                                    cTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTcc                                                          <  2:2370916/93‑1 (MQ=255)
                                      tATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCtt                                                        >  2:3158921/1‑93 (MQ=255)
                                      tATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCtt                                                        <  2:4362942/93‑1 (MQ=255)
                                      tATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCtt                                                        <  1:1278786/93‑1 (MQ=255)
                                      tATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCtt                                                        >  1:1964571/1‑93 (MQ=255)
                                       aTCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTg                                                       >  2:2177686/1‑93 (MQ=255)
                                        tCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGa                                                      <  1:4648027/93‑1 (MQ=255)
                                        tCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGa                                                      >  1:5720631/1‑93 (MQ=255)
                                         cGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGAt                                                     <  1:3227794/93‑1 (MQ=255)
                                            tttGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGAttga                                                  <  2:4256814/93‑3 (MQ=255)
                                            tttGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTaa                                                  >  2:3368372/1‑93 (MQ=255)
                                             ttGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAac                                                 <  1:2672286/93‑1 (MQ=255)
                                              tGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAaca                                                <  1:4576622/93‑1 (MQ=255)
                                               gATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAacac                                               >  2:7331448/1‑93 (MQ=255)
                                                aTAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGGTGTCCTTGATTAACAcc                                              >  1:6636428/1‑93 (MQ=255)
                                                aTAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACAcc                                              >  2:226544/1‑93 (MQ=255)
                                                 tAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCa                                             >  1:6144969/1‑93 (MQ=255)
                                                 tAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCa                                             >  1:5008449/1‑93 (MQ=255)
                                                  aaTGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCaa                                            >  1:3093570/1‑93 (MQ=255)
                                                  aaTGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCaa                                            >  2:7227649/1‑93 (MQ=255)
                                                  aaTGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCaa                                            <  2:6263061/93‑1 (MQ=255)
                                                    tGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCaaaa                                          <  2:1841465/93‑1 (MQ=255)
                                                    tGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCaaaa                                          <  1:1422796/93‑1 (MQ=255)
                                                    tGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCaaaa                                          <  2:2743709/93‑1 (MQ=255)
                                                     ggAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAAt                                         <  1:6837927/93‑1 (MQ=255)
                                                     ggAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAAt                                         <  1:3307878/93‑1 (MQ=255)
                                                      gAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAAtt                                        <  1:5123918/93‑1 (MQ=255)
                                                          cGCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAAc                                    <  1:189290/93‑1 (MQ=255)
                                                           gCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAAcc                                   >  2:3664160/1‑93 (MQ=255)
                                                           gCATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAAcc                                   >  2:869156/1‑93 (MQ=255)
                                                            cATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCt                                  >  1:4750020/1‑93 (MQ=255)
                                                            cATTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCt                                  >  2:2957883/1‑93 (MQ=255)
                                                             nntAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCtt                                 <  1:3583530/91‑1 (MQ=255)
                                                             aTTAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCtt                                 >  2:5805056/1‑93 (MQ=255)
                                                               tAGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCtttt                               <  1:4518173/93‑1 (MQ=255)
                                                                aGCCGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTa                              >  1:4096538/1‑93 (MQ=255)
                                                                  ccGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTaaa                            <  2:7400698/93‑1 (MQ=255)
                                                                   cGAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTaaaa                           >  1:4026963/1‑93 (MQ=255)
                                                                    gAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTaaaaa                          >  1:7019358/1‑93 (MQ=255)
                                                                    gAATCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTaaaaa                          <  1:340660/93‑1 (MQ=255)
                                                                      aTCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAAcc                        <  1:7351732/93‑1 (MQ=255)
                                                                      aTCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAAcc                        >  2:6048856/1‑93 (MQ=255)
                                                                       tCGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCa                       >  1:233396/1‑93 (MQ=255)
                                                                        cGGCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAg                      <  1:4929574/93‑1 (MQ=255)
                                                                         ggCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAgg                     <  2:2615778/93‑1 (MQ=255)
                                                                          gCAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAGGc                    >  1:5381507/1‑93 (MQ=255)
                                                                           cAAAAATTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAGGCa                   <  2:937202/93‑1 (MQ=255)
                                                                             aaaaTTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAGGCAtt                 <  1:7122247/93‑1 (MQ=255)
                                                                              aaaTTGGTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAGGCATTc                <  1:1002346/93‑1 (MQ=255)
                                                                                  tGGTTACCTTACATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAGGCATTCaaaa            >  2:6488020/1‑93 (MQ=255)
                                                                                   ggTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAGGCATTCaaaaa           <  1:4121993/93‑1 (MQ=255)
                                                                                    gTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAGGCATTCAAAAAc          <  2:399192/93‑1 (MQ=255)
                                                                                    gTTACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAGGCATTCAAAAAc          <  2:3279435/93‑1 (MQ=255)
                                                                                     ttACCTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAGGCATTCAAAAACg         <  1:2934011/93‑1 (MQ=255)
                                                                                     ttACCTTACATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAGGCATTCAAAAACg         <  1:7243767/93‑1 (MQ=255)
                                                                                        ccTTGCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAGGCATTCAAAAACGGCg      >  2:692791/1‑93 (MQ=255)
                                                                                            aCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAGGCATTCAAAAACGGCGAAtt  >  2:4403630/1‑93 (MQ=255)
                                                                                            |                                                                                            
TTCACAAAGTTGTCTTACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTACATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAGGCATTCAAAAACGGCGAATT  >  NC_000913/638045‑638229

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