Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,469,675 A→G S136S (AGT→AGC yfdL ← putative RNA polymerase beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,469,6750AG86.9% 604.2 / 71.9 245S136S (AGT→AGCyfdLputative RNA polymerase beta
Reads supporting (aligned to +/- strand):  ref base A (18/14);  new base G (98/115);  total (116/129)
Fisher's exact test for biased strand distribution p-value = 3.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.87e-01

Only 100 of 247 total aligned reads displayed.

AATGAAGCACCGATACCGGGTACGCCTGCTCTGGAAATAAATTTCACCGGGTCCTGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCACTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAATGCAATCGGAAAATAGGTTACTGTGTCCTGAGACGTCAG  >  NC_000913/2469583‑2469767
                                                                                            |                                                                                            
 aTGAAGCACCGATACCGGGTACGCCTGCTCTGGAAATAAATTTCACCGGGTCCTGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCAc                                                                                             >  2:614923/1‑93 (MQ=255)
  tGAAGCACCGATACCGGGTACGCCTGCTCTGGAAATAAATTTCACCGGGTCCTGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCGCt                                                                                            >  2:5346709/1‑93 (MQ=255)
     aGCACCGATACCGGGTACGCCTGCTCTGGAAATAAATTTCACCGGGTCCTGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCGCTGcc                                                                                         >  1:6264233/1‑93 (MQ=255)
       cACCGATACCGGGTACGCCTGCTCTGGAAATAAATTTCACCGGGTCCTGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCGCTGCCgg                                                                                       >  2:6217339/1‑93 (MQ=255)
        aCCGATACCGGGTACGCCTGCTCTGGAAATAAATTTCACCGGGTCCTGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCGCTGCCggg                                                                                      >  1:4905434/1‑93 (MQ=255)
         ccGATACCGGGTACGCCTGCTCTGGAAATAAATTTCACCGGGTCCTGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCGCTGCCGGGa                                                                                     >  2:3328659/1‑93 (MQ=255)
           gATACCGGGTACGCCTGCTCTGGAAATAAATTTCACCGGGTCCTGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCGCTGCCGGGAgt                                                                                   >  1:5461890/1‑93 (MQ=255)
           gATACCGGGTACGCCTGCTCTGGAAATAAATTTCACCGGGTCCTGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCGCTGCCGGGAgt                                                                                   <  2:5575768/93‑1 (MQ=255)
            aTACCGGGTACGCCTGCTCTGGAAATAAATTTCACCGGGTCCTGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCGCTGCCGGGAgtg                                                                                  >  2:406819/1‑93 (MQ=255)
             tACCGGGTACGCCTGCTCTGGAAATAAATTTCACCGGGTCCTGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCGCTGCCGGGAgtgt                                                                                 >  1:5945901/1‑93 (MQ=255)
              aCCGGGTACGCCTGCTCTGGAAATAAATTTCACCGGGTCCTGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTa                                                                                <  2:2137113/93‑1 (MQ=255)
                 gggTACGCCTGCTCTGGAAATAAATTTCACCGGGTCCTGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATc                                                                             <  1:4741778/93‑1 (MQ=255)
                 gggTACGCCTGCTCTGGAAATAAATTTCACCGGGTCCTGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATc                                                                             <  1:1883708/93‑1 (MQ=255)
                  ggTACGCCTGCTCTGGAAATAAATTTCACCGGGTCCTGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATcc                                                                            <  2:4953195/93‑1 (MQ=255)
                  ggTACGCCTGCTCTGGAAATAAATTTCACCGGGTCCTGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATcc                                                                            <  1:1403294/93‑1 (MQ=255)
                       gCCTGCTCTGGAAATAAATTTCACCGGGTCCTGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGt                                                                       >  1:801284/1‑93 (MQ=255)
                       gCCTGCTCTGGAAATAAATTTCACCGGGTCCTGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCACTGCCGGGAGTGTAATCCTGAGt                                                                       >  2:7312021/1‑93 (MQ=255)
                        ccTGCTCTGGAAATAAATTTCACCGGGTCCTGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTc                                                                      <  1:3970811/93‑1 (MQ=255)
                         cTGCTCTGGAAATAAATTTCACCGGGTCCTGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCg                                                                     >  1:509145/1‑93 (MQ=255)
                         cTGCTCTGGAAATAAATTTCACCGGGTCCTGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCg                                                                     >  2:3814763/1‑93 (MQ=255)
                             tctGGAAATAAATTTCACCGGGTCCTGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCaaa                                                                 <  2:376273/93‑1 (MQ=255)
                              ctGGAAATAAATTTCACCGGGTCCTGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAAt                                                                <  2:660418/93‑1 (MQ=255)
                              ctGGAAATAAATTTCACCGGGTCCTGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAAt                                                                >  2:2969717/1‑93 (MQ=255)
                               tGGAAATAAATTTCACCGGGTCCTGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATa                                                               >  2:4293183/1‑93 (MQ=255)
                                ggAAATAAATTTCACCGGGTCCTGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATAc                                                              >  1:7075854/1‑93 (MQ=255)
                                 gAAATAAATTTCACCGGGTCCTGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATAcc                                                             <  1:6036056/93‑1 (MQ=255)
                                   aataaatTTCACCGGGTCCTGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCg                                                           >  2:4997291/1‑93 (MQ=255)
                                     taaatTTCACCGGGTCCTGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCa                                                         >  1:4212568/1‑93 (MQ=255)
                                     taaatTTCACCGGGTCCTGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCa                                                         <  2:5644274/93‑1 (MQ=255)
                                        atTTCACCGGGTCCTGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGAc                                                      >  2:2783590/1‑93 (MQ=255)
                                        atTTCACCGGGTCCTGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGAc                                                      >  2:1517275/1‑93 (MQ=255)
                                         tttCACCGGGTCCTGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACa                                                     <  2:467923/93‑1 (MQ=255)
                                             aCCGGGTCCTGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTa                                                 <  1:466914/93‑1 (MQ=255)
                                              ccGGGTCCTGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTAt                                                <  1:6444614/93‑1 (MQ=255)
                                               cGGGTCCTGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTAtt                                               <  2:5770499/93‑1 (MQ=255)
                                               cGGGTCCTGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCACTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTAtt                                               <  2:1602333/93‑1 (MQ=255)
                                                gggTCCTGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTc                                              <  2:6980733/93‑1 (MQ=255)
                                                gggTCCTGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTc                                              <  1:349162/93‑1 (MQ=255)
                                                gggTCCTGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCACTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTc                                              <  1:748656/93‑1 (MQ=255)
                                                gggTCCTGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCACTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTc                                              >  2:2361533/1‑93 (MQ=255)
                                                 ggTCCTGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCg                                             <  1:6447052/93‑1 (MQ=255)
                                                   tCCTGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCACTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGa                                           >  1:7219213/1‑93 (MQ=255)
                                                    ccTGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGaa                                          >  1:1781485/1‑93 (MQ=255)
                                                     cTGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGaaa                                         <  2:5276915/93‑1 (MQ=255)
                                                     cTGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGaaa                                         >  2:261444/1‑93 (MQ=255)
                                                     cTGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGaaa                                         <  1:1096291/93‑1 (MQ=255)
                                                      tGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATCCCCGCAGAAATCTATTCGGAAAt                                        <  2:6458529/93‑1 (MQ=255)
                                                      tGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAAt                                        <  2:830432/93‑1 (MQ=255)
                                                      tGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAAt                                        <  2:6864149/93‑1 (MQ=255)
                                                      tGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAAt                                        <  1:4080934/93‑1 (MQ=255)
                                                      tGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCACTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAAt                                        <  1:3701124/93‑1 (MQ=255)
                                                         ttATAACCGGCACAAGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAATGCa                                     >  2:7292083/1‑93 (MQ=255)
                                                         ttATAACCGGCACAAGCTATATAGCCAACATTTGCACTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAATGCa                                     >  2:1808350/1‑93 (MQ=255)
                                                           ataACCGGCACAAGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAATGCAAt                                   >  2:7167151/1‑93 (MQ=255)
                                                            taACCGGCACAAGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAATGCAATc                                  >  2:2085996/1‑93 (MQ=255)
                                                             aaCCGGCACAAGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAATGCAATCg                                 <  1:2128210/93‑1 (MQ=255)
                                                             aaCCGGCACAAGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAATGCAATCg                                 >  1:1609555/1‑93 (MQ=255)
                                                              aCCGGCACAAGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAATGCAATCgg                                >  1:4422061/1‑93 (MQ=255)
                                                               ccGGCACAAGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAATGCAATCGGa                               >  2:5965821/1‑93 (MQ=255)
                                                                 ggCACAAGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAATGCAATCGGaaa                             >  2:2910842/1‑93 (MQ=255)
                                                                  gCACAAGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAATGCAATCGGaaaa                            <  1:865563/93‑1 (MQ=255)
                                                                    acaAGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAATGCAATCGGAAAATa                          <  1:3948387/93‑1 (MQ=255)
                                                                    acaAGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAATGCAATCGGAAAATa                          <  2:5472661/93‑1 (MQ=255)
                                                                    acaAGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAATGCAATCGGAAAATa                          <  1:5735348/93‑1 (MQ=255)
                                                                    acaAGCTATATAGCCAACATTTGCACTGCCGGGAGTGTAATCCTGAGTCGCAGATACCCGCAGACATCTATTCGGAAATGCAATCGGAAAATa                          <  2:3681113/93‑1 (MQ=255)
                                                                     caAGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAATGCAATCGGAAAATAg                         >  2:6380508/1‑93 (MQ=255)
                                                                     caAGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAATGCAATCGGAAAATAg                         >  1:5702516/1‑93 (MQ=255)
                                                                      aaGCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAATGCAATCGGAAAATAgg                        <  1:2563945/93‑1 (MQ=255)
                                                                       aGCTATATAGCCAACATTTGCACTGCCGGGAGGGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAATGCAATCGGAAAATAGGt                       >  1:7400843/1‑93 (MQ=255)
                                                                        gCTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAATGCAATCGGAAAATAGGtt                      <  1:3224175/93‑1 (MQ=255)
                                                                        gCTATATAGCCAACATTTGCACTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAATGCAATCGGAAAATAGGtt                      <  2:4356344/93‑1 (MQ=255)
                                                                         cTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAATGCAATCGGAAAATAGGTTa                     <  2:7085779/93‑1 (MQ=255)
                                                                         cTATATAGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAATGCAATCGGAAAATAGGTTa                     >  1:3397483/1‑93 (MQ=255)
                                                                          tatataGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAATGCAATCGGAAAATAGGTTAc                    <  1:1102095/93‑1 (MQ=255)
                                                                           atataGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAATGCAATCGGAAAATAGGTTACt                   <  1:2713502/93‑1 (MQ=255)
                                                                           atataGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAATGCAATCGGAAAATAGGTTACt                   >  1:1129360/1‑93 (MQ=255)
                                                                           atataGCCAACATTTGCACTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAATGCAATCGGAAAATAGGTTACt                   >  2:3205816/1‑93 (MQ=255)
                                                                            tataGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAATGCAATCGGAAAATAGGTTACtg                  >  1:6693165/1‑93 (MQ=255)
                                                                            tataGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAATGCAATCGGAAAATAGGTTACtg                  >  2:7496972/1‑93 (MQ=255)
                                                                             ataGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAATGCAATCGGAAAATAGGTTACtgt                 <  2:1365418/93‑1 (MQ=255)
                                                                              taGCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAATGCAATCGGAAAATAGGTTACtgtg                <  2:1268074/93‑1 (MQ=255)
                                                                                gCCAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAATGCAATCGGAAAATAGGTTACTGTGTc              <  1:3746199/93‑1 (MQ=255)
                                                                                gCCAACATTTGCACTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAATGCAATCGGAAAATAGGTTACTGTGTc              <  1:5374054/93‑1 (MQ=255)
                                                                                 ccAACATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAATGCAATCGGAAAATAGGTTACTGTGTcc             <  2:2242437/93‑1 (MQ=255)
                                                                                    aCATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAATGCAATCGGAAAATAGGTTACTGTGTCCTga          <  1:4534094/93‑1 (MQ=255)
                                                                                    aCATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAATGCAATCGGAAAATAGGTTACTGTGTCCTga          <  1:5614828/93‑1 (MQ=255)
                                                                                    aCATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAATGCAATCGGAAAATAGGTTACTGTGTCCTga          >  2:5107098/1‑93 (MQ=255)
                                                                                    aCATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAATGCAATCGGAAAATAGGTTACTGTGTCCTga          >  2:3155777/1‑93 (MQ=255)
                                                                                     cATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAATGCAATCGGAAAATAGGTTACTGTGTCCTgag         >  1:3481968/1‑93 (MQ=255)
                                                                                     cATTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAATGCAATCGGAAAATAGGTTACTGTGTCCTgag         <  1:1706559/93‑1 (MQ=255)
                                                                                      aTTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAATGCAATCGGAAAATAGGTTACTGTGTCCTgaga        <  2:1881839/93‑1 (MQ=255)
                                                                                      aTTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAATGCAATCGGAAAATAGGTTACTGTGTCCTgaga        <  2:1411221/93‑1 (MQ=255)
                                                                                      aTTTGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAATGCAATCGGAAAATAGGTTACTGTGTCCTgaga        <  1:6295073/93‑1 (MQ=255)
                                                                                       tttGCGCTGCCGGGAGTGTAATTCTGAGTCGCAAATTCCGGGAGATATCTATTCGGAAATGCAATCGGAAAATAGGTTACTGTGTCCTGAGAc       >  1:2243323/1‑93 (MQ=255)
                                                                                       tttGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAATGCAATCGGAAAATAGGTTACTGTGTCCTGAGAc       <  1:6134395/93‑1 (MQ=255)
                                                                                       tttGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAATGCAATCGGAAAATAGGTTACTGTGTCCTGAGAc       <  1:3591454/93‑1 (MQ=255)
                                                                                        ttGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAATGCAATCGGAAAATAGGTTACTGTGTCCTGAGACg      <  1:5914646/93‑1 (MQ=255)
                                                                                         tGCGCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAATGCAATCGGAAAATAGGTTACTGTGTCCTGAGACGt     >  2:679702/1‑93 (MQ=255)
                                                                                         tGCACTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAATGCAATCGGAAAATAGGTTACTGTGTCCTGAGACGt     <  1:7416950/93‑1 (MQ=255)
                                                                                            aCTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAATGCAATCGGAAAATAGGTTACTGTGTCCTGAGACGTCAg  >  1:2212641/1‑93 (MQ=255)
                                                                                            |                                                                                            
AATGAAGCACCGATACCGGGTACGCCTGCTCTGGAAATAAATTTCACCGGGTCCTGGTTATAACCGGCACAAGCTATATAGCCAACATTTGCACTGCCGGGAGTGTAATCCTGAGTCGCAAATACCCGCAGACATCTATTCGGAAATGCAATCGGAAAATAGGTTACTGTGTCCTGAGACGTCAG  >  NC_000913/2469583‑2469767

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