Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,725,180 G→A I70I (ATC→ATT yiaB ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,725,1800GA99.5% 2791.9 / ‑3.0 762I70I (ATC→ATTyiaBhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (1/3);  new base A (334/424);  total (335/427)
Fisher's exact test for biased strand distribution p-value = 6.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

Only 100 of 769 total aligned reads displayed.

GGGATAGATCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAGATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAACACGCCAAG  >  NC_000913/3725090‑3725269
                                                                                          |                                                                                         
gggATAGATCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAat                                                                                         >  1:4219446/1‑93 (MQ=255)
 ggATAGATCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACATCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAata                                                                                        <  2:1857802/93‑1 (MQ=255)
 ggATAGATCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAata                                                                                        <  1:7747583/93‑1 (MQ=255)
  gATAGATCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAatat                                                                                       <  1:5193061/93‑1 (MQ=255)
   aTAGATCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATg                                                                                      <  1:3406835/93‑1 (MQ=255)
   aTAGATCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATg                                                                                      <  2:3263219/93‑1 (MQ=255)
    tAGATCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATgg                                                                                     <  2:8033371/93‑1 (MQ=255)
    tAGATCACCATTTCATAAGTATTGATAGGAGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATgg                                                                                     >  1:9375752/1‑93 (MQ=255)
      gATCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGt                                                                                   >  2:5927759/1‑93 (MQ=255)
       aTCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTa                                                                                  >  2:8330463/1‑93 (MQ=255)
       aTCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTa                                                                                  >  1:9060868/1‑93 (MQ=255)
        tCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTaa                                                                                 <  2:7654371/93‑1 (MQ=255)
         cACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTaaa                                                                                >  1:5928604/1‑93 (MQ=255)
         cACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTaaa                                                                                <  2:2155201/93‑1 (MQ=255)
          aCCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGATaaaa                                                                               >  2:5227124/1‑93 (MQ=255)
           ccATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTaaaaa                                                                              <  2:9300696/93‑1 (MQ=255)
            cATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTaaaaaa                                                                             <  1:6236755/93‑1 (MQ=255)
            cATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTaaaaaa                                                                             >  2:510252/1‑93 (MQ=255)
              tttCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAat                                                                           >  1:4615671/1‑93 (MQ=255)
              tttCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAat                                                                           >  2:7145273/1‑93 (MQ=255)
                  aTAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAatcatc                                                                       <  1:5363060/93‑1 (MQ=255)
                  aTAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAatcatc                                                                       <  2:8646718/93‑1 (MQ=255)
                   tAAGTATTGATAGGTTCGTTTAAGACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATcc                                                                      <  1:1521356/93‑1 (MQ=255)
                   tAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATcc                                                                      <  2:6257944/93‑1 (MQ=255)
                   tAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATcc                                                                      <  1:4941877/93‑1 (MQ=255)
                   tAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATcc                                                                      >  2:2981211/1‑93 (MQ=255)
                   tAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATcc                                                                      >  1:7354358/1‑93 (MQ=255)
                    aaGTATTGATAGGTGCGTTTANAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCa                                                                     <  1:1185574/93‑1 (MQ=255)
                     aGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAg                                                                    <  1:3628058/93‑1 (MQ=255)
                      gTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGt                                                                   >  2:6657309/1‑93 (MQ=255)
                      gTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGt                                                                   >  1:5077736/1‑93 (MQ=255)
                      gTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGt                                                                   >  2:4985629/1‑93 (MQ=255)
                       tATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGtt                                                                  <  1:8804454/93‑1 (MQ=255)
                       tATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGtt                                                                  <  1:9604568/93‑1 (MQ=255)
                       tATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGtt                                                                  >  2:8134235/1‑93 (MQ=255)
                        aTTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTg                                                                 >  2:899917/1‑93 (MQ=255)
                        aTTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTg                                                                 >  1:7983678/1‑93 (MQ=255)
                         ttGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGc                                                                >  2:250066/1‑93 (MQ=255)
                           gATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGccc                                                              <  2:9696204/93‑1 (MQ=255)
                           gATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGccc                                                              >  1:247582/1‑93 (MQ=255)
                             tAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGa                                                            <  1:10012906/93‑1 (MQ=255)
                               ggTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGt                                                          >  1:6533031/1‑93 (MQ=255)
                               ggTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGt                                                          >  2:451868/1‑93 (MQ=255)
                                gTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGtt                                                         >  2:4099569/1‑93 (MQ=255)
                                   cGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAGATATGGGTAAAAAAATCATCCAGTTGCCCGAGttttt                                                      <  1:8621380/93‑1 (MQ=255)
                                   cGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGttttt                                                      <  2:3175976/93‑1 (MQ=255)
                                   cGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGttttt                                                      >  2:6794593/1‑93 (MQ=255)
                                   cGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTCGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGttttt                                                      >  2:5340181/1‑93 (MQ=255)
                                     tttAAAACACCTATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCt                                                    <  2:8498856/93‑1 (MQ=255)
                                     tttAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCt                                                    >  2:6140875/1‑93 (MQ=255)
                                     tttAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCt                                                    >  2:9417616/1‑93 (MQ=255)
                                      ttAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTg                                                   >  1:4379862/1‑93 (MQ=255)
                                      ttAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTg                                                   <  2:6411582/93‑1 (MQ=255)
                                        aaaaCACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTgcg                                                 >  2:3862825/1‑93 (MQ=255)
                                           acacCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGa                                              <  2:5444929/93‑1 (MQ=255)
                                            cacCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGaa                                             <  2:1324764/93‑1 (MQ=255)
                                            cacCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGaa                                             <  1:6407995/93‑1 (MQ=255)
                                            cacCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGaa                                             <  1:6930325/93‑1 (MQ=255)
                                             acCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAg                                            >  1:2604874/1‑93 (MQ=255)
                                               cGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGat                                          <  1:7601116/93‑1 (MQ=255)
                                                gATAAACAGGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAATAATCATCCAGTTGCCCGAGTTTTTCAACGCGAAGata                                         <  2:8937134/93‑1 (MQ=255)
                                                gATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGata                                         <  2:686854/93‑1 (MQ=255)
                                                 aTAAACAAGAGACCAATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGatat                                        <  1:1940415/93‑1 (MQ=255)
                                                  tAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATg                                       <  2:957689/93‑1 (MQ=255)
                                                   aaaCAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGc                                      >  2:4259539/1‑93 (MQ=255)
                                                       aaGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAg                                  >  2:1648433/1‑93 (MQ=255)
                                                            aCCGATAGGGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCTGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGcc                             >  2:3540994/1‑93 (MQ=255)
                                                              cGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCaa                           >  1:2928293/1‑93 (MQ=255)
                                                                aTAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCaaaa                         >  1:7094096/1‑93 (MQ=255)
                                                                aTAGTGATTAACGCAACTAACTGGCAAATATGGGGAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCaaaa                         >  2:2646852/1‑93 (MQ=255)
                                                                  aGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGt                       <  1:8492042/93‑1 (MQ=255)
                                                                  aGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGt                       >  1:480770/1‑93 (MQ=255)
                                                                   gTGATTAACGCAACTAACTGGCAGATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGtt                      <  1:2799528/93‑1 (MQ=255)
                                                                   attaTTAACACAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGtt                      <  2:7418335/90‑1 (MQ=255)
                                                                     gATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGc                    <  2:5849928/93‑1 (MQ=255)
                                                                       ttAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTg                  <  2:10025636/93‑1 (MQ=255)
                                                                       ttAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTg                  >  1:3217071/1‑93 (MQ=255)
                                                                         accGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTc                <  2:483810/91‑1 (MQ=255)
                                                                         aaCGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTc                <  1:8280961/93‑1 (MQ=255)
                                                                          aCGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCa               <  1:1362511/93‑1 (MQ=255)
                                                                          aCGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCa               >  1:1680103/1‑93 (MQ=255)
                                                                           cGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCAt              >  1:980036/1‑93 (MQ=255)
                                                                            gCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCAtt             >  1:6148261/1‑93 (MQ=255)
                                                                              aactaactGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTaa           <  2:2301241/93‑1 (MQ=255)
                                                                              aactaactGGCAAATATGGGTAAAAAAATCATCCAGATGCCCGAGTTTTTCTGCGCGAAGATATGCATAATTGCCAAAAGTTGCTGTCATTaa           >  2:5610556/1‑93 (MQ=255)
                                                                               actaactGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAac          >  2:7863453/1‑93 (MQ=255)
                                                                               actaactGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAac          >  1:8372367/1‑93 (MQ=255)
                                                                                ctaactGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAaca         <  1:6294780/93‑1 (MQ=255)
                                                                                ctaactGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAaca         <  1:5716289/93‑1 (MQ=255)
                                                                                ctaactGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAaca         <  2:4694491/93‑1 (MQ=255)
                                                                                ctaactGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAaca         <  2:1539387/93‑1 (MQ=255)
                                                                                 taactGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAacac        <  1:5522770/93‑1 (MQ=255)
                                                                                 taactGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAacac        <  2:3708571/93‑1 (MQ=255)
                                                                                   actGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAACACGc      <  1:2260154/93‑1 (MQ=255)
                                                                                    ctGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAACACGcc     <  2:9860824/93‑1 (MQ=255)
                                                                                    ctGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAACACGcc     <  2:1192028/93‑1 (MQ=255)
                                                                                    ctGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAACACGcc     >  1:9503756/1‑93 (MQ=255)
                                                                                      ggCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAACACGCCaa   <  1:6832391/93‑1 (MQ=255)
                                                                                      ggCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAACACGCCaa   <  2:6562738/93‑1 (MQ=255)
                                                                                       gCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAACACGCCAAg  >  1:4010369/1‑93 (MQ=255)
                                                                                          |                                                                                         
GGGATAGATCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAGATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAACACGCCAAG  >  NC_000913/3725090‑3725269

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: