Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,165,630 G→A 5.1% intergenic (+109/‑38) yegP → / → yegQ UPF0339 family protein/putative peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,165,6300GA 5.1% 258.2 / 3.4 78intergenic (+109/‑38)yegP/yegQUPF0339 family protein/putative peptidase
Reads supporting (aligned to +/- strand):  ref base G (35/39);  new base A (2/2);  total (37/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01

AGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATGCTGTTTATGCGGGCCAGCCGCGTTATTCC  >  NC_000913/2165482‑2165772
                                                                                                                                                    |                                                                                                                                              
aGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAgg                                                                                                                                                >  4:480088/1‑149 (MQ=255)
 gCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGt                                                                                                                                               >  3:782153/1‑149 (MQ=255)
  ccAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGt                                                                                                                                               >  1:396779/1‑148 (MQ=255)
   cAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTCAAGTGGGGGCACTCCCCTAACCGCTTCATTAGGTGa                                                                                                                                             >  3:602034/1‑149 (MQ=255)
   cAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTg                                                                                                                                              >  1:304259/1‑148 (MQ=255)
       ggCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCACCTGATGAAGc                                                                                                                                          <  1:682084/148‑1 (MQ=255)
       ggCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCACCTGATGAAGCg                                                                                                                                         <  4:184298/149‑1 (MQ=255)
              gCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCTTCAGGTGAAGCGGATCTGa                                                                                                                                  <  3:532052/149‑1 (MQ=255)
              gCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGa                                                                                                                                  <  3:304743/149‑1 (MQ=255)
               ccAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGa                                                                                                                                  >  2:243481/1‑148 (MQ=255)
                    cGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGtc                                                                                                                            <  4:459443/149‑1 (MQ=255)
                    cGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGtc                                                                                                                            <  3:730048/149‑1 (MQ=255)
                    cGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGt                                                                                                                             <  1:92523/148‑1 (MQ=255)
                      gTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGtca                                                                                                                           <  1:646378/148‑1 (MQ=255)
                           aGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGaa                                                                                                                     <  4:843953/149‑1 (MQ=255)
                           aGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGaa                                                                                                                     >  3:840968/1‑149 (MQ=255)
                           aGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGa                                                                                                                      <  1:295529/148‑1 (MQ=255)
                           aGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGa                                                                                                                      >  2:200516/1‑148 (MQ=255)
                            gACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGaa                                                                                                                     <  1:877734/148‑1 (MQ=255)
                              cAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACga                                                                                                                  <  3:704374/149‑1 (MQ=255)
                                 tACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACgagag                                                                                                               >  3:522286/1‑149 (MQ=255)
                                     tAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAAt                                                                                                            <  2:304259/148‑1 (MQ=255)
                                     tAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCACCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAAAACGAGAGAAtt                                                                                                           <  3:480088/149‑1 (MQ=255)
                                      aaTTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTa                                                                                                          >  3:789038/1‑149 (MQ=255)
                                         tACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGt                                                                                                       <  3:81280/149‑1 (MQ=255)
                                            gATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTa                                                                                                    >  4:602991/1‑149 (MQ=255)
                                              tGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTaa                                                                                                   <  2:386389/148‑1 (MQ=255)
                                                ccGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAAcc                                                                                                <  3:602991/149‑1 (MQ=255)
                                                           ttgttgGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCttt                                                                                     >  4:715030/1‑149 (MQ=255)
                                                             gttgGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCATGTGAAGCGGAt                                                                                                                                      >  4:655498/1‑98 (MQ=255)
                                                             gttgGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGAt                                                                                                                                      <  3:655498/98‑1 (MQ=255)
                                                               tgGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTccc                                                                                  <  2:826693/148‑1 (MQ=255)
                                                                    tCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCcggcgg                                                                            <  4:522286/149‑1 (MQ=255)
                                                                    tCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCcggcg                                                                             >  1:518664/1‑148 (MQ=255)
                                                                          gcTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGc                                                                      >  4:336318/1‑149 (MQ=255)
                                                                            tttACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTg                                                                    >  4:842744/1‑149 (MQ=255)
                                                                                      ccccccACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGc                                                          <  3:800327/149‑1 (MQ=255)
                                                                                      ccccccACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGc                                                          <  4:782153/149‑1 (MQ=255)
                                                                                         cccACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGt                                                        <  2:285841/148‑1 (MQ=255)
                                                                                           cacaGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTAc                                                     <  3:715030/149‑1 (MQ=255)
                                                                                               gcgcTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCt                                                  >  1:678775/1‑148 (MQ=255)
                                                                                                cgcTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCtt                                                 <  1:155095/148‑1 (MQ=255)
                                                                                                   tACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGc                                             >  3:473455/1‑149 (MQ=255)
                                                                                                    aCAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGc                                             >  1:406036/1‑148 (MQ=255)
                                                                                                       aTGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTa                                          >  2:718044/1‑148 (MQ=255)
                                                                                                        tGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTAt                                         <  2:678775/148‑1 (MQ=255)
                                                                                                        tGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATg                                        <  4:473455/149‑1 (MQ=255)
                                                                                                          cccgcccTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGgc                                      >  3:835404/1‑149 (MQ=255)
                                                                                                           ccgcccTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACgaga                                                                                                                >  3:454787/1‑74 (MQ=255)
                                                                                                           ccgcccTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACgaga                                                                                                                <  4:454787/74‑1 (MQ=255)
                                                                                                           ccgcccTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGtt                                                                                                      <  3:647221/84‑1 (MQ=255)
                                                                                                           ccgcccTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGtt                                                                                                      >  4:647221/1‑84 (MQ=255)
                                                                                                           ccgcccTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGgcg                                     >  4:215643/1‑149 (MQ=255)
                                                                                                               ccTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAAc                                                                                                 <  3:54901/85‑1 (MQ=255)
                                                                                                               ccTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAAc                                                                                                 >  4:54901/1‑85 (MQ=255)
                                                                                                                cTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGAt                                <  4:602034/149‑1 (MQ=255)
                                                                                                                 ttAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATg                               >  3:650084/1‑149 (MQ=255)
                                                                                                                 ttAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATg                               >  3:298875/1‑149 (MQ=255)
                                                                                                                  tAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATg                               >  1:829175/1‑148 (MQ=255)
                                                                                                                   aaaGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATGCt                             <  4:227469/149‑1 (MQ=255)
                                                                                                                     aGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATGCTg                            <  1:184268/148‑1 (MQ=255)
                                                                                                                                tCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATGCTGTTTATGCGGGcc                >  3:107018/1‑149 (MQ=255)
                                                                                                                                tCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCATTCGCTTATGGCGCAGATGCTGTTTATGCGGGc                 >  1:705708/1‑148 (MQ=255)
                                                                                                                                  cccTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTc                                               >  4:435865/1‑116 (MQ=255)
                                                                                                                                  cccTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTc                                               <  3:435865/116‑1 (MQ=255)
                                                                                                                                   ccTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATGCTGTTTATGCGGGCCAg              >  1:791700/1‑148 (MQ=255)
                                                                                                                                   ccTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATGCTGTTTATGCGGGCCAg              >  2:913648/1‑148 (MQ=255)
                                                                                                                                   ccTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATGCTGTTTATGCGGGCCAg              >  2:453575/1‑148 (MQ=255)
                                                                                                                                       tccGCTTCACCTGATGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATGCTGTTTATGCGGGCCAGCcgcg         >  4:29391/2‑149 (MQ=255)
                                                                                                                                       tccGCTTCACCTGATGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATGCTGTTTATGCGGGCCAGCcgcg         >  3:94969/2‑149 (MQ=255)
                                                                                                                                         cGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATGCTGTTTATGCGGGCCAGCCGCGtt       >  3:250270/1‑149 (MQ=255)
                                                                                                                                         cGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATGCTGTTTATGCGGGCCAGCCGCGtt       <  4:835404/149‑1 (MQ=255)
                                                                                                                                          gCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATGCTGTTTATGCGGGCCAGCCGCGTTa      <  4:298875/149‑1 (MQ=255)
                                                                                                                                          gCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATGCTGTTTATGCGGGCCAGCCGCGTTa      <  4:650084/149‑1 (MQ=255)
                                                                                                                                              catcaGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATGCTGTTTATGCGGGCCAGCCGCGTTATTcc  <  3:215643/149‑1 (MQ=255)
                                                                                                                                               atcaGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCttt                                                                                     >  4:586686/1‑65 (MQ=255)
                                                                                                                                               atcaGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCttt                                                                                     <  3:586686/65‑1 (MQ=255)
                                                                                                                                               atcaGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATGCTGTTTATGCGGGCCAGCCGCGTTATTcc  <  2:829175/148‑1 (MQ=255)
                                                                                                                                               atcaGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATGCTGTTTATGCGGGCCAGCCGCGTTATTcc  <  1:200516/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
AGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATGCTGTTTATGCGGGCCAGCCGCGTTATTCC  >  NC_000913/2165482‑2165772

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