Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,458,811 C→T A51V (GCG→GTG)  hofG → putative general protein secretion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,458,8110CT82.3% 677.7 / 103.0 273A51V (GCG→GTG) hofGputative general protein secretion protein
Reads supporting (aligned to +/- strand):  ref base C (33/15);  new base T (114/110);  total (148/125)
Fisher's exact test for biased strand distribution p-value = 2.61e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.74e-01

Only 100 of 279 total aligned reads displayed.

TTATCGGCGTACTTGCCAGCCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGCGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAGTAGAAGCGCCGACAC  >  NC_000913/3458721‑3458903
                                                                                          |                                                                                            
ttATCGGCGTACTTGCCAGCCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGgtgc                                                                                            >  1:5002562/1‑93 (MQ=255)
ttATCGGCGTACTTGCCAGCCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGgtgc                                                                                            >  1:7060420/1‑93 (MQ=255)
 tATCGGCGTACTTGCCAGCCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGCGCt                                                                                           >  1:7296427/1‑93 (MQ=255)
  aTCGGCGTACTTGCCAGCCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTg                                                                                          >  2:5856936/1‑93 (MQ=255)
      gCGTACTTGCCAGCCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGaaa                                                                                      >  2:751789/1‑93 (MQ=255)
      gCGTACTTGCCAGCCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGCGCTGGaaa                                                                                      <  2:4560777/93‑1 (MQ=255)
       cGTACTTGCCAGCCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGaaaa                                                                                     >  2:7337679/1‑93 (MQ=255)
        gTACTTGCCAGCCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAAc                                                                                    <  2:5944072/93‑1 (MQ=255)
         tACTTGCCAGCCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACg                                                                                   >  2:2429844/1‑93 (MQ=255)
          aCTTGCCAGCCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGc                                                                                  >  1:6178942/1‑93 (MQ=255)
           cTTGCCAGCCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGcc                                                                                 <  1:1041982/93‑1 (MQ=255)
           cTTGCCAGCCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGcc                                                                                 >  1:1810528/1‑93 (MQ=255)
            ttGCCAGCCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGccc                                                                                >  2:2204005/1‑93 (MQ=255)
            ttGCCAGCCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGccc                                                                                >  1:4640621/1‑93 (MQ=255)
              gCCAGCCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCtt                                                                              <  2:937766/93‑1 (MQ=255)
                cAGCGTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGa                                                                            <  2:3632639/93‑1 (MQ=255)
                cAGCCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGa                                                                            <  1:73847/93‑1 (MQ=255)
                 aGCCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGat                                                                           <  1:168692/93‑1 (MQ=255)
                 aGCCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGat                                                                           >  1:1465703/1‑93 (MQ=255)
                  gCCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGata                                                                          <  2:2837564/93‑1 (MQ=255)
                  gCCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGata                                                                          >  2:325677/1‑93 (MQ=255)
                    cTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATg                                                                        >  2:4847368/1‑93 (MQ=255)
                     tggtggTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGt                                                                       >  1:1109331/1‑93 (MQ=255)
                     tggtggTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGt                                                                       <  1:6545417/93‑1 (MQ=255)
                           ttCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAAc                                                                 >  2:5960395/1‑93 (MQ=255)
                               taacctaaTGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGa                                                             >  2:3812609/1‑93 (MQ=255)
                                  cctaaTGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACaa                                                          <  2:4807972/93‑1 (MQ=255)
                                  cctaaTGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACaa                                                          >  2:2054109/1‑93 (MQ=255)
                                  cctaaTGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACaa                                                          <  1:1633414/93‑1 (MQ=255)
                                   ctaaTGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAat                                                         <  2:4115011/93‑1 (MQ=255)
                                   ctaaTGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAat                                                         >  1:2052285/1‑93 (MQ=255)
                                   ctaaTGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGCGCTGGAAAACGCCCTTGATATGTACAAACTCGACAat                                                         >  1:6905621/1‑93 (MQ=255)
                                    taaTGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAGCTCGACAatc                                                        <  2:29244/93‑1 (MQ=255)
                                    taaTGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGCGCTGGAAAACGCCCTTGATATGTACAAACTCGACAatc                                                        <  1:2540691/93‑1 (MQ=255)
                                      aTGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGNNNNNGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAatcat                                                      <  2:6063739/93‑1 (MQ=255)
                                        gggCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAAtcatca                                                    <  1:6148519/93‑1 (MQ=255)
                                           cAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGCGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTa                                                 <  2:6645602/93‑1 (MQ=255)
                                            aaTAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTAt                                                >  1:2784244/1‑93 (MQ=255)
                                             aTAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATc                                               <  2:3101441/93‑1 (MQ=255)
                                             aTAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATc                                               >  1:804902/1‑93 (MQ=255)
                                              tAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGCGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATcc                                              >  2:5380337/1‑93 (MQ=255)
                                              tAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGCGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATcc                                              >  2:4345253/1‑93 (MQ=255)
                                              tAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGCGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATcc                                              >  1:4080234/1‑93 (MQ=255)
                                                aaGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCaa                                            >  1:2287162/1‑93 (MQ=255)
                                                   aaaaGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACca                                         >  2:6874498/1‑93 (MQ=255)
                                                    aaaGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACcac                                        >  1:5514717/1‑93 (MQ=255)
                                                    aaaGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACcac                                        >  2:3723188/1‑93 (MQ=255)
                                                     aaGGCGGATAAGCAAAAAGCCGTCAACGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACcaca                                       <  2:737891/93‑1 (MQ=255)
                                                       ggCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACaaa                                     >  1:6466453/1‑93 (MQ=255)
                                                       ggCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACaaa                                     <  2:6577218/93‑1 (MQ=255)
                                                       ggCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACaaa                                     >  2:1079938/1‑93 (MQ=255)
                                                        gCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAAt                                    <  1:380495/93‑1 (MQ=255)
                                                         cGGATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATc                                   <  1:3107874/93‑1 (MQ=255)
                                                         cGGATAAGCAAAAAGCCGTCAGCGATATTGTGGCGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATc                                   <  1:1519497/93‑1 (MQ=255)
                                                          ggATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCa                                  >  2:1151688/1‑93 (MQ=255)
                                                          ggATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCa                                  <  2:1484011/93‑1 (MQ=255)
                                                           gATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAg                                 >  2:6299371/1‑93 (MQ=255)
                                                           gATAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAg                                 >  1:625837/1‑93 (MQ=255)
                                                            aTAAGCAAAAAGCCGTCAGCGATATTGTGGCGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAgg                                >  2:2615474/1‑93 (MQ=255)
                                                             tAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAggg                               <  1:3155922/93‑1 (MQ=255)
                                                             tAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAggg                               >  1:6754772/1‑93 (MQ=255)
                                                             tAAGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAggg                               <  1:1259184/93‑1 (MQ=255)
                                                              aaGCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAgggg                              <  1:4193571/93‑1 (MQ=255)
                                                                gCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCt                            <  1:5509216/93‑1 (MQ=255)
                                                                gCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCt                            >  2:3513426/1‑93 (MQ=255)
                                                                gCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCt                            >  2:514433/1‑93 (MQ=255)
                                                                gCAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAAACAaggggt                            >  2:5537436/1‑91 (MQ=255)
                                                                 cAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCtt                           <  2:1860313/93‑1 (MQ=255)
                                                                 cAAAAAGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCtt                           <  1:2137064/93‑1 (MQ=255)
                                                                  aaaaaGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTg                          >  2:7116171/1‑93 (MQ=255)
                                                                  aaaaaGCCGTCAGCGATATTGTGGCGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTg                          >  2:494510/1‑93 (MQ=255)
                                                                   aaaaGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGa                         >  1:2835803/1‑93 (MQ=255)
                                                                   aaaaGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGa                         <  2:972535/93‑1 (MQ=255)
                                                                   aaaaGCCGTCAGCGATATTGTGGCGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGa                         >  1:3773228/1‑93 (MQ=255)
                                                                     aaGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAAt                       <  1:1345802/93‑1 (MQ=255)
                                                                      aGCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATc                      >  1:3425573/1‑93 (MQ=255)
                                                                       gCCGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATcc                     >  1:3881130/1‑93 (MQ=255)
                                                                         cGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCtt                   <  2:163548/93‑1 (MQ=255)
                                                                         cGTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCtt                   >  2:5052602/1‑93 (MQ=255)
                                                                         cGTCAGCGATATTGTGGCGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCtt                   <  2:5679430/93‑1 (MQ=255)
                                                                          gTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTta                  <  1:5782680/93‑1 (MQ=255)
                                                                          gTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTta                  <  1:4422937/93‑1 (MQ=255)
                                                                          gTCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTta                  <  2:6472929/93‑1 (MQ=255)
                                                                           tCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTtag                 >  1:7007005/1‑93 (MQ=255)
                                                                           tCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTtag                 <  1:5263052/93‑1 (MQ=255)
                                                                           tCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTtag                 <  1:5217258/93‑1 (MQ=255)
                                                                           tCAGCGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTtag                 >  1:2717498/1‑93 (MQ=255)
                                                                           tCAGCGATATTGTGGCGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTtag                 >  2:6434041/1‑93 (MQ=255)
                                                                               cGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAGTAGa             <  1:93701/93‑1 (MQ=255)
                                                                               cGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAGTAGa             <  1:5954968/93‑1 (MQ=255)
                                                                               cGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAGTAGa             >  1:482645/1‑93 (MQ=255)
                                                                               cGATATTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAGTAGa             <  1:4686117/93‑1 (MQ=255)
                                                                                 atatTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAGTAGAAg           >  1:4897102/1‑93 (MQ=255)
                                                                                  tatTGTGGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAGTAGAAgc          <  2:2330568/93‑1 (MQ=255)
                                                                                  tatTGTGGCGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAGTAGAAgc          >  1:6063676/1‑93 (MQ=255)
                                                                                  tatTGTGGCGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAGTAGAAgc          >  2:3373340/1‑93 (MQ=255)
                                                                                   atTGTGGCGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAGTAGAAgcg         >  1:7478034/1‑93 (MQ=255)
                                                                                    ttGTGGCGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAGTAGAAgcgc        >  2:205696/1‑93 (MQ=255)
                                                                                     tgtgGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAGTAGAAGCGcc       <  2:2520750/93‑1 (MQ=255)
                                                                                       tgGTGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAGTAGAAGCGCCGa     >  1:5350924/1‑93 (MQ=255)
                                                                                          |                                                                                            
TTATCGGCGTACTTGCCAGCCTGGTGGTTCCTAACCTAATGGGCAATAAAGAAAAGGCGGATAAGCAAAAAGCCGTCAGCGATATTGTGGCGCTGGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAGTAGAAGCGCCGACAC  >  NC_000913/3458721‑3458903

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