Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,814,809 A→G K169K (AAA→AAG yicC → putative alpha helix protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,814,8090AG98.9% 560.8 / ‑4.3 174K169K (AAA→AAGyicCputative alpha helix protein
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (92/80);  total (92/82)
Fisher's exact test for biased strand distribution p-value = 2.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 182 total aligned reads displayed.

GACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGAAGATGCTCAGGTGCAACTGGAAAACAA  >  NC_000913/3814720‑3814898
                                                                                         |                                                                                         
gacTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTc                                                                                        <  1:3665868/93‑1 (MQ=255)
gacTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTc                                                                                        >  2:573583/1‑93 (MQ=255)
   tttATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGc                                                                                     <  2:3281053/93‑1 (MQ=255)
    ttATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCAGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCt                                                                                    >  2:4628792/1‑93 (MQ=255)
      aTTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTcc                                                                                  <  2:3205809/93‑1 (MQ=255)
      aTTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTcc                                                                                  <  2:4276850/93‑1 (MQ=255)
      aTTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTcc                                                                                  <  2:277531/93‑1 (MQ=255)
      aTTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTcc                                                                                  <  1:1113274/93‑1 (MQ=255)
          tCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCata                                                                              >  2:2556988/1‑93 (MQ=255)
          tCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCata                                                                              <  1:4482016/93‑1 (MQ=255)
           cgcgcgcgAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCatat                                                                             <  1:1247060/93‑1 (MQ=255)
              gcgcgAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCACTATGcc                                                                          >  2:2438835/1‑93 (MQ=255)
               cgcgAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCg                                                                         >  2:3368441/1‑93 (MQ=255)
                 cgAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGTCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGa                                                                       >  1:4167118/1‑93 (MQ=255)
                  gAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGaa                                                                      <  1:3675686/93‑1 (MQ=255)
                  gAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGaa                                                                      <  1:3238952/93‑1 (MQ=255)
                   aaaCCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCTCCTCCATTTGCCGGaaa                                                                     >  2:2392510/1‑93 (MQ=255)
                   aaaCCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGaaa                                                                     >  1:2540763/1‑93 (MQ=255)
                   aaaCCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGaaa                                                                     >  1:2313553/1‑93 (MQ=255)
                    aaCCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGaaac                                                                    >  2:3631640/1‑92 (MQ=255)
                    aaCCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAAt                                                                    <  2:3044532/93‑1 (MQ=255)
                    aaCCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAAt                                                                    <  1:2430849/93‑1 (MQ=255)
                    aaCCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAAt                                                                    >  1:4240851/1‑93 (MQ=255)
                      ccGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATcc                                                                  <  1:2580022/93‑1 (MQ=255)
                       cGAAGGTCAGGCGCTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCt                                                                 <  1:3585503/93‑1 (MQ=255)
                       cGAAGGTCAGGCACTGAAAGCATTGATCGATCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCt                                                                 <  1:4308696/93‑1 (MQ=255)
                       cGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCt                                                                 >  1:2746252/1‑93 (MQ=255)
                         aaGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGc                                                               <  2:4075976/93‑1 (MQ=255)
                         aaGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGc                                                               >  1:3130808/1‑93 (MQ=255)
                         aaGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGc                                                               >  1:1424045/1‑93 (MQ=255)
                           ggTCAGGCACTGAAAGCATTAATCGAGCAGCGTCTGTGAAGCGTCACCGCCGAAGTGGTCAAGGTCGGGTCCCATATGCCGGAAATCCTGCaa                                                             >  2:2946972/1‑93 (MQ=255)
                            gTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAAt                                                            <  1:4513820/93‑1 (MQ=255)
                            gTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAAt                                                            <  2:1356622/93‑1 (MQ=255)
                            gTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAAt                                                            <  2:2696118/93‑1 (MQ=255)
                            gTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAAt                                                            <  2:2413245/93‑1 (MQ=255)
                               aGGCACTGAAAGCATTGATCGATCAGCGTCTGGAATGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGc                                                         >  1:851136/1‑93 (MQ=255)
                               aGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGc                                                         >  1:1027468/1‑93 (MQ=255)
                               aGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGc                                                         <  2:539590/93‑1 (MQ=255)
                               aGGCACTGAAAGCATTGATCGAGCAGCGTCTCGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAACACTGCAATGGc                                                         >  2:2780742/1‑93 (MQ=255)
                                ggCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCa                                                        <  1:592700/93‑1 (MQ=255)
                                ggCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCa                                                        <  1:1983592/93‑1 (MQ=255)
                                ggCACGGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCa                                                        <  2:186254/93‑1 (MQ=255)
                                 gCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAg                                                       <  2:2152187/93‑1 (MQ=255)
                                   aCTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCg                                                     <  1:3458006/93‑1 (MQ=255)
                                   aCTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCg                                                     >  1:2839140/1‑93 (MQ=255)
                                   aCTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCg                                                     >  2:920623/1‑93 (MQ=255)
                                    cTGAAAGCATTGATCGAGCAGCGTCTGTAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGt                                                    <  2:67384/93‑1 (MQ=255)
                                    cTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCgg                                                    >  1:2631654/1‑92 (MQ=255)
                                    cTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGt                                                    >  2:2053912/1‑93 (MQ=255)
                                     tGAAAGCATTGATCGAGCATCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTg                                                   >  2:2176198/1‑93 (MQ=255)
                                     tGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTg                                                   >  2:3829914/1‑93 (MQ=255)
                                     tGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTg                                                   >  1:4373340/1‑93 (MQ=255)
                                     tGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTg                                                   <  1:4562501/93‑1 (MQ=255)
                                      gAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGa                                                  <  2:655439/93‑1 (MQ=255)
                                      gAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGa                                                  >  2:1937384/1‑93 (MQ=255)
                                       aaaGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAg                                                 <  1:1080862/93‑1 (MQ=255)
                                       aaaGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAg                                                 >  2:832827/1‑93 (MQ=255)
                                            aTTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCt                                            >  2:2377201/1‑93 (MQ=255)
                                            aTTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCt                                            >  2:3065986/1‑93 (MQ=255)
                                             ttGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCATCGTGAGCGTCTg                                           <  2:2406931/93‑1 (MQ=255)
                                             ttGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTg                                           <  1:3019894/93‑1 (MQ=255)
                                             ttGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTg                                           <  2:3130808/93‑1 (MQ=255)
                                             ttGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGGCTg                                           >  2:3463532/1‑93 (MQ=255)
                                              tGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTgg                                          <  1:1007562/93‑1 (MQ=255)
                                                aTCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTc                                        >  1:4400543/1‑93 (MQ=255)
                                                   gAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTcgcg                                     >  1:2910821/1‑93 (MQ=255)
                                                     gcagcGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGaa                                   <  2:2604032/93‑1 (MQ=255)
                                                       agcGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGc                                 <  1:1912657/93‑1 (MQ=255)
                                                       agcGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGc                                 >  2:3957853/1‑93 (MQ=255)
                                                        gcGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCt                                <  1:710467/93‑1 (MQ=255)
                                                        gcGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCt                                <  2:3338974/93‑1 (MQ=255)
                                                          gTCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTgg                              >  1:2198732/1‑93 (MQ=255)
                                                           tCTGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGa                             <  2:4376564/93‑1 (MQ=255)
                                                             tGGAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGAAg                           >  2:797684/1‑93 (MQ=255)
                                                              ggAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGAAGa                          <  1:530339/93‑1 (MQ=255)
                                                               gAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGAAGAt                         >  2:949290/1‑93 (MQ=255)
                                                               gAAGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGAAGAt                         <  2:1820312/93‑1 (MQ=255)
                                                                aaGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGAAGATg                        <  2:4047473/93‑1 (MQ=255)
                                                                 aGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGAAGATGc                       >  2:1422932/1‑93 (MQ=255)
                                                                 aGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGAAGATGc                       >  2:3789913/1‑93 (MQ=255)
                                                                 aGGCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGAAGATGc                       >  2:1711874/1‑93 (MQ=255)
                                                                  ggCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGAAGATGCt                      >  2:4652338/1‑93 (MQ=255)
                                                                  ggCGTCACCGCCGAAGTGGTCAAGGTCAGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGAAGATGCt                      >  1:1616139/1‑93 (MQ=255)
                                                                   gCGTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGAAGATGCTc                     >  2:1949001/1‑93 (MQ=255)
                                                                     gTCACCGCCGAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGAAGATGCTCAg                   <  2:3486436/93‑1 (MQ=255)
                                                                         ccgccgAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGAAGATGCTCAGGTGc               >  1:821912/1‑93 (MQ=255)
                                                                         ccgccgAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGAAGATGCTCAGGTGc               >  2:1360481/1‑93 (MQ=255)
                                                                          cgccgAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGAAGATGCTCAGGTGCa              <  2:492378/93‑1 (MQ=255)
                                                                             cgAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGAAGATGCTCAGGTGCAAcg           >  1:1342192/1‑92 (MQ=255)
                                                                             cgAAGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGAAGATGCTCAGGTGCAACt           <  2:1059721/93‑1 (MQ=255)
                                                                                aGTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGAAGATGCTCAGGTGCAACTGGa        >  2:4504582/1‑93 (MQ=255)
                                                                                 gTGGTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGAAGATGCTCAGGTGCAACTGGaa       >  1:225628/1‑93 (MQ=255)
                                                                                   ggTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGAAGATGCTCAGGTGCAACTGGaaaa     >  2:699313/1‑93 (MQ=255)
                                                                                   ggTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGAAGATGCTCAGGTGCAACTGGaaaa     <  1:2888380/93‑1 (MQ=255)
                                                                                   ggTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGAAGATGCTCAGGTGCAACTGGaaaa     >  1:3081563/1‑93 (MQ=255)
                                                                                    gTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGAAGATGCTCAGGTGCAACTGGAAaac    >  2:2343359/1‑93 (MQ=255)
                                                                                    gTCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGAAGATGCTCAGGTGCAACTGGAAaac    >  2:3896264/1‑93 (MQ=255)
                                                                                     tCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGAAGATGCTCAGGTGCAACTGGAAaaca   >  1:403919/1‑93 (MQ=255)
                                                                                     tCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGAAGATGCTCAGGTGCAACTGGAAaaca   >  2:3743148/1‑93 (MQ=255)
                                                                                     tCAAGGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGAAGATGCTCAGGTGCAACTGGAAaaca   >  2:3743147/1‑93 (MQ=255)
                                                                                         |                                                                                         
GACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGAAGATGCTCAGGTGCAACTGGAAAACAA  >  NC_000913/3814720‑3814898

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: