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.7% 1011.4 / NA 308Q161R (CAG→CGG) rfaFADP‑heptose‑‑lps heptosyltransferase II
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (129/178);  minor base C (1/0);  total (130/178)
Fisher's exact test for biased strand distribution p-value = 4.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 316 total aligned reads displayed.

TGGTGGAACGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCGGGTGCGGAGT  >  NC_000913/3792947‑3793126
                                                                                          |                                                                                         
  gtggAACGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTg                                                                                       <  2:1741827/93‑1 (MQ=255)
   tggAACGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGa                                                                                      <  1:6710964/93‑1 (MQ=255)
   tggAACGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGa                                                                                      >  1:1597086/1‑93 (MQ=255)
   tggAACGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGa                                                                                      <  1:4672194/93‑1 (MQ=255)
    ggAACGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAg                                                                                     >  2:3873161/1‑93 (MQ=255)
      aaCGCTATATAGCGCTGGCCTATGACAAGGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCg                                                                                   <  1:1025814/93‑1 (MQ=255)
          cTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAgg                                                                               <  1:4310750/93‑1 (MQ=255)
           tatataGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGt                                                                              <  1:457579/93‑1 (MQ=255)
            atataGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCTCAGTTGCGGGTGAGCGAAGGTg                                                                             >  2:3081217/1‑93 (MQ=255)
            atataGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTg                                                                             <  2:144083/93‑1 (MQ=255)
            atataGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTg                                                                             <  2:3407551/93‑1 (MQ=255)
            atataGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTg                                                                             >  2:3814218/1‑93 (MQ=255)
            atataGCGCTGGCCTATGACAAAGGCATTAGGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTg                                                                             <  1:1811510/93‑1 (MQ=255)
             tataGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGa                                                                            >  1:5928459/1‑93 (MQ=255)
             tataGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGa                                                                            <  1:6219898/93‑1 (MQ=255)
             tataGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGa                                                                            >  1:4514404/1‑93 (MQ=255)
             tataGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGa                                                                            <  1:637705/93‑1 (MQ=255)
              ataGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGaa                                                                           <  2:2555010/93‑1 (MQ=255)
               taGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGaaa                                                                          <  1:1508345/93‑1 (MQ=255)
                aGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGaaaa                                                                         <  1:812920/93‑1 (MQ=255)
                aGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGaaaa                                                                         <  2:2743460/93‑1 (MQ=255)
                 gcgcTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGaaaaa                                                                        >  2:3513860/1‑93 (MQ=255)
                  cgcTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAAt                                                                       <  2:830772/93‑1 (MQ=255)
                     tGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCat                                                                    >  2:2895793/1‑93 (MQ=255)
                     tGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCat                                                                    >  2:2371219/1‑93 (MQ=255)
                     tGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCat                                                                    <  2:2322706/93‑1 (MQ=255)
                      ggCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCata                                                                   <  1:3020161/93‑1 (MQ=255)
                       gCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCatat                                                                  >  2:4570133/1‑93 (MQ=255)
                        ccgaTGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAATATCAtata                                                                 <  1:6274126/90‑1 (MQ=255)
                         cTATGACAAAGGCATTATGCGCACAGCACAAGATCTGTCGCAGCCATTGTTATGGCCGCAGTTGCGGGGGAGCGAAGGTGAAAAATCATATAc                                                                >  1:2451572/1‑93 (MQ=255)
                          tATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATAcc                                                               >  1:163771/1‑93 (MQ=255)
                           aTGACAAAGGCATTATGGCCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATCCCt                                                              <  2:4849836/93‑1 (MQ=255)
                           aTGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCt                                                              >  1:3555549/1‑93 (MQ=255)
                              aCAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTa                                                           <  2:3301375/93‑1 (MQ=255)
                              aCAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTa                                                           >  1:2569007/1‑93 (MQ=255)
                                aaaGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAAt                                                         >  1:5331846/1‑93 (MQ=255)
                                 aaGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATc                                                        <  2:1542945/93‑1 (MQ=255)
                                  aGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCa                                                       >  1:5802547/1‑93 (MQ=255)
                                   ggCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCaa                                                      <  1:4087958/93‑1 (MQ=255)
                                   ggCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCaa                                                      <  2:2882678/93‑1 (MQ=255)
                                   ggCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCCGGCGAGCGAAGGTGAAAAATCATATACCTGTAATCaa                                                      >  1:5484332/1‑93 (MQ=255)
                                    gCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAAt                                                     <  1:3227870/93‑1 (MQ=255)
                                    gCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAAt                                                     <  2:5671930/93‑1 (MQ=255)
                                      aTTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAAttt                                                   >  1:956034/1‑93 (MQ=255)
                                       ttATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAAtttt                                                  >  2:1342138/1‑93 (MQ=255)
                                        tATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTc                                                 <  2:5594800/93‑1 (MQ=255)
                                        tATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTc                                                 <  1:2675672/93‑1 (MQ=255)
                                           gcgcACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCt                                              <  2:5114463/93‑1 (MQ=255)
                                             gcacagcacaAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCttt                                            <  2:3152021/93‑1 (MQ=255)
                                              cacagcacaAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTtc                                           >  2:5314760/1‑93 (MQ=255)
                                              cacagcacaAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACATGTAATCAATTTTCGCTTtc                                           >  2:4438872/1‑93 (MQ=255)
                                              cacagcacaAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGGGAGCGAAGGTGAAAAATCATATACATGTAATCAATTTTCGCTTtc                                           >  1:49449/1‑93 (MQ=255)
                                               ccagcacaAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTtca                                          <  2:1179550/92‑1 (MQ=255)
                                               acagcacaAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTtca                                          <  2:4791910/93‑1 (MQ=255)
                                               acagcacaAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTtca                                          <  1:3317562/93‑1 (MQ=255)
                                               acagcacaAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTtca                                          <  2:4966769/93‑1 (MQ=255)
                                               acagcacaAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTtca                                          <  1:5547511/93‑1 (MQ=255)
                                                cagcacaAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGATAAATCATATACCTGTAATCAATTTTCGCTTtcat                                         >  1:5679032/1‑93 (MQ=255)
                                                cagcacaAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTtcat                                         >  1:2075057/1‑93 (MQ=255)
                                                cagcacaAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTtcat                                         <  1:5587643/93‑1 (MQ=255)
                                                 agcacaAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTtcatc                                        <  1:4255621/93‑1 (MQ=255)
                                                 agcacaAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTtcatc                                        <  1:2347333/93‑1 (MQ=255)
                                                  gcacaAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTtcatca                                       >  2:3747223/1‑93 (MQ=255)
                                                  gcacaAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTtcatca                                       <  2:1818482/93‑1 (MQ=255)
                                                  gcacaAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTtcatca                                       >  1:5428212/1‑93 (MQ=255)
                                                    acaAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGa                                     >  2:4315502/1‑93 (MQ=255)
                                                    acaAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGa                                     >  2:3197143/1‑93 (MQ=255)
                                                     caAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGaa                                    >  1:4840650/1‑93 (MQ=255)
                                                     caAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGaa                                    >  2:4131149/1‑93 (MQ=255)
                                                     caAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGaa                                    <  2:5046218/93‑1 (MQ=255)
                                                      aaGATCTGCCGCAGCTATTGTTATGGCAGCAGTTGCGGGTGAGCAAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAAc                                   <  1:6486652/93‑1 (MQ=255)
                                                      aaGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAAc                                   >  1:6584048/1‑93 (MQ=255)
                                                      aaGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAAc                                   <  2:6545235/93‑1 (MQ=255)
                                                      aaGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAAc                                   >  1:4432840/1‑93 (MQ=255)
                                                      aaGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAAc                                   <  2:4257068/93‑1 (MQ=255)
                                                        gATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATGCCTGTAATCAATTTTCGCTTTCATCAGAACGt                                 <  2:300018/93‑1 (MQ=255)
                                                        gATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGt                                 <  2:548688/93‑1 (MQ=255)
                                                        gATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGt                                 <  1:3962189/93‑1 (MQ=255)
                                                         aTCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTc                                <  2:5624494/93‑1 (MQ=255)
                                                          tCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTcc                               <  1:2797294/93‑1 (MQ=255)
                                                            tGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCga                             <  1:2637234/93‑1 (MQ=255)
                                                            tGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCga                             >  2:2171372/1‑93 (MQ=255)
                                                             gCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCgat                            <  2:4222343/93‑1 (MQ=255)
                                                                gCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCgatgat                         <  2:4825963/93‑1 (MQ=255)
                                                                 cAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGAtt                        <  2:611804/93‑1 (MQ=255)
                                                                 cAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGAtt                        >  2:5133275/1‑93 (MQ=255)
                                                                   gCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTgg                      >  2:6630455/1‑93 (MQ=255)
                                                                   gCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTgg                      <  1:5123334/93‑1 (MQ=255)
                                                                    ccATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTgct                     <  1:3468772/93‑3 (MQ=255)
                                                                    ccATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGt                     <  2:5940564/93‑1 (MQ=255)
                                                                    ccATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGt                     <  1:1172988/93‑1 (MQ=255)
                                                                       ttGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGtttt                  <  2:1946338/93‑1 (MQ=255)
                                                                       ttGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGtttt                  >  2:200560/1‑93 (MQ=255)
                                                                          ttATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGc               >  2:4028974/1‑93 (MQ=255)
                                                                             tGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCg            <  1:1301417/93‑1 (MQ=255)
                                                                                     gTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCGGGTGCGGa    <  1:1408265/93‑1 (MQ=255)
                                                                                     gTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCGGGTGCGGa    <  2:6187112/93‑1 (MQ=255)
                                                                                     gTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCGGGTGCGGa    >  1:2763218/1‑93 (MQ=255)
                                                                                      ttGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCGGGTGCGGAg   >  1:6022277/1‑93 (MQ=255)
                                                                                      ttGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCGGGTGCGGAg   >  2:2462971/1‑93 (MQ=255)
                                                                                          |                                                                                         
TGGTGGAACGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCGGGTGCGGAGT  >  NC_000913/3792947‑3793126

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