Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,165,626 T→C 5.3% intergenic (+105/‑42) 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,6260TC 5.3% 251.8 / 3.2 75intergenic (+105/‑42)yegP/yegQUPF0339 family protein/putative peptidase
Reads supporting (aligned to +/- strand):  ref base T (35/36);  new base C (2/2);  total (37/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.34e-01

TAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATGCTGTTTATGCGGGCCAGCCGCGTTATTCC  >  NC_000913/2165479‑2165772
                                                                                                                                                   |                                                                                                                                                  
tAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTtcat                                                                                                                                                    <  1:765559/148‑1 (MQ=255)
 aaaGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTtcatca                                                                                                                                                  <  3:642578/149‑1 (MQ=255)
   aGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAgg                                                                                                                                                >  4:480088/1‑149 (MQ=255)
   aGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAg                                                                                                                                                 <  1:512996/148‑1 (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)
                                                                                                                                                   |                                                                                                                                                  
TAAAGCCAATGGCACAAGCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATGCTGTTTATGCGGGCCAGCCGCGTTATTCC  >  NC_000913/2165479‑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: