Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,791,583 T→C intergenic (‑258/‑31) htrL ← / → rfaD HtrL/ADP‑L‑glycero‑D‑mannoheptose‑6‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,791,5830TC96.4% 497.0 / 5.9 165intergenic (‑258/‑31)htrL/rfaDHtrL/ADP‑L‑glycero‑D‑mannoheptose‑6‑epimerase
Reads supporting (aligned to +/- strand):  ref base T (5/1);  new base C (107/52);  total (112/53)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.89e-01

Only 100 of 169 total aligned reads displayed.

GTTGCAATTAGCATCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGTCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAAGCCCTGAA  >  NC_000913/3791494‑3791672
                                                                                         |                                                                                         
 ttgCAATTAGCATCCTTGCGCCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCctct                                                                                       <  2:3307848/93‑1 (MQ=255)
 ttgCAATTAGCATCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCctct                                                                                       <  2:715427/93‑1 (MQ=255)
   gCAATTAGCATCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGa                                                                                     >  2:2290206/1‑93 (MQ=255)
     aaTTAGCATCCTTGCACCTCTATGTAAAGGGCTGACGGGATTCGGATGGGCTGGTATGATTACAGCCCTTCGTGTCTGAGATTGCCTCTGagc                                                                                   >  1:2843067/1‑91 (MQ=255)
     aaTTAGCATCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACt                                                                                   >  1:3625782/1‑93 (MQ=255)
     aaTTAGCATCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACt                                                                                   >  2:779309/1‑93 (MQ=255)
      aTTAGCATCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTc                                                                                  >  1:913009/1‑93 (MQ=255)
      aTTAGCATCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTc                                                                                  >  2:2044015/1‑93 (MQ=255)
       ttAGCATCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTcc                                                                                 >  2:1095300/1‑93 (MQ=255)
          gCATCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATa                                                                              >  1:3063694/1‑93 (MQ=255)
          gCATCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATa                                                                              >  2:2365924/1‑93 (MQ=255)
            aTCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTTGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAAt                                                                            >  1:726954/1‑93 (MQ=255)
            aTCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAAt                                                                            >  2:1005177/1‑93 (MQ=255)
            aTCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAAt                                                                            >  2:4735424/1‑93 (MQ=255)
            aTCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAAt                                                                            <  2:31921/93‑1 (MQ=255)
             tCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAAtt                                                                           >  2:1376374/1‑93 (MQ=255)
             tCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGAGTGCCTCTGACTCCATAAtt                                                                           >  2:3618730/1‑93 (MQ=255)
               cTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCg                                                                         >  1:312101/1‑93 (MQ=255)
                  gCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAg                                                                      >  1:3314914/1‑93 (MQ=255)
                   cACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAgg                                                                     >  2:2626465/1‑93 (MQ=255)
                   cACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAgg                                                                     >  1:4578896/1‑93 (MQ=255)
                   cACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAgg                                                                     >  1:1560813/1‑93 (MQ=255)
                   cACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAgg                                                                     >  1:323511/1‑93 (MQ=255)
                    aCCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGt                                                                    >  2:1186965/1‑93 (MQ=255)
                    aCCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGt                                                                    >  2:826806/1‑93 (MQ=255)
                    aCCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGt                                                                    >  1:1790551/1‑93 (MQ=255)
                     ccTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGtt                                                                   >  1:3908786/1‑93 (MQ=255)
                      ctctATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTa                                                                  <  1:473934/93‑1 (MQ=255)
                       tctATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCT‑‑GATTGCCTCTGACTCCATAATTCGAAGGTTACAg                                                               >  1:1007173/1‑93 (MQ=255)
                        ctATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGTCTCTGACTCCATAATTCGAAGGTTACa                                                                >  1:4015206/1‑93 (MQ=255)
                         tATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAg                                                               >  1:889451/1‑93 (MQ=255)
                           tGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGtt                                                             >  1:1685798/1‑93 (MQ=255)
                           tGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGtt                                                             <  2:285582/93‑1 (MQ=255)
                           tGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGtt                                                             >  1:275615/1‑93 (MQ=255)
                            gTAAAGGGCTGAAGGGATTCGGCTGTGATGGGATGATTAAAGACATTCGAGTATGAGCTTGACTCTGACTCCATAATTCGAGGGTTTCAGttt                                                            >  2:2396596/1‑92 (MQ=255)
                            gTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTa                                                            <  1:1737074/93‑1 (MQ=255)
                             tAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTAt                                                           >  2:3570295/1‑93 (MQ=255)
                             tAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTAt                                                           <  1:1397120/93‑1 (MQ=255)
                              aaaGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATg                                                          >  2:991448/1‑93 (MQ=255)
                                aGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGat                                                        >  1:800097/1‑93 (MQ=255)
                                aGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGat                                                        >  1:3899325/1‑93 (MQ=255)
                                 gggCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGatc                                                       >  1:2659902/1‑93 (MQ=255)
                                  ggCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGatca                                                      <  2:961423/93‑1 (MQ=255)
                                    cTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGatcatc                                                    <  1:2680489/93‑1 (MQ=255)
                                     tGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTTAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCg                                                   >  2:4314471/1‑93 (MQ=255)
                                     tGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCg                                                   >  2:1798712/1‑93 (MQ=255)
                                     tGAAGAGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCg                                                   >  1:66960/1‑93 (MQ=255)
                                      gAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGt                                                  >  1:2199901/1‑93 (MQ=255)
                                      gAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTGTCAGTTATGATCATCGt                                                  >  1:4225656/1‑93 (MQ=255)
                                       aaGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGtt                                                 >  2:3356357/1‑93 (MQ=255)
                                         gggATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTAc                                               <  1:4382762/93‑1 (MQ=255)
                                         gggATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTAc                                               <  2:4154709/93‑1 (MQ=255)
                                         gggATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTAc                                               >  1:2744095/1‑93 (MQ=255)
                                          ggATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTAcc                                              <  2:2060779/93‑1 (MQ=255)
                                          ggATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTAcc                                              >  1:3682941/1‑93 (MQ=255)
                                           gATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACcg                                             >  2:4491167/1‑93 (MQ=255)
                                            atgcGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACcgg                                            <  1:1642973/90‑1 (MQ=255)
                                            aTTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACcgg                                            <  1:1642974/93‑1 (MQ=255)
                                            aTTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACcgg                                            >  1:3796981/1‑93 (MQ=255)
                                             ttCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACcggc                                           <  1:4185419/93‑1 (MQ=255)
                                                 gATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGgcggcg                                       >  1:3596294/1‑93 (MQ=255)
                                                 gATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGgcggcg                                       >  2:491417/1‑93 (MQ=255)
                                                  aTGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGgcgc                                      >  2:1708120/1‑93 (MQ=255)
                                                   tgtgATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGcgcg                                     >  2:1600295/1‑93 (MQ=255)
                                                   tgtgATGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGcgcg                                     <  2:452980/93‑1 (MQ=255)
                                                       aTGGTATGATTACCGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGCAGGTTACAGTTATGATCATCGTTACCGGTGGCGCGGGCt                                 >  1:2396556/1‑93 (MQ=255)
                                                       aTGGTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCt                                 <  2:383884/93‑1 (MQ=255)
                                                         ggTATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCttt                               >  2:1420714/1‑93 (MQ=255)
                                                           tATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTAt                             <  2:3494661/93‑1 (MQ=255)
                                                           tATGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTAt                             >  1:4624560/1‑93 (MQ=255)
                                                            aTGATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATc                            <  1:2130654/93‑1 (MQ=255)
                                                              gATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCgg                          <  2:2550466/93‑1 (MQ=255)
                                                              gATTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCgg                          <  2:2034535/93‑1 (MQ=255)
                                                               aTTACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGgc                         >  2:539261/1‑93 (MQ=255)
                                                                ttACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGgca                        >  2:303369/1‑93 (MQ=255)
                                                                 tACAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGGTTTATCGgcag                       >  1:1610640/1‑93 (MQ=255)
                                                                  aCAGACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGgcagc                      <  2:3994922/93‑1 (MQ=255)
                                                                     gACATTCGTGTCTGAGATTGTCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAAc                   >  1:2343193/1‑93 (MQ=255)
                                                                     gACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCCAc                   >  1:3032860/1‑93 (MQ=255)
                                                                     gACATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAAc                   >  2:2934245/1‑93 (MQ=255)
                                                                      aCATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACa                  <  1:27228/93‑1 (MQ=255)
                                                                       cATTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACAt                 >  1:3508865/1‑93 (MQ=255)
                                                                        aTTCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATc                <  2:1865465/93‑1 (MQ=255)
                                                                         ttCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCg               >  1:1013875/1‑93 (MQ=255)
                                                                          tCGTGTCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGt              <  2:3141510/93‑1 (MQ=255)
                                                                            gtgtCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTa            >  2:852699/1‑93 (MQ=255)
                                                                             tgtCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTaa           >  2:3815192/1‑93 (MQ=255)
                                                                              gtCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCCACATCGTTaaa          >  2:1966885/1‑93 (MQ=255)
                                                                              gtCTGAGATGTCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTaaa          <  2:2530028/93‑1 (MQ=255)
                                                                               tCTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAAg         >  2:4713994/1‑93 (MQ=255)
                                                                                cTGAGATTGTCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAAGc        >  2:4506798/1‑93 (MQ=255)
                                                                                cTGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAAGc        >  2:1251321/1‑93 (MQ=255)
                                                                                 tGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAAGcc       <  2:2133009/93‑1 (MQ=255)
                                                                                 tGAGATTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAAGcc       <  1:2057076/93‑1 (MQ=255)
                                                                                  gagaTTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGGGGGCTTTATCGGCAGAAACATCGTTAACGccc      >  2:1646673/1‑93 (MQ=255)
                                                                                  gagaTTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAAGccc      >  2:1646674/1‑93 (MQ=255)
                                                                                    gaTTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAAGCCCTg    >  2:2845712/1‑93 (MQ=255)
                                                                                     aTTGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAAGCCCTGa   >  2:3257694/1‑93 (MQ=255)
                                                                                      ttGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAAGCCCTGaa  <  1:837750/93‑1 (MQ=255)
                                                                                      ttGCCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAAGCCCTGaa  <  1:1058151/93‑1 (MQ=255)
                                                                                         |                                                                                         
GTTGCAATTAGCATCCTTGCACCTCTATGTAAAGGGCTGAAGGGATTCGGATGTGATGGTATGATTACAGACATTCGTGTCTGAGATTGTCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAAGCCCTGAA  >  NC_000913/3791494‑3791672

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