Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,583,503 G→A Q428* (CAA→TAA)  hsdR ← endonuclease R Type I restriction enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,583,5030GA99.6% 2056.5 / NA 533Q428* (CAA→TAA) hsdRendonuclease R Type I restriction enzyme
Reads supporting (aligned to +/- strand):  ref base G (0/2);  major base A (259/269);  minor base C (0/2);  total (259/274)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 571 total aligned reads displayed.

GATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAG  >  NC_000913/4583405‑4583600
                                                                                                  |                                                                                                 
gATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGGTAAGAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGCCAAACCACTGATTCTGTGTTTACg                                                                                                 >  1:4704346/1‑101 (MQ=255)
 aTTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACgg                                                                                                <  1:13425885/101‑1 (MQ=255)
   tGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGtt                                                                                              >  2:1258926/1‑101 (MQ=255)
     ccTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGtnnn                                                                                            >  2:10031634/1‑98 (MQ=255)
      cTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGc                                                                                           <  1:13194145/101‑1 (MQ=255)
        tttCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTg                                                                                         <  2:265102/101‑1 (MQ=255)
         nnnaaCCGCGCTGAAGGCATATTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCCCTGATTCTGTTTTTACGGTTCGCTGc                                                                                        <  2:13917832/98‑1 (MQ=255)
          tCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGcc                                                                                       <  1:1619402/101‑1 (MQ=255)
          nnnaCCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGcc                                                                                       <  2:7263060/98‑1 (MQ=255)
            aaCCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGa                                                                                     >  1:3374856/1‑101 (MQ=255)
            aaCCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGa                                                                                     <  2:5520169/101‑1 (MQ=255)
                nnncgGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCAt                                                                                 <  2:7732462/98‑1 (MQ=255)
                 cgcgGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACTGTTCGCTGCCGAGCAtt                                                                                >  2:6863604/1‑101 (MQ=255)
                  gcgGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCAttt                                                                               <  1:5134403/101‑1 (MQ=255)
                  gcgGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCAttt                                                                               <  1:1854780/101‑1 (MQ=255)
                   cgGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTc                                                                              >  1:8587190/1‑101 (MQ=255)
                    nnncGGCATCTTCCTGATAATAACGCCGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCAttttc                                                                             <  2:12108030/98‑3 (MQ=255)
                     nnnggCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCca                                                                            <  2:8900655/98‑1 (MQ=255)
                     gACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCca                                                                            <  2:6389994/101‑1 (MQ=255)
                     gACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCca                                                                            <  1:6371075/101‑1 (MQ=255)
                       cGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTccnnn                                                                          >  2:9794797/1‑98 (MQ=255)
                        ggCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCATACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCcagca                                                                         <  1:13274150/101‑1 (MQ=255)
                         gCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCcagnnn                                                                        >  2:12384725/1‑98 (MQ=255)
                         gCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCcagcag                                                                        >  2:1172418/1‑101 (MQ=255)
                          nnncTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCagcagc                                                                       <  2:8510801/98‑1 (MQ=255)
                          cATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCagcagc                                                                       >  1:11174837/1‑101 (MQ=255)
                          cATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCagcagc                                                                       >  1:6526011/1‑101 (MQ=255)
                           nnnttCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGct                                                                      <  2:14708020/98‑1 (MQ=255)
                           aTCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCcagcannn                                                                      >  2:8169993/1‑98 (MQ=255)
                           aTCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGct                                                                      <  1:14567522/101‑1 (MQ=255)
                           aTCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGct                                                                      <  1:12379606/101‑1 (MQ=255)
                              nnncTGATAATCACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCtt                                                                   <  2:8351014/98‑1 (MQ=255)
                                 cTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCg                                                                >  2:2142129/1‑101 (MQ=255)
                                  nnntaataaCGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCgg                                                               <  2:9473284/98‑1 (MQ=255)
                                   gATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCggg                                                              >  1:6170723/1‑101 (MQ=255)
                                    nnnataaCGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGc                                                             <  2:13705677/98‑1 (MQ=255)
                                    ataataACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGc                                                             >  1:3145693/1‑101 (MQ=255)
                                    ataataACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGc                                                             >  1:13900910/1‑101 (MQ=255)
                                    ataataACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGc                                                             <  1:8055423/101‑1 (MQ=255)
                                     taataaCGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCg                                                            >  1:10016613/1‑101 (MQ=255)
                                       nnnaCGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGt                                                          <  2:11502617/98‑1 (MQ=255)
                                         nnngCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGc                                                        <  2:12891355/98‑1 (MQ=255)
                                          aCGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGcc                                                       >  1:14163251/1‑101 (MQ=255)
                                           cGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGGGCCa                                                      >  1:7655055/1‑101 (MQ=255)
                                             cAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGccnnn                                                    >  2:13304063/1‑98 (MQ=255)
                                             cAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCAcc                                                    >  1:2378359/1‑100 (MQ=255)
                                              aGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTc                                                   <  2:6071891/101‑1 (MQ=255)
                                              aGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTc                                                   <  1:1184989/101‑1 (MQ=255)
                                               nnnccATCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCg                                                  <  2:10718631/98‑1 (MQ=255)
                                               ggCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCg                                                  <  1:13771592/101‑1 (MQ=255)
                                                 cccAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCggg                                                <  1:5589990/101‑1 (MQ=255)
                                                  ccAGCTCGCTCATGCCAGGGTTATCGTCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGt                                               <  1:3983147/101‑1 (MQ=255)
                                                  ccAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGt                                               <  2:4324587/101‑1 (MQ=255)
                                                  ccAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGt                                               >  1:1012018/1‑101 (MQ=255)
                                                  ccAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGt                                               <  1:4862669/101‑1 (MQ=255)
                                                   cAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCggnnn                                              >  2:9117635/1‑98 (MQ=255)
                                                    cgctcgctcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTaa                                             <  2:2498949/100‑1 (MQ=255)
                                                      ctcgctcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAgg                                           >  1:10367564/1‑101 (MQ=255)
                                                       tcgctcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTaannn                                          >  2:13488749/1‑98 (MQ=255)
                                                       tcgctcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGc                                          >  1:7357527/1‑101 (MQ=255)
                                                       nnnctcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGc                                          <  2:14034775/98‑1 (MQ=255)
                                                        cgctcATGCCAGGGTTATCGGCAAACCCCTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCt                                         <  1:8696309/101‑1 (MQ=255)
                                                         gctcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAggnnn                                        >  2:7516520/1‑98 (MQ=255)
                                                         gctcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCtt                                        >  1:12872350/1‑101 (MQ=255)
                                                          ctcATTCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCttt                                       <  2:6642463/101‑1 (MQ=255)
                                                          ctcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGCAAGGcnnn                                       >  2:8630828/1‑98 (MQ=255)
                                                           tcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTc                                      >  2:5794446/1‑101 (MQ=255)
                                                           tcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTc                                      <  2:3920735/101‑1 (MQ=255)
                                                              tGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGAc                                   <  1:9401974/101‑1 (MQ=255)
                                                                cccgggTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACat                                 <  2:4072464/98‑1 (MQ=255)
                                                                 cAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACata                                <  2:3568293/101‑1 (MQ=255)
                                                                  aGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGCTGCCACTCGGGTAAGGCTTTCGACannn                               >  2:10394603/1‑98 (MQ=255)
                                                                     gctATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACatatnnn                            >  2:10942061/3‑98 (MQ=255)
                                                                     gTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGc                            >  1:2822418/1‑101 (MQ=255)
                                                                     gTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGc                            <  1:12085326/101‑1 (MQ=255)
                                                                     gTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGc                            <  1:7209366/101‑1 (MQ=255)
                                                                        nnnggCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCggg                         <  2:7355685/98‑1 (MQ=255)
                                                                           ggCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCgggnnn                      >  2:14322906/1‑98 (MQ=255)
                                                                           ggCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGGGCCACTCGGGTAAGGCTTTCGACATATTGCAGGTAt                      >  1:208101/1‑101 (MQ=255)
                                                                             cAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATca                    <  2:3374856/101‑1 (MQ=255)
                                                                               aaCCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGGAAGGCTTTCGACATATTGCGGGTATcacg                  >  1:4512948/1‑101 (MQ=255)
                                                                                aCCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATcannn                 >  2:8749260/1‑98 (MQ=255)
                                                                                    cTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGt             >  1:9758974/1‑101 (MQ=255)
                                                                                     tGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCCCTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTc            <  2:2932569/101‑1 (MQ=255)
                                                                                     tGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTc            <  1:1468069/101‑1 (MQ=255)
                                                                                     tGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGGAAGGCTTTCGACATATTGCGGGTATCACGCACGTc            >  2:4713356/1‑101 (MQ=255)
                                                                                     tGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGGAAGGCTTTCGACATATTGCGGGTATCACGCACGTc            >  1:11356790/1‑101 (MQ=255)
                                                                                      gATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTcg           >  1:8866014/1‑101 (MQ=255)
                                                                                      gATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTcg           <  1:9117635/101‑1 (MQ=255)
                                                                                       aTTCTGTTTTTACGGTTCGCTGCCGAGCATTTCTAGCAGCTCTTCCGGGCGGGGCCGCTCGGGTAAGGCTTTCGGCATATTGCGGACATCACGCAACTcgc          >  2:1561994/1‑101 (MQ=255)
                                                                                        ttCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTcgcg         >  1:5879168/1‑101 (MQ=255)
                                                                                        ttCTGCTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTcgcg         <  2:7014146/101‑1 (MQ=255)
                                                                                         tCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGCGCCACTCGGGTAAGGCTTTCGACATATTTCGGGGATCACGCACGTCGCGa        >  2:5102164/1‑101 (MQ=255)
                                                                                           tGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATa      <  1:2142129/101‑1 (MQ=255)
                                                                                            gTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATAc     <  1:8305027/101‑1 (MQ=255)
                                                                                            gTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACGCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCAcgccgctgttt     <  2:4469368/101‑6 (MQ=255)
                                                                                              ttttACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATannn   >  2:7920805/1‑98 (MQ=255)
                                                                                              ttttACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATannn   >  2:10192687/1‑98 (MQ=255)
                                                                                              ttttACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCa   <  1:3647070/101‑1 (MQ=255)
                                                                                              ttttACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCa   >  1:705069/1‑101 (MQ=255)
                                                                                                  |                                                                                                 
GATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAG  >  NC_000913/4583405‑4583600

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