Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,090,046 G→A *17* (TAG→TAA hisL → his operon leader peptide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,090,0460GA99.5% 623.3 / NA 190*17* (TAG→TAAhisLhis operon leader peptide
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (93/95);  minor base T (0/1);  total (94/96)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 190 total aligned reads displayed.

TATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAG  >  NC_000913/2089915‑2090179
                                                                                                                                   |                                                                                                                                     
tataAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCa                                                                                                                                >  2:876433/1‑139 (MQ=255)
tataAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCa                                                                                                                                >  2:134266/1‑139 (MQ=255)
 ataAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAg                                                                                                                               <  1:461860/139‑1 (MQ=255)
     aaaaGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGa                                                                                                                           >  2:48002/1‑139 (MQ=255)
      aaaGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAt                                                                                                                          <  2:897785/139‑1 (MQ=255)
      aaaGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAt                                                                                                                          <  2:413241/139‑1 (MQ=255)
       aaGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAtg                                                                                                                         <  1:966516/139‑1 (MQ=255)
        aGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAt                                                                                                                          >  2:188150/1‑137 (MQ=255)
         gCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAtgtg                                                                                                                       <  2:246492/139‑1 (MQ=255)
           ccTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAtgtgtg                                                                                                                     >  2:55871/1‑139 (MQ=255)
             ttGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCt                                                                                                                   >  2:719736/1‑139 (MQ=255)
               gCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTgg                                                                                                                 <  1:732232/139‑1 (MQ=255)
               gCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTgg                                                                                                                 >  2:671564/1‑139 (MQ=255)
                 tttCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGaa                                                                                                               <  2:36641/139‑1 (MQ=255)
                   tCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGa                                                                                                             >  2:334715/1‑139 (MQ=255)
                        cGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTc                                                                                                        <  1:945615/139‑1 (MQ=255)
                            aaaGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGaa                                                                                                               <  2:986720/128‑1 (MQ=255)
                             aaGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATc                                                                                                   <  1:140772/139‑1 (MQ=255)
                                   tttAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAg                                                                                             <  1:318376/139‑1 (MQ=255)
                                   tttAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAg                                                                                             <  1:98391/139‑1 (MQ=255)
                                   tttAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAg                                                                                             <  2:468873/139‑1 (MQ=255)
                                     tAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAtg                                                                                                                         >  2:131388/1‑109 (MQ=255)
                                          tAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAt                                                                                                                          <  2:299244/103‑1 (MQ=255)
                                          tAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAt                                                                                                                          >  1:299244/1‑103 (MQ=255)
                                           aaaaGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCAt                                                                                     <  1:671564/139‑1 (MQ=255)
                                            aaaGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATGTTTCAGGCGAt                                                                                                                          <  2:567975/101‑1 (MQ=255)
                                            aaaGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAt                                                                                                                          >  1:567975/1‑101 (MQ=255)
                                            aaaGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGaa                                                                                                               <  2:213662/112‑1 (MQ=255)
                                            aaaGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAg                                                                                             >  2:228149/1‑130 (MQ=255)
                                            aaaGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAg                                                                                             <  1:228149/130‑1 (MQ=255)
                                            aaaGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATg                                                                                    >  2:557171/1‑139 (MQ=255)
                                             aaGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGa                                                                                                     >  2:699579/1‑121 (MQ=255)
                                             aaGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGa                                                                                                     <  1:699579/121‑1 (MQ=255)
                                             aaGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGa                                                                                   <  1:557171/139‑1 (MQ=255)
                                             aaGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGa                                                                                   <  1:742429/139‑1 (MQ=255)
                                              aGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAg                                                                                             <  2:325628/128‑1 (MQ=255)
                                                  aTCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCg                                                                                                                            >  2:650350/1‑93 (MQ=255)
                                                  aTCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAt                                                                                                                          <  2:771994/95‑1 (MQ=255)
                                                  aTCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAt                                                                                                                          >  1:771994/1‑95 (MQ=255)
                                                  aTCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCa                                                                              >  1:42175/1‑139 (MQ=255)
                                                  aTCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAAAGCa                                                                              >  1:203618/1‑139 (MQ=255)
                                                  aTCAGTTGAATAAACATCCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCg                                                                                                                            <  1:650350/93‑1 (MQ=255)
                                                       ttGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCt                                                                                                                   <  1:270035/97‑1 (MQ=255)
                                                       ttGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATgaga                                                                         <  2:434567/139‑1 (MQ=255)
                                                       ttGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATgaga                                                                         <  2:42175/139‑1 (MQ=255)
                                                          aaTAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAg                                                                                             >  1:471335/1‑116 (MQ=255)
                                                          aaTAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAg                                                                                             >  1:624309/1‑116 (MQ=255)
                                                          aaTAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAg                                                                      <  1:54320/139‑1 (MQ=255)
                                                              aaCATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGcccc                                                                  <  2:476299/139‑1 (MQ=255)
                                                               aCATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGccccc                                                                 <  2:263119/139‑1 (MQ=255)
                                                               aCATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGAGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGccccc                                                                 >  1:830964/1‑139 (MQ=255)
                                                                 aTTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCg                                                                                                                            >  2:175184/1‑78 (MQ=255)
                                                                 aTTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGaa                                                                                                               <  1:260542/91‑1 (MQ=255)
                                                                 aTTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCgg                                                               >  2:453633/1‑139 (MQ=255)
                                                                 aTTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCgg                                                               >  2:742411/1‑139 (MQ=255)
                                                                   tCACAGAGACTTTTATAACCCGCGTACAATTTAAACACCACCATCATCACCATCATCCTGACTACTCTTTCAGGCGAACTGTGCTACACGACATTCAGCTCTTCCAGTGGTTCATGAACGCATGAGAAAACCACCGTaa                                                             >  1:573925/1‑139 (MQ=255)
                                                                    cacaGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGTGATGTGTGCTGGAAg                                                                                                              >  1:332055/1‑89 (MQ=255)
                                                                    cacaGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAg                                                                                             <  2:525739/106‑1 (MQ=255)
                                                                    cacaGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAg                                                                                             >  1:525739/1‑106 (MQ=255)
                                                                         agaCTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAt                                                                                                                          <  1:386094/72‑1 (MQ=255)
                                                                         agaCTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCa                                                                                                       >  1:1004368/1‑91 (MQ=255)
                                                                          gaCTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCg                                                                                                                            >  2:130080/1‑69 (MQ=255)
                                                                          gaCTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCg                                                                                                                            <  1:130080/69‑1 (MQ=255)
                                                                           aCTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGaa                                                                                  <  2:950911/110‑1 (MQ=255)
                                                                               ttATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCttt                                                                                                                                  >  1:912728/1‑58 (MQ=255)
                                                                                   gACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGa                                                              >  2:817613/1‑122 (MQ=255)
                                                                                   gACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGa                                                              <  1:817613/122‑1 (MQ=255)
                                                                                       cgcgTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGa                                                                                                                <  1:186625/68‑1 (MQ=255)
                                                                                       cgcgTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAg                                                                                             >  2:934062/1‑87 (MQ=255)
                                                                                       cgcgTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGccccc                                                                 >  1:401516/1‑115 (MQ=255)
                                                                                       cgcgTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGcccc                                                                  <  2:166018/114‑1 (MQ=255)
                                                                                       cgcgTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGcccc                                                                  >  1:166018/1‑114 (MQ=255)
                                                                                       cgcgTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATATGTGCTGGAAGACATTCAGATCTTCCAg                                                                                             <  2:625627/87‑1 (MQ=255)
                                                                                         cgTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAg                                                                      >  1:626068/1‑108 (MQ=255)
                                                                                          gTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAt                                                                                                                          >  1:988363/1‑55 (MQ=255)
                                                                                                  tAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTc                                                                                                        >  1:132841/1‑65 (MQ=255)
                                                                                                        ccaccaTCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATgaga                                                                         <  1:477226/90‑1 (MQ=255)
                                                                                                        ccaccaTCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGAt                        <  1:498731/139‑1 (MQ=255)
                                                                                                        ccaccaTCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGAt                        >  2:445836/1‑139 (MQ=255)
                                                                                                         caccatcatcaccatcatCCTGACTAATCTTTCAg                                                                                                                               >  1:506256/1‑35 (MQ=255)
                                                                                                         caccatcatcaccatcatCCTGACTAATCTTTCAg                                                                                                                               <  2:506256/35‑1 (MQ=255)
                                                                                                          accatcatcaccatcatcCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATaa                      >  2:971117/1‑139 (MQ=255)
                                                                                                              tcatcaccatcatcCTGACTAATCTTTCAGGCGATGTGTGCTGGaa                                                                                                               >  2:82877/1‑46 (MQ=255)
                                                                                                                   accatcatcCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGACCGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAgc             >  1:220463/1‑138 (MQ=255)
                                                                                                                   accatcatcCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGa             >  2:310467/1‑139 (MQ=255)
                                                                                                                   accatcatcCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGa             <  1:340856/139‑1 (MQ=255)
                                                                                                                   accatcatcCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGa             <  2:754164/139‑1 (MQ=255)
                                                                                                                   accatcatcCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGa             <  1:445836/139‑1 (MQ=255)
                                                                                                                    ccatcatcCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGAc            <  1:310467/139‑1 (MQ=255)
                                                                                                                    ccatcatcCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGATGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGAc            <  1:237017/139‑1 (MQ=255)
                                                                                                                     catcatcCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCAt                                                                                     <  2:589989/65‑1 (MQ=255)
                                                                                                                     catcatcCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCAt                                                                                     >  1:589989/1‑65 (MQ=255)
                                                                                                                     catcatcCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACa           >  2:542927/1‑139 (MQ=255)
                                                                                                                      atcatcCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAg          >  2:393157/1‑139 (MQ=255)
                                                                                                                      atcatcCTGACTAATCTTTCAAGCGGTGTGGGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAg          <  2:565784/139‑1 (MQ=255)
                                                                                                                        catcCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGt        >  1:445701/1‑139 (MQ=255)
                                                                                                                          tcCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGttt      >  1:858313/1‑139 (MQ=255)
                                                                                                                            cTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTaa    >  2:603910/1‑139 (MQ=255)
                                                                                                                              gACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAag  <  1:603910/139‑1 (MQ=255)
                                                                                                                               aCTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTcc                                                                                               >  2:297492/1‑45 (MQ=255)
                                                                                                                                   |                                                                                                                                     
TATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAG  >  NC_000913/2089915‑2090179

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