Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,139,436 G→A intergenic (‑1/‑1) ptsA ← / → yijI PEP‑protein phosphotransferase system enzyme I/IIA/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,139,4360GA85.6% 913.2 / 95.8 319intergenic (‑1/‑1)ptsA/yijIPEP‑protein phosphotransferase system enzyme I/IIA/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (26/20);  new base A (146/127);  total (172/147)
Fisher's exact test for biased strand distribution p-value = 7.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 324 total aligned reads displayed.

AACGCCTTTGGCTGCGGGAAGATTACCCAGCGCATTGAGATCTAAAGAAGAGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACGGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCAGTGGTTCCAGTTCGTCAGATTTAACTTCCGCCAGCGGCGCATC  >  NC_000913/4139346‑4139528
                                                                                          |                                                                                             
aaCGCCTTTGGCTGCGGGAAGATTACCCAGCGCATTGAGATCTAAAGAAGAGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACGGt                                                                                             <  1:3638978/93‑1 (MQ=255)
  cGCCTTTGGCTGCGGGAAGATTACCCAGCGCATTGAGATCTAAAGAAGAGTTCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACAGTGc                                                                                           >  2:5245299/1‑93 (MQ=255)
  cGCCTTTGGCTGCGGGAAGATTACCCAGCGCATTGAGATCTAAAGAAGAGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACAGTGc                                                                                           <  1:4079540/93‑1 (MQ=255)
   gCCTTTGGCTGCGGGAAGATTACCCAGCGCATTGAGATCTAAAGAAGAGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACAGTGCg                                                                                          <  1:1124235/93‑1 (MQ=255)
   gCCTTTGGCTGCGGGAAGATTACCCAGCGCATTGAGATCTAAAGAAGAGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACAGTGCg                                                                                          >  1:870862/1‑93 (MQ=255)
    ccTTTGGCTGCGGGAAGATTACCCAGCGCATTGAGATCTAAAGAAGAGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACAGTGCgg                                                                                         >  1:5357879/1‑93 (MQ=255)
    ccTTTGGCTGCGGGAAGATTACCCAGCGCATTGAGATCTAAAGAAGAGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACAGTGCgg                                                                                         <  2:1735576/93‑1 (MQ=255)
     cTTTGGCTGCGGGAAGATTACCCAGCGCATTGAGATCTAAAGAAGAGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACGGTGCggg                                                                                        >  2:1344251/1‑93 (MQ=255)
      tttGGCTGCGGGAAGATTACCCAGCGCATTGAGATCTAAAGAAGAGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACAGTGCGGgc                                                                                       >  2:5674864/1‑93 (MQ=255)
       ttGGCTGCGGGAAGATTACCCAGCGCATTGAGATCTAAAGAAGAGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACAGTGCGGgcg                                                                                      <  1:489660/93‑1 (MQ=255)
            tGCGGGAAGATTACCCAGCGCATTGAGATCTAAAGAAGAGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACAGTGCGGGCGCGGAt                                                                                >  2:5473193/1‑93 (MQ=255)
            tGCGGGAAGATTACCCAGCGCATTGAGATCTAAAGAAGAGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACAGTGCGGGCGCGGAt                                                                                >  1:4520575/1‑93 (MQ=255)
             gCGGGAAGATTACCCAGCGCATTGAGATCTAAAGAAGAGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACAGTGCGGGCGCGGATa                                                                               >  1:4739627/1‑93 (MQ=255)
                gggaGATTACCCAGCGCATTGAGATCTAAAGAAGAGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACAGTGCGGGCGCGGATAAtt                                                                            <  2:4339916/90‑1 (MQ=255)
                   aGATTACCCAGCGCATTGAGATCTAAAGAAGAGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGc                                                                         >  1:2672464/1‑93 (MQ=255)
                   aGATTACCCAGCGCATTGAGATCTAAAGAAGAGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGc                                                                         <  2:6794788/93‑1 (MQ=255)
                    gATTACCCAGCGCATTGAGATCTAAAGAAGAGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCg                                                                        <  2:2645760/93‑1 (MQ=255)
                     aTTACCCAGCGCATTGAGATCTAAAGAAGAGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACGGTGCGGGCGCGGATAATTTGCgg                                                                       >  2:2766010/1‑93 (MQ=255)
                      ttACCCAGCGCATTGAGATCTAAAGAAGAGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGa                                                                      <  1:7438924/93‑1 (MQ=255)
                      ttACCCAGCGCATTGAGATCTAAAGAAGAGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGa                                                                      >  1:6395939/1‑93 (MQ=255)
                      ttACCCAGCGCATTGAGATCTAAAGAAGAGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGa                                                                      <  1:4807730/93‑1 (MQ=255)
                          ccAGCGCTTTGAGATCTAAAGAAGAGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACGGTGCGGGCGCGGATAATTTGCGGATTCa                                                                  >  2:2365387/1‑93 (MQ=255)
                            aGCGCATTGAGATCTAAAGAAGAGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGa                                                                <  2:410436/93‑1 (MQ=255)
                            aGCGCATTGAGATCTAAAGAAGAGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGa                                                                >  1:7182271/1‑93 (MQ=255)
                             gcgcATTGAGATCTAAAGAAGAGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGAt                                                               >  2:6915104/1‑93 (MQ=255)
                             gcgcATTGAGATCTAAAGAAGAGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGAt                                                               <  1:1889095/93‑1 (MQ=255)
                              cgcATTGAGATCTAAAGAAGAGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACGGTGCGGGCGCGGATAATTTGCGGATTCAGAtt                                                              <  1:4595925/93‑1 (MQ=255)
                               gcATTGAGATCTAAAGAAGAGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTg                                                             <  2:228810/93‑1 (MQ=255)
                                cATTGAGATCTAAAGAAGAGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACGGTGCGGGCGCGGATAATTTGCGGATTCAGATTgg                                                            <  1:1046609/93‑1 (MQ=255)
                                cATTGAGATCTAAAGAAGAGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTgg                                                            >  1:4284746/1‑93 (MQ=255)
                                 aTTGAGATCTAAAGAAGAGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGt                                                           <  2:4836484/93‑1 (MQ=255)
                                  ttGAGATCTAAAGAAGAGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTc                                                          >  1:5088225/1‑93 (MQ=255)
                                  ttGAGATCTAAAGAAGAGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTc                                                          >  1:6100827/1‑93 (MQ=255)
                                   tGAGATCTAAAGAAGAGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCa                                                         <  2:3420243/93‑1 (MQ=255)
                                      gaTCTAAAGAAGAGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACGGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTg                                                      >  1:1711935/1‑93 (MQ=255)
                                       aTCTAAAGAAGAGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGa                                                     <  2:6340673/93‑1 (MQ=255)
                                       aTCTAAAGAAGAGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGa                                                     >  2:3150086/1‑93 (MQ=255)
                                       aTCTAAAGAAGAGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGa                                                     <  2:3797543/93‑1 (MQ=255)
                                        tCTAAAGAAGAGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGaa                                                    <  2:2093198/93‑1 (MQ=255)
                                        tCTAAAGAAGAGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGaa                                                    >  2:1559616/1‑93 (MQ=255)
                                          tAAAGAAGAGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACGGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAAc                                                  >  2:1048399/1‑93 (MQ=255)
                                          tAAAGAAGAGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACGGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAAc                                                  <  2:7328162/93‑1 (MQ=255)
                                          tAAAGAAGAGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAAc                                                  <  1:3745926/93‑1 (MQ=255)
                                            aagaagAGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACGGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCg                                                >  2:1197475/1‑93 (MQ=255)
                                              gaagaGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACGGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGc                                              >  1:3291558/1‑93 (MQ=255)
                                              gaagaGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGc                                              >  2:1868775/1‑93 (MQ=255)
                                              gaagaGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGc                                              >  2:4451109/1‑93 (MQ=255)
                                              gaagaGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGc                                              >  2:6295808/1‑93 (MQ=255)
                                               aagaGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCa                                             <  1:5562239/93‑1 (MQ=255)
                                               aagaGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCa                                             <  2:5842460/93‑1 (MQ=255)
                                               aagaGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGACa                                             <  2:4732854/93‑1 (MQ=255)
                                                agagATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACGGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCAg                                            >  1:3112842/1‑93 (MQ=255)
                                                 gagaTCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCAGt                                           <  1:1004703/93‑1 (MQ=255)
                                                 gagaTCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCAGt                                           <  2:2921262/93‑1 (MQ=255)
                                                  agaTCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCAGTg                                          <  2:4040834/93‑1 (MQ=255)
                                                   gaTCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCAGTgg                                         >  2:2965727/1‑93 (MQ=255)
                                                    aTCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACGGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCAGTGGt                                        <  2:5293527/93‑1 (MQ=255)
                                                    aTCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCAGTGGt                                        >  2:2181371/1‑93 (MQ=255)
                                                      cGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCAGTGGTTc                                      <  1:1335807/93‑1 (MQ=255)
                                                       ggCGTCAGAATGCCGCCTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCAGTGGTTcc                                     <  2:7474617/93‑1 (MQ=255)
                                                        gCGTCAGAATGCCGCCTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCAGTGGTTCCa                                    <  1:6712710/93‑1 (MQ=255)
                                                           tCAGAATGCCGCCTGCACTACCGCTGCACACGGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCAGTGGTTCCAGtt                                 >  2:3100858/1‑93 (MQ=255)
                                                           tCAGAATGCCGCCTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCAGTGGTTCCAGtt                                 >  2:1448447/1‑93 (MQ=255)
                                                           tCAGAATGCCGCCTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCAGTGGTTCCAGtt                                 >  1:1603082/1‑93 (MQ=255)
                                                           tCAGAATGCCGCCTCCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCAGTGGTTCCAGtt                                 >  2:3350247/1‑93 (MQ=255)
                                                            cAGAATGCCGCCTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCAGTGGTTCCAGTTc                                <  2:6892837/93‑1 (MQ=255)
                                                            cAGAATGCCGCCTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCAGTGGTTCCAGTTc                                <  1:5221240/93‑1 (MQ=255)
                                                              gAATGCCGCCTGCACTACCGCTGCACACGGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCAGTGGTTCCAGTTCGt                              <  2:6024184/93‑1 (MQ=255)
                                                              gAATACCGCCTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCAGTGGTTCCAGTTCGt                              >  1:5998706/1‑93 (MQ=255)
                                                                aggccgccTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCAGTGGTTCCAGTTCGTCa                            <  2:6760188/91‑1 (MQ=255)
                                                                aTGCCGCCTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCAGTGGTTCCAGTTCGTCa                            <  2:7038083/93‑1 (MQ=255)
                                                                 tGCCGCCTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTGAGATTGGTCAGTGAAACCGGCAGTGGTTCCAGTTCGTCAg                           >  1:3494457/1‑93 (MQ=255)
                                                                 tGCCGCCTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGAGTCAGATTGGTCAGTGAAACCGGCAGTGGTTCCAGTTCGTCAg                           >  1:7373980/1‑93 (MQ=255)
                                                                   ccgccTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCAGTGGTTCCAGTTCGTCAGAt                         >  2:799098/1‑93 (MQ=255)
                                                                    cgccTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCAGTGGTTCCAGTTCGTCAGAtt                        >  1:4996998/1‑93 (MQ=255)
                                                                    cgccTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCAGTGGTTCCAGTTCGTCAGAtt                        <  2:4611893/93‑1 (MQ=255)
                                                                      ccTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCAGTGGTTCCAGTTCGTCAGATTTa                      >  1:2887228/1‑93 (MQ=255)
                                                                       cTGCACTACCGCTGCACACGGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCAGTGGTTCCAGTTCGTCAGATTTaa                     >  2:2867171/1‑93 (MQ=255)
                                                                       cTGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCAGTGGTTCCAGTTCGTCAGATTTaa                     >  1:577761/1‑93 (MQ=255)
                                                                        tGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCAGTGGTTCCAGTTCGTCAGATTTAAc                    <  1:6509459/93‑1 (MQ=255)
                                                                        tGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCAGTGGTTCCAGTTCGTCAGATTTAAc                    <  2:3227306/93‑1 (MQ=255)
                                                                        tGCACTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCAGTGGTTCCAGTTCGTCAGATTTAAc                    >  1:6964760/1‑93 (MQ=255)
                                                                           aCTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCAGTGGTTCCAGTTCGTCAGATTTAACTTc                 <  2:1083772/93‑1 (MQ=255)
                                                                           aCTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCAGTGGTTCCAGTTCGTCAGATTTAACTTc                 >  1:2080108/1‑93 (MQ=255)
                                                                           aCTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCAGTGGTTCCAGTTCGTCAGATTTAACTTc                 >  1:3851147/1‑93 (MQ=255)
                                                                           aCTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCAGTGGTTCCAGTTCGTCAGATTTAACTTc                 >  2:3496614/1‑93 (MQ=255)
                                                                           aCTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCAGTGGTTCCAGTTCGTCAGATTTAACTTc                 >  1:3394364/1‑93 (MQ=255)
                                                                           aCTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCAGTGGTTCCAGTTCCTCAGATTTAACTTc                 <  2:6445461/93‑1 (MQ=255)
                                                                            cTACCGCTGCACACGGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCAGTGGTTCCAGTTCGTCAGATTTAACTTcc                >  2:5596351/1‑93 (MQ=255)
                                                                            cTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCAGTGGTTCCAGTTCGTCAGATTTAACTTcc                >  1:5632750/1‑93 (MQ=255)
                                                                            cTACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCAGTGGTTCCAGTTCGTCAGATTTAACTTcc                >  2:4500335/1‑93 (MQ=255)
                                                                             tACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCAGTGGTTCCAGTTCGTCAGATTTAACTTCCg               >  2:3845522/1‑93 (MQ=255)
                                                                             tACCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCAGTGGTTCCAGTTCGTCAGATTTAACTTCCg               >  1:5886468/1‑93 (MQ=255)
                                                                              aCCGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCAGTGGTTCCAGTTCGTCAGATTTAACTTCCGc              >  1:2799458/1‑93 (MQ=255)
                                                                                cGCTGCACACGGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCAGTGGTTCCAGTTCGTCAGATTTAACTTCCGCCa            <  2:2598132/93‑1 (MQ=255)
                                                                                cGCTGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCAGTGGTTCCAGTTCGTCAGATTTAACTTCCGCCa            <  1:4031102/93‑1 (MQ=255)
                                                                                  cTGCACACAGTGCGGGCGCGGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCAGTGGTTCCAGTTCGTCAGATTTAACTTCCGCCAg           >  1:2344559/1‑93 (MQ=255)
                                                                                   tGCACACGGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCAGTGGTTCCAGTTCGTCAGATTTAACTTCCGCCAgcg         >  2:503742/1‑93 (MQ=255)
                                                                                   tGCACACAGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCAGTGGTTCCAGTTCGTCAGATTTAACTTCCGCCAgcg         >  1:5458865/1‑93 (MQ=255)
                                                                                     cacacaGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCAGTGGTTCCAGTTCGTCAGATTTAACTTCCGCCAgcggc       <  2:1634304/93‑1 (MQ=255)
                                                                                          |                                                                                             
AACGCCTTTGGCTGCGGGAAGATTACCCAGCGCATTGAGATCTAAAGAAGAGATCGGCGTCAGAATGCCGCCTGCACTACCGCTGCACACGGTGCGGGCGCGGATAATTTGCGGATTCAGATTGGTCAGTGAAACCGGCAGTGGTTCCAGTTCGTCAGATTTAACTTCCGCCAGCGGCGCATC  >  NC_000913/4139346‑4139528

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