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,0370AG98.0% 632.1 / 0.7 202Q161R (CAG→CGG) rfaFADP‑heptose‑‑lps heptosyltransferase II
Reads supporting (aligned to +/- strand):  ref base A (1/3);  new base G (94/104);  total (95/107)
Fisher's exact test for biased strand distribution p-value = 6.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.91e-01

Only 100 of 213 total aligned reads displayed.

TGGTGGAACGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCGGGTGCGGAG  >  NC_000913/3792947‑3793125
                                                                                          |                                                                                        
tggtggAACGCTATATAGCGCTGTACTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCggg                                                                                        >  1:2556803/1‑93 (MQ=255)
  gtggAACGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTg                                                                                      >  2:1179117/1‑93 (MQ=255)
  gtggAACGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTg                                                                                      <  1:2023127/93‑1 (MQ=255)
  gtggAACGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTg                                                                                      <  2:3171621/93‑1 (MQ=255)
    ggAACGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAg                                                                                    >  1:3265346/1‑93 (MQ=255)
     gAACGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGc                                                                                   <  1:1001010/93‑1 (MQ=255)
     gAACGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGc                                                                                   <  2:3926604/93‑1 (MQ=255)
      aaCGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGCTGAGCg                                                                                  <  2:2233681/93‑1 (MQ=255)
       ttgCTATATAGCGCTGGTCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTCGCGGGTTAGCGa                                                                                 >  2:4221423/3‑93 (MQ=255)
        cGCTATATAGCGCTGGCCTATGACAAAGGCACTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGaa                                                                                <  2:434926/93‑1 (MQ=255)
         gCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAg                                                                               >  1:3106823/1‑93 (MQ=255)
          cTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAgg                                                                              <  1:3894553/93‑1 (MQ=255)
           tatataGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGt                                                                             <  2:3324464/93‑1 (MQ=255)
            atataGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTg                                                                            >  1:1541772/1‑93 (MQ=255)
            atataGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTg                                                                            <  1:3547149/93‑1 (MQ=255)
             tataGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGa                                                                           >  1:405187/1‑93 (MQ=255)
             tataGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGa                                                                           >  2:4555278/1‑93 (MQ=255)
                aGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAATTGCGGGTGAGCGAAGGTGaaaa                                                                        >  2:3387122/1‑93 (MQ=255)
                 gcgcTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGaaaaa                                                                       <  1:1086213/93‑1 (MQ=255)
                 gcgcTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGaaaaa                                                                       >  2:4010735/1‑93 (MQ=255)
                   gcTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATc                                                                     <  2:1817303/93‑1 (MQ=255)
                    cTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCa                                                                    <  2:1688155/93‑1 (MQ=255)
                     tGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCat                                                                   <  2:361886/93‑1 (MQ=255)
                      ggCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCata                                                                  >  2:1754105/1‑93 (MQ=255)
                      ggCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCata                                                                  <  2:3015201/93‑1 (MQ=255)
                       gCCTATGACAAGGGCATTATGCGCACAGCACAGGATCTGCGGCAGCCATTGTTATGGCCGCAGTTGAGGGTGAGCGAAGGTGAAAAATCatat                                                                 <  2:2656182/93‑1 (MQ=255)
                       gCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCatat                                                                 <  1:3832455/93‑1 (MQ=255)
                       gCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCatat                                                                 <  2:122937/93‑1 (MQ=255)
                        ccTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCAtata                                                                <  2:517394/93‑1 (MQ=255)
                          tATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATAcc                                                              >  1:2489517/1‑93 (MQ=255)
                          tATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATAcc                                                              <  1:452481/93‑1 (MQ=255)
                           aTGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCt                                                             <  1:3391841/93‑1 (MQ=255)
                           aTGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCt                                                             <  2:874439/93‑1 (MQ=255)
                            tGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTg                                                            <  2:4452757/93‑1 (MQ=255)
                             gACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGt                                                           >  1:2743257/1‑93 (MQ=255)
                              aCAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTa                                                          >  2:3059377/1‑93 (MQ=255)
                                 aaGCCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATc                                                       <  1:624806/93‑1 (MQ=255)
                                   ggCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCaa                                                     >  2:530558/1‑93 (MQ=255)
                                   ggCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCaa                                                     >  1:231459/1‑93 (MQ=255)
                                   ggCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCaa                                                     >  2:3534559/1‑93 (MQ=255)
                                   ggCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCaa                                                     <  1:3688594/93‑1 (MQ=255)
                                   ggCAGTATGCGCACAGGACAGGATCTGCCGCAGCCATGGTTAGGGCCGAAGTTGGGGGTGAGCGAAGGTGAAAAATCATATGCCTGTAATCaa                                                     <  1:4278390/93‑1 (MQ=255)
                                      aTTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAAttt                                                  >  1:1455168/1‑93 (MQ=255)
                                      aTTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAAttt                                                  >  2:1779603/1‑93 (MQ=255)
                                      aTTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAAttt                                                  <  1:3998332/93‑1 (MQ=255)
                                       ttATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAAtttt                                                 <  1:4700368/93‑1 (MQ=255)
                                       ttATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAAtttt                                                 <  1:596992/93‑1 (MQ=255)
                                       ttAGGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAAtttt                                                 <  2:2425087/93‑1 (MQ=255)
                                        tATGCGCACAGCCCACGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTc                                                <  2:4212679/93‑1 (MQ=255)
                                        tATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTc                                                <  1:4189158/93‑1 (MQ=255)
                                        tATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTc                                                <  1:1983991/93‑1 (MQ=255)
                                          tGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGc                                              <  1:1363537/93‑1 (MQ=255)
                                           gcgcACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCt                                             <  2:3652233/93‑1 (MQ=255)
                                            cgcACAGCACAAGATCTGCCGCATCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCtt                                            <  2:2581433/93‑1 (MQ=255)
                                            cgcACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCtt                                            <  2:2033012/93‑1 (MQ=255)
                                             gcacagcacaAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCttt                                           >  1:2924731/1‑93 (MQ=255)
                                             gcacagcacaAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCttt                                           >  2:1125112/1‑93 (MQ=255)
                                             gcacagcacaAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCttt                                           <  2:1654135/93‑1 (MQ=255)
                                               acagcacaAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTtca                                         <  2:121361/93‑1 (MQ=255)
                                                  gcacaAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTtcatca                                      <  2:4075624/93‑1 (MQ=255)
                                                   cacaAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAg                                     >  1:2652080/1‑93 (MQ=255)
                                                    acaAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGa                                    >  1:1744421/1‑93 (MQ=255)
                                                     caAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGaa                                   >  1:3422714/1‑93 (MQ=255)
                                                     caAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGaa                                   <  1:1099373/93‑1 (MQ=255)
                                                     aaaGACCGGCCGCAGCCATTGTTATAGCCGCAGTTGCGGGGTAGTGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGaa                                   <  1:4367888/92‑1 (MQ=255)
                                                      aaGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAAc                                  >  2:3230312/1‑93 (MQ=255)
                                                      aaGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAAc                                  >  2:1566508/1‑93 (MQ=255)
                                                       aGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACg                                 >  2:342853/1‑93 (MQ=255)
                                                       aGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACg                                 >  2:2472865/1‑93 (MQ=255)
                                                       aGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCCTTCATCAGAACg                                 >  2:1974706/1‑93 (MQ=255)
                                                        gATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGGAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGt                                <  2:3842289/93‑1 (MQ=255)
                                                        gATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGt                                <  1:3305612/93‑1 (MQ=255)
                                                        gATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGt                                >  1:2285042/1‑93 (MQ=255)
                                                         aTCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTc                               >  2:1377599/1‑93 (MQ=255)
                                                           cTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCg                             <  1:2855224/93‑1 (MQ=255)
                                                             gCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCgat                           <  2:709675/93‑1 (MQ=255)
                                                              ccGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCgatg                          <  2:2118549/93‑1 (MQ=255)
                                                              ccGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCgatg                          <  1:2727537/93‑1 (MQ=255)
                                                                 cAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGAtt                       <  2:2816300/93‑1 (MQ=255)
                                                                   gCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTgt                     >  1:1172020/1‑92 (MQ=255)
                                                                   gCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTgg                     >  1:439724/1‑93 (MQ=255)
                                                                    ccATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGt                    >  2:4427096/1‑93 (MQ=255)
                                                                     cATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGtt                   >  1:3760726/1‑93 (MQ=255)
                                                                     cATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGtt                   >  1:1284723/1‑93 (MQ=255)
                                                                     cATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGtt                   <  1:2581976/93‑1 (MQ=255)
                                                                     cATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGtt                   >  2:3300334/1‑93 (MQ=255)
                                                                      aTTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGttt                  <  2:1384782/93‑1 (MQ=255)
                                                                      aTTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGttt                  <  2:4066584/93‑1 (MQ=255)
                                                                        tGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGttttt                <  2:3487886/93‑1 (MQ=255)
                                                                          ttATGGGCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGc              <  1:1163695/93‑1 (MQ=255)
                                                                           tATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGcc             <  2:1340513/93‑1 (MQ=255)
                                                                            aTGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGccc            >  1:3054247/1‑93 (MQ=255)
                                                                             tGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCg           <  1:926882/93‑1 (MQ=255)
                                                                               gCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCggg         <  1:3175212/93‑1 (MQ=255)
                                                                                ccGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCGGGt        <  1:3219841/93‑1 (MQ=255)
                                                                                ccGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCGGGt        <  1:1855749/93‑1 (MQ=255)
                                                                                 cGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCGGGTg       >  2:232911/1‑93 (MQ=255)
                                                                                     gTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCGGGTGCGGa   <  2:3434947/93‑1 (MQ=255)
                                                                                     gTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCGGGTGCGGa   >  2:486653/1‑93 (MQ=255)
                                                                                      ttGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCGGGTGCGGAg  >  2:605033/1‑93 (MQ=255)
                                                                                          |                                                                                        
TGGTGGAACGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCGGGTGCGGAG  >  NC_000913/3792947‑3793125

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