Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,793,037 A→G Q161R (CAG→CGG)  rfaF → ADP‑heptose‑‑lps heptosyltransferase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,793,0370AG99.8% 1432.0 / ‑5.4 454Q161R (CAG→CGG) rfaFADP‑heptose‑‑lps heptosyltransferase II
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (220/233);  total (221/233)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.14e-01

Only 100 of 471 total aligned reads displayed.

TGGTGGAACGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCGGGTGCGGAGT  >  NC_000913/3792947‑3793126
                                                                                          |                                                                                         
tggtggAGCGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCggg                                                                                         >  2:7759875/1‑93 (MQ=255)
tggtggAACGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCggg                                                                                         <  2:4460910/93‑1 (MQ=255)
tggtggAACGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCggg                                                                                         <  2:2257602/93‑1 (MQ=255)
 ggtggAACGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGt                                                                                        >  2:7353568/1‑93 (MQ=255)
  gtggAACGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTg                                                                                       <  2:3253962/93‑1 (MQ=255)
     gAACGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGc                                                                                    <  1:1903464/93‑1 (MQ=255)
     gAACGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGc                                                                                    <  1:4553145/93‑1 (MQ=255)
      aaCGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAAATCTGCCACAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCg                                                                                   >  1:3747075/1‑93 (MQ=255)
       aCGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGa                                                                                  <  2:3655113/93‑1 (MQ=255)
       aCGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGa                                                                                  <  1:2385743/93‑1 (MQ=255)
        cGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGaa                                                                                 <  2:3111716/93‑1 (MQ=255)
         gCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAg                                                                                <  1:6951945/93‑1 (MQ=255)
         gCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAg                                                                                >  1:1711803/1‑93 (MQ=255)
          cTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAgg                                                                               >  2:4778030/1‑93 (MQ=255)
          cTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAgg                                                                               <  2:2644690/93‑1 (MQ=255)
           tatataGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAggg                                                                              >  2:3403179/1‑92 (MQ=255)
           tatataGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAggg                                                                              >  1:2207066/1‑92 (MQ=255)
            atataGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTg                                                                             <  2:2329441/93‑1 (MQ=255)
             tataGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGa                                                                            >  1:1337733/1‑93 (MQ=255)
             tataGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGa                                                                            <  1:3304020/93‑1 (MQ=255)
                 gcgcGGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGaaaaa                                                                        <  1:3663171/93‑1 (MQ=255)
                   gcTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATc                                                                      <  1:2550173/93‑1 (MQ=255)
                     tGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCat                                                                    >  1:1098975/1‑93 (MQ=255)
                     tGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCat                                                                    >  2:5325367/1‑93 (MQ=255)
                      ggCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCata                                                                   >  1:1182568/1‑93 (MQ=255)
                      ggCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCata                                                                   >  1:3938825/1‑93 (MQ=255)
                        cttATGACAAAGACATTATGAGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCatatg                                                                 <  2:1294519/91‑2 (MQ=255)
                        ccTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCAtata                                                                 <  1:345435/93‑1 (MQ=255)
                         cTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATAc                                                                >  2:810655/1‑93 (MQ=255)
                         cTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATAc                                                                <  1:5970682/93‑1 (MQ=255)
                          tATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATAcc                                                               >  1:5484299/1‑93 (MQ=255)
                          tATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATAcc                                                               >  2:5264386/1‑93 (MQ=255)
                           aTGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCt                                                              >  2:6134670/1‑93 (MQ=255)
                           aTGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCt                                                              >  2:5570540/1‑93 (MQ=255)
                            tGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTg                                                             >  2:6702351/1‑93 (MQ=255)
                               cAAAGGCATTATGCGCAGAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGTGCTAAGGTGAAAAATCATATACCTGTaa                                                          >  2:5446082/1‑93 (MQ=255)
                               cAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTaa                                                          >  1:4497164/1‑93 (MQ=255)
                               cAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTaa                                                          >  2:1839985/1‑93 (MQ=255)
                                aaaGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAAt                                                         >  2:5728725/1‑93 (MQ=255)
                                aaaGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAAt                                                         >  1:3801012/1‑93 (MQ=255)
                                aaaGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAAt                                                         <  2:2832120/93‑1 (MQ=255)
                                aaaGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAAt                                                         >  2:7216275/1‑93 (MQ=255)
                                aaaGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAAt                                                         <  1:4170249/93‑1 (MQ=255)
                                 aaGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATc                                                        >  1:1381421/1‑93 (MQ=255)
                                   ggCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGGGAAAAATCATATACCTGTAATCaa                                                      >  2:5060127/1‑93 (MQ=255)
                                     cATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAAtt                                                    <  2:2862960/93‑1 (MQ=255)
                                      aTTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAAttt                                                   >  1:6813189/1‑93 (MQ=255)
                                      aTTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAAttt                                                   >  2:1579697/1‑93 (MQ=255)
                                       ttATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAAtttt                                                  >  2:437772/1‑93 (MQ=255)
                                       ttATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAAtttt                                                  >  1:7400785/1‑93 (MQ=255)
                                        tATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTc                                                 >  2:6522256/1‑93 (MQ=255)
                                        tATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTc                                                 <  1:6064508/93‑1 (MQ=255)
                                           gcgcACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCt                                              <  2:871241/93‑1 (MQ=255)
                                           gcgcACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCt                                              >  2:7557765/1‑93 (MQ=255)
                                            cgcACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCtt                                             <  2:4228875/93‑1 (MQ=255)
                                            cgcACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCtt                                             <  2:4053007/93‑1 (MQ=255)
                                             gcacagcacaAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCttt                                            <  1:838269/93‑1 (MQ=255)
                                             gcacagcacaAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCttt                                            >  2:4550931/1‑93 (MQ=255)
                                               acagcacaAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTtca                                          >  2:3914625/1‑93 (MQ=255)
                                                cagcacaAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTtcat                                         <  1:4381348/93‑1 (MQ=255)
                                                 agcacaAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGTGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTtcatc                                        >  2:4673028/1‑93 (MQ=255)
                                                  gcacaAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTtcatca                                       <  1:2940204/93‑1 (MQ=255)
                                                  gcacaAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTtcatca                                       <  1:723203/93‑1 (MQ=255)
                                                     caAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGaa                                    >  2:3526591/1‑93 (MQ=255)
                                                     caAGATCTGCCGCAGCCATTGTTATGGCCGCAGCTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGaa                                    >  1:2880166/1‑93 (MQ=255)
                                                      aaGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAAc                                   <  1:3055948/93‑1 (MQ=255)
                                                       aGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACg                                  >  2:1019992/1‑93 (MQ=255)
                                                       aGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACg                                  <  1:7097558/93‑1 (MQ=255)
                                                         aTCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTc                                <  1:2104857/93‑1 (MQ=255)
                                                         aTCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTc                                <  1:6492719/93‑1 (MQ=255)
                                                           cTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCg                              >  1:1483696/1‑93 (MQ=255)
                                                             gCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCgat                            >  1:5422402/1‑93 (MQ=255)
                                                             gCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCgat                            >  1:6154529/1‑93 (MQ=255)
                                                             gCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGAGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCgat                            >  1:175835/1‑93 (MQ=255)
                                                                gCAGCCATTGTGATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCgatgat                         <  1:6367887/93‑1 (MQ=255)
                                                                 cAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGAtt                        <  2:137386/93‑1 (MQ=255)
                                                                 cAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGAtt                        <  1:5220740/93‑1 (MQ=255)
                                                                 cAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGAtt                        <  1:1543773/93‑1 (MQ=255)
                                                                  aGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTg                       <  2:7096192/93‑1 (MQ=255)
                                                                   gCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTgg                      >  1:4052788/1‑93 (MQ=255)
                                                                     cATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGtt                    >  1:5637824/1‑93 (MQ=255)
                                                                     cATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGGAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGtt                    >  2:1651967/1‑93 (MQ=255)
                                                                      aTTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGttt                   <  1:4861943/93‑1 (MQ=255)
                                                                      aTTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGttt                   <  2:6307925/93‑1 (MQ=255)
                                                                       ttGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGtttt                  >  2:4967556/1‑93 (MQ=255)
                                                                       ttGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGtttt                  <  1:5013261/93‑1 (MQ=255)
                                                                       ttGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTGTCATCAGAACGTCCGATGATTGGtttt                  >  1:7569840/1‑93 (MQ=255)
                                                                          ttATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGc               >  1:5763759/1‑93 (MQ=255)
                                                                          ttATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGc               >  2:2431213/1‑93 (MQ=255)
                                                                           tATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGcc              <  1:2614881/93‑1 (MQ=255)
                                                                            aTGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGccc             >  1:2257388/1‑93 (MQ=255)
                                                                             tGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCg            <  2:6376532/93‑1 (MQ=255)
                                                                              ggCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCgg           >  1:3014580/1‑93 (MQ=255)
                                                                              ggCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCgg           <  1:312892/93‑1 (MQ=255)
                                                                               gCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCggg          <  2:587199/93‑1 (MQ=255)
                                                                               gCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCggg          >  1:6631479/1‑93 (MQ=255)
                                                                               gCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCggg          >  2:1140844/1‑93 (MQ=255)
                                                                                ccGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCGGGt         <  1:4287146/93‑1 (MQ=255)
                                                                                   cAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCGGGTGCg      >  1:3225335/1‑93 (MQ=255)
                                                                                      ttGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCGGGTGCGGAg   <  1:785918/93‑1 (MQ=255)
                                                                                          |                                                                                         
TGGTGGAACGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCGGGTGCGGAGT  >  NC_000913/3792947‑3793126

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: ‑