Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,486,776 T→A 7.1% intergenic (+24/‑15) crp → / → yhfK cAMP‑activated global transcription factor, mediator of catabolite repression/putative transporter, FUSC superfamily inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,486,7760TA 7.1% 340.0 / 13.6 99intergenic (+24/‑15)crp/yhfKcAMP‑activated global transcription factor, mediator of catabolite repression/putative transporter, FUSC superfamily inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (46/46);  new base A (4/3);  total (50/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.22e-01

Only 100 of 103 total aligned reads displayed.

CAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGTTAATGCTTGGCGAATTACGATTCGGTCTGCTCTTC  >  NC_000913/3486630‑3486922
                                                                                                                                                  |                                                                                                                                                  
cAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGCTACCAgg                                                                                                                                                  >  2:133014/1‑149 (MQ=255)
   gAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAg                                                                                                                                               >  2:239483/1‑149 (MQ=255)
      aTTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAgcgc                                                                                                                                            <  2:638545/149‑1 (MQ=255)
      aTTGGTCAGATTGTCGGCTGTTATCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGAGCTCCGCACACGGCAAAACCATCGTCGTTTACGGCACCCGTTAATCCCGTTGGAGTGGCCCGTTATCTGGTAgca                                                                                                                                             <  4:507083/148‑2 (MQ=255)
        tGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGCTACCAGGTAacgcgc                                                                                                                                          <  1:712979/149‑1 (MQ=255)
        tGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGCTACCAGGTAacgcgc                                                                                                                                          <  1:133014/149‑1 (MQ=255)
         cgTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCCGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAgcgcgc                                                                                                                                          <  3:94346/147‑1 (MQ=255)
             aGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCAtt                                                                                                                                      >  3:681653/1‑148 (MQ=255)
              gATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCAtttt                                                                                                                                    <  1:656420/149‑1 (MQ=255)
                ttGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTg                                                                                                                                   <  3:416402/148‑1 (MQ=255)
                 tGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGt                                                                                                                                  <  4:1004044/148‑1 (MQ=255)
                    cGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTc                                                                                                                               >  4:209501/1‑148 (MQ=255)
                     ggCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTcc                                                                                                                              <  4:37890/148‑1 (MQ=255)
                      gCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTccc                                                                                                                             >  3:854415/1‑148 (MQ=255)
                      gCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTccc                                                                                                                             >  3:1021949/1‑148 (MQ=255)
                          ttCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGAt                                                                                                                        <  1:320170/149‑1 (MQ=255)
                          ttCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCATTACCTGGTAGCGCGCCATTTTGTTTCCCCCGAt                                                                                                                        <  2:758569/149‑1 (MQ=255)
                                tGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGgc                                                                                                                   >  4:114970/1‑148 (MQ=255)
                                 gAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCt                                                                                                                                                    >  3:793729/1‑114 (MQ=255)
                                 gAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCt                                                                                                                                                    <  4:793729/114‑1 (MQ=255)
                                     ccGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGAc                                                                                                             >  1:739278/1‑149 (MQ=255)
                                       gTGGGAAGCAGTCTGAAGATGCTGGAAGATCAGAACCGGATCTCCGCACACGGTAAAACCATCGTCGTGTACTGCACTCGTTACTCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTg                                                                                                           <  2:865375/149‑1 (MQ=255)
                                        tGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTCCCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGAAGACTg                                                                                                           <  3:553585/148‑1 (MQ=255)
                                         gggACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGAt                                                                                                         <  2:1005321/149‑1 (MQ=255)
                                          ggACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGAtt                                                                                                        <  1:197555/149‑1 (MQ=255)
                                              gCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTAt                                                                                                     <  4:865178/148‑1 (MQ=255)
                                               cATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATc                                                                                                    <  3:762817/148‑1 (MQ=255)
                                                 ttCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGa                                                                                                          >  1:69228/1‑140 (MQ=255)
                                                 ttCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGa                                                                                                          <  2:69228/140‑1 (MQ=255)
                                                 ttCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCAc                                                                                                  <  3:114970/148‑1 (MQ=255)
                                                             tGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGa                                                                                                          <  1:80317/128‑1 (MQ=255)
                                                             tGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGa                                                                                                          >  2:80317/1‑128 (MQ=255)
                                                               gAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCt                                                                                                                                                    >  1:851226/1‑84 (MQ=255)
                                                               gAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCt                                                                                                                                                    <  2:851226/84‑1 (MQ=255)
                                                                   aTCAGAACCTGATCGCCGCGCACGGTAAAACCATCGTCGTTTACGGCACTCGTTACTCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCAc                                                                               <  2:707375/149‑1 (MQ=255)
                                                                    tCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACt                                                                              >  2:483542/1‑149 (MQ=255)
                                                                      aGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACtt                                                                             >  3:356265/1‑148 (MQ=255)
                                                                           cTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACaa                                                                       >  1:378676/1‑149 (MQ=255)
                                                                           cTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACaa                                                                       >  2:361325/1‑149 (MQ=255)
                                                                           cTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACa                                                                        >  4:881388/1‑148 (MQ=255)
                                                                             gATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAAc                                                                     <  1:483542/149‑1 (MQ=255)
                                                                                 tcCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTg                                                                 <  2:40680/149‑1 (MQ=255)
                                                                                  ccGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTgg                                                                >  1:195389/1‑149 (MQ=255)
                                                                                  ccGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTg                                                                 >  3:756559/1‑148 (MQ=255)
                                                                                     cacaCGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGa                                                                                                                         <  3:206301/89‑1 (MQ=255)
                                                                                     cacaCGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGa                                                                                                                         >  4:206301/1‑89 (MQ=255)
                                                                                        acGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTAt                                                          >  1:544974/1‑149 (MQ=255)
                                                                                            tAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGttt                                                      >  1:848806/1‑149 (MQ=255)
                                                                                             aaaaCCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGttt                                                      >  3:1041743/1‑148 (MQ=255)
                                                                                              aaaCCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTg                                                     <  3:555634/148‑1 (MQ=255)
                                                                                               aaCCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGc                                                    >  3:562097/1‑148 (MQ=255)
                                                                                                aCCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGccc                                                  <  1:408878/149‑1 (MQ=255)
                                                                                                  cATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGAtt                                                                                                        >  2:773694/1‑93 (MQ=255)
                                                                                                  cATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGAtt                                                                                                        <  1:773694/93‑1 (MQ=255)
                                                                                                  cATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCt                                                            >  1:302452/1‑137 (MQ=255)
                                                                                                  cATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCt                                                            <  2:302452/137‑1 (MQ=255)
                                                                                                   aTCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGt                                                <  3:585413/148‑1 (MQ=255)
                                                                                                   aTCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTg                                               <  1:812934/149‑1 (MQ=255)
                                                                                                    tcgtcgTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGCAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTg                                               <  3:129422/148‑1 (MQ=255)
                                                                                                      gtcgtTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGcc                                            >  2:762026/1‑149 (MQ=255)
                                                                                                         gtTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGtt                                         <  2:623871/149‑1 (MQ=255)
                                                                                                         gtTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGt                                          <  4:756559/148‑1 (MQ=255)
                                                                                                          tttACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTTGCCGtt                                         <  4:527092/148‑1 (MQ=255)
                                                                                                          tttACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGtt                                         <  3:171699/148‑1 (MQ=255)
                                                                                                           ttACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTg                                        >  3:913262/1‑148 (MQ=255)
                                                                                                             aCGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGt                                     <  2:901356/149‑1 (MQ=255)
                                                                                                                gCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTg                                                                 >  3:370802/1‑118 (MQ=255)
                                                                                                                gCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTg                                                                 <  4:370802/118‑1 (MQ=255)
                                                                                                                gCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGTTaa                                  >  1:678361/1‑149 (MQ=255)
                                                                                                                 cACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGcc                                                   <  1:283096/131‑1 (MQ=255)
                                                                                                                 cACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGcc                                                   >  2:283096/1‑131 (MQ=255)
                                                                                                                   cTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGTTAATg                                >  3:329906/1‑148 (MQ=255)
                                                                                                                    tCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGTTAATGCt                              >  1:819459/1‑149 (MQ=255)
                                                                                                                     cGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGTTAATGCt                              >  4:550308/1‑148 (MQ=255)
                                                                                                                      gTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTg                                        <  2:766905/137‑1 (MQ=255)
                                                                                                                      gTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTg                                        >  1:766905/1‑137 (MQ=255)
                                                                                                                      gTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGTTAATGCTTg                            >  1:455505/1‑149 (MQ=255)
                                                                                                                       ttAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTAt                                                                                   >  2:760833/1‑93 (MQ=255)
                                                                                                                       ttAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTAt                                                                                   <  1:760833/93‑1 (MQ=255)
                                                                                                                       ttAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGTTAATGCTTg                            <  3:867925/148‑1 (MQ=255)
                                                                                                                         aaTCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCTGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTg                                               >  1:290989/1‑127 (MQ=255)
                                                                                                                            cccGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGTTAATGCTTGGCGaa                       >  4:252197/1‑148 (MQ=255)
                                                                                                                              cGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGTTAATGCTTGGCGAAtt                     <  4:854415/148‑1 (MQ=255)
                                                                                                                               gTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGTTAATGCTTGGCGAATTa                    <  4:329906/148‑1 (MQ=255)
                                                                                                                               gTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGTTAATGCTTGGCGAATTa                    >  3:163618/1‑148 (MQ=255)
                                                                                                                                   gAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGTTAATGCTTGGCGAATTACGAt                <  3:209501/148‑1 (MQ=255)
                                                                                                                                   aaatGGCGCGCTACCAGGTAACGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGTTAATGCTTGGCGAATTACGAt                >  4:881311/4‑148 (MQ=255)
                                                                                                                                    aGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGTTAATGCTTGGCGAATTACGAtt               >  3:576861/1‑148 (MQ=255)
                                                                                                                                     atGGCGCGCTACCAGGTAACGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGTTAATGCTTGGCGAATTACGATTc              >  3:303908/2‑148 (MQ=255)
                                                                                                                                        gcgcgcTACCAGGTAACGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGTTAATGCTTGGCGAATTACGATTCGGt           <  4:303908/148‑1 (MQ=255)
                                                                                                                                        gcgcGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGTTAATGCTTGGCGAATTACGATTCGGTc          >  1:246176/1‑149 (MQ=255)
                                                                                                                                        gcgcGTTACCTGGGAGCGCGCCATTTTGTCTCCCCCGATGTGGCGCAGACTGATTTATCACCCCAATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGTTAATGCTTGGCGAATTACGATTCGGTc          <  1:642219/149‑1 (MQ=255)
                                                                                                                                         cgcgCTACCAGGTAACGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGTTAATGCTTGGCGAATTACGATTCGGTCt         >  1:418124/1‑149 (MQ=255)
                                                                                                                                          gcgTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCGCCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGTTAATGCTTGGCGAATTACGATTCGGTCTg        >  2:878885/1‑149 (MQ=255)
                                                                                                                                          gcgTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGTTAATGCTTGGCGAATTACGATTCGGTCTg        >  1:336268/1‑149 (MQ=255)
                                                                                                                                           cgTTAACTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGCTGCTATGTTTGCCCGTGGCCGTTGGGTTAATGCTTGGCGAATTACGATTCGGTCTg        >  3:931937/1‑148 (MQ=255)
                                                                                                                                             ttACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGTTAATGCTTGGCGAATTACGATTCGGTCTGct      <  3:43180/148‑1 (MQ=255)
                                                                                                                                                 cTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGTTAATGCTTGGc                          <  1:38676/124‑1 (MQ=255)
                                                                                                                                                 cTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGTTAATGCTTGGc                          >  2:38676/1‑124 (MQ=255)
                                                                                                                                                 cTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGTTAATGCTTGGCGAATTACGATTCGGTCTGCtcttc  >  3:698917/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGTTAATGCTTGGCGAATTACGATTCGGTCTGCTCTTC  >  NC_000913/3486630‑3486922

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