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

Only 100 of 289 total aligned reads displayed.

TGGTGGAACGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCGGGTGCGGAGT  >  NC_000913/3792947‑3793126
                                                                                          |                                                                                         
tggtggAACGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCggg                                                                                         >  2:299491/1‑93 (MQ=255)
tggtggAACGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCggg                                                                                         >  2:557589/1‑93 (MQ=255)
 ggtggAACGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGt                                                                                        <  1:1419043/93‑1 (MQ=255)
    ggAACGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAg                                                                                     >  1:3356455/1‑93 (MQ=255)
    ggAACGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAg                                                                                     >  1:2869123/1‑93 (MQ=255)
    ggAACGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAg                                                                                     >  1:2165756/1‑93 (MQ=255)
     gAACGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGc                                                                                    <  2:2931949/93‑1 (MQ=255)
      aaCGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGca                                                                                   >  2:1983172/1‑92 (MQ=255)
        cGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGaa                                                                                 <  1:3719086/93‑1 (MQ=255)
         gCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAg                                                                                <  1:2280978/93‑1 (MQ=255)
          cTATATAGCGCTGGCCTATGACAAAGGCATTAAGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAgg                                                                               <  2:463601/93‑1 (MQ=255)
          cTATATAGCGCTGGCCTATGACAAAGCCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAgg                                                                               <  2:1907763/93‑1 (MQ=255)
           tatataGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGt                                                                              <  2:2849347/93‑1 (MQ=255)
            atataGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTg                                                                             >  2:2233129/1‑93 (MQ=255)
            atataGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTg                                                                             >  2:23536/1‑93 (MQ=255)
             tataGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGa                                                                            <  1:1099907/93‑1 (MQ=255)
             tataGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGa                                                                            <  1:1663350/93‑1 (MQ=255)
              ataGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGaa                                                                           >  1:1164183/1‑93 (MQ=255)
              ataGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGaa                                                                           >  2:3082930/1‑93 (MQ=255)
              ataGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGaa                                                                           <  1:3277937/93‑1 (MQ=255)
               taGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGaaa                                                                          >  2:3426504/1‑93 (MQ=255)
               taGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGaaa                                                                          <  1:1711055/93‑1 (MQ=255)
                aGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGaaaa                                                                         <  1:19024/93‑1 (MQ=255)
                aGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGaaaa                                                                         <  1:593840/93‑1 (MQ=255)
                aGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGaaaa                                                                         <  2:3847004/93‑1 (MQ=255)
                 gcgcTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGaaaaa                                                                        >  1:3237142/1‑93 (MQ=255)
                 gcgcTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGaaaaa                                                                        <  1:4508422/93‑1 (MQ=255)
                  cgcTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAAt                                                                       <  1:150187/93‑1 (MQ=255)
                   gcTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATc                                                                      <  2:869999/93‑1 (MQ=255)
                    cTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCa                                                                     <  1:2002110/93‑1 (MQ=255)
                     tGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGAGAGCGAAGGTGAAAaataat                                                                    >  1:935723/1‑93 (MQ=255)
                       gCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCatat                                                                  <  2:2675069/93‑1 (MQ=255)
                       gCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCatat                                                                  <  1:741130/93‑1 (MQ=255)
                        ccTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCAtata                                                                 <  2:2586665/93‑1 (MQ=255)
                        ccTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCAtata                                                                 <  2:460898/93‑1 (MQ=255)
                        ccTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAAACAtata                                                                 >  2:2462103/1‑93 (MQ=255)
                          tATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATAcc                                                               <  2:3468936/93‑1 (MQ=255)
                          tATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATAcc                                                               >  2:307142/1‑93 (MQ=255)
                          tATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATAcc                                                               <  2:3604139/93‑1 (MQ=255)
                          tATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATAcc                                                               <  1:3641652/93‑1 (MQ=255)
                          tATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGCTTCGGGTGAGCGAAGGTGAAAAGTCATATAcc                                                               >  2:1190082/1‑93 (MQ=255)
                            tGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTg                                                             >  1:3845942/1‑93 (MQ=255)
                            tGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTg                                                             <  2:3571143/93‑1 (MQ=255)
                              aCAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTa                                                           >  1:3933896/1‑93 (MQ=255)
                              aCAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTa                                                           >  1:2725489/1‑93 (MQ=255)
                               cAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTaa                                                          >  2:3896089/1‑93 (MQ=255)
                               cAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTaa                                                          <  2:1630620/93‑1 (MQ=255)
                                 aaGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGGAAAATCATATACCTGTAAtt                                                        >  2:1720024/1‑92 (MQ=255)
                                 aaGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATc                                                        <  2:4135590/93‑1 (MQ=255)
                                  aGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCa                                                       >  2:3948070/1‑93 (MQ=255)
                                   ggCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCaa                                                      <  2:3692941/93‑1 (MQ=255)
                                   ggCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCaa                                                      >  1:301612/1‑93 (MQ=255)
                                    gCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAAt                                                     <  2:3813323/93‑1 (MQ=255)
                                     cATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAAtt                                                    <  2:805913/93‑1 (MQ=255)
                                      aTTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAAttt                                                   >  2:1119787/1‑93 (MQ=255)
                                      aTTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAAttt                                                   <  2:949392/93‑1 (MQ=255)
                                        tATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCCGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTc                                                 >  1:1194248/1‑93 (MQ=255)
                                          tGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGc                                               <  1:4391772/93‑1 (MQ=255)
                                            cgcACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCtt                                             <  2:638938/93‑1 (MQ=255)
                                            cgcACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCtt                                             <  2:1754572/93‑1 (MQ=255)
                                            cgcACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCtt                                             <  2:1845614/93‑1 (MQ=255)
                                            cgcACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCtt                                             <  1:1853078/93‑1 (MQ=255)
                                             gcacagcacaAGATCTGCCGCAGCCATTGTTATGGCCGCCGTTGCGGGTGAGCGAAGGTGGAAACTCATATACCTGTAATCAATTTTCGCttt                                            >  2:1006337/1‑93 (MQ=255)
                                             gcacagcacaAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCttt                                            >  1:2063159/1‑93 (MQ=255)
                                              cacagcacaAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTtc                                           >  1:352799/1‑93 (MQ=255)
                                              cacagcacaAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTtc                                           <  1:498566/93‑1 (MQ=255)
                                               acagcacaAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTtca                                          <  1:1278051/93‑1 (MQ=255)
                                                    acaAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGa                                     >  2:1624696/1‑93 (MQ=255)
                                                     caAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGaa                                    <  1:2398630/93‑1 (MQ=255)
                                                       aGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACg                                  <  1:4164608/93‑1 (MQ=255)
                                                       aGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACg                                  >  1:2652428/1‑93 (MQ=255)
                                                        aaccTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGt                                 <  2:720156/90‑1 (MQ=255)
                                                          tCTGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTcc                               >  1:491648/1‑93 (MQ=255)
                                                            tGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCga                             <  2:2756015/93‑1 (MQ=255)
                                                            tGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCga                             <  1:2496505/93‑1 (MQ=255)
                                                            tGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCga                             <  2:3177910/93‑1 (MQ=255)
                                                            tGCCGCAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCga                             <  1:3991255/93‑1 (MQ=255)
                                                                 cAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGAtt                        <  2:4477850/93‑1 (MQ=255)
                                                                 cAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGAtt                        >  1:2705422/1‑93 (MQ=255)
                                                                 cAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGAtt                        <  1:3411330/93‑1 (MQ=255)
                                                                 cAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGAtt                        <  1:1047522/93‑1 (MQ=255)
                                                                 cAGCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGAtt                        >  1:4272323/1‑93 (MQ=255)
                                                                   gCCATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTgt                      >  1:1794924/1‑92 (MQ=255)
                                                                    ccATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGt                     <  1:1037513/93‑1 (MQ=255)
                                                                     cATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGtt                    <  2:2166566/93‑1 (MQ=255)
                                                                     cATTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGtt                    <  2:3132273/93‑1 (MQ=255)
                                                                      aTTGTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGttt                   >  1:2973440/1‑93 (MQ=255)
                                                                         gTTATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTg                <  2:4424625/93‑1 (MQ=255)
                                                                          ttATGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGc               >  1:3850266/1‑93 (MQ=255)
                                                                             tGGCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCg            >  2:3237445/1‑93 (MQ=255)
                                                                               gCCGCAGTTGCGGGTGATCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCggg          <  1:2932290/93‑1 (MQ=255)
                                                                               gCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCggg          >  1:256847/1‑93 (MQ=255)
                                                                               gCCGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCggg          <  2:136793/93‑1 (MQ=255)
                                                                                ccGCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCGGGt         >  2:3396826/1‑93 (MQ=255)
                                                                                  gCAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCGGGTGc       >  2:4305380/1‑93 (MQ=255)
                                                                                   cAGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCGGGTGCg      >  1:1930369/1‑93 (MQ=255)
                                                                                    aGTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCGGGTGCgg     >  2:3642657/1‑93 (MQ=255)
                                                                                     gTTGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCGGGTGCGGa    >  1:4461794/1‑93 (MQ=255)
                                                                                     gTTGCGGGTGAGCGAAGGCGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTCGCCGGGTGCGGa    >  2:2607176/1‑93 (MQ=255)
                                                                                      ttGCGGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCGGGTGCGGAg   <  1:1333136/93‑1 (MQ=255)
                                                                                          |                                                                                         
TGGTGGAACGCTATATAGCGCTGGCCTATGACAAAGGCATTATGCGCACAGCACAAGATCTGCCGCAGCCATTGTTATGGCCGCAGTTGCAGGTGAGCGAAGGTGAAAAATCATATACCTGTAATCAATTTTCGCTTTCATCAGAACGTCCGATGATTGGTTTTTGCCCGGGTGCGGAGT  >  NC_000913/3792947‑3793126

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