Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,486,771 T→C 5.8% intergenic (+19/‑20) 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,7710TC 5.8% 333.0 / 8.0 104intergenic (+19/‑20)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 (47/51);  new base C (4/2);  total (51/53)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.50e-01

Only 100 of 106 total aligned reads displayed.

CCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGTTAATGCTTGGCGAATTACGATTCGGTCTGCT  >  NC_000913/3486626‑3486918
                                                                                                                                                 |                                                                                                                                                   
ccGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTa                                                                                                                                                   <  4:630129/148‑1 (MQ=255)
 cGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTAc                                                                                                                                                  <  3:269487/148‑1 (MQ=255)
 cGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTAc                                                                                                                                                  <  4:55413/148‑1 (MQ=255)
    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:133014/149‑1 (MQ=255)
            tGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGCTACCAGGTAacgcgc                                                                                                                                      <  1:712979/149‑1 (MQ=255)
             cgTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCCGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAgcgcgc                                                                                                                                      <  3:94346/147‑1 (MQ=255)
                cAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCAt                                                                                                                                   >  3:483797/1‑148 (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)
                              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)
                                                                               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                                                                                                                     >  4:206301/1‑89 (MQ=255)
                                                                                         cacaCGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGa                                                                                                                     <  3:206301/89‑1 (MQ=255)
                                                                                            acGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTAt                                                      >  1:544974/1‑149 (MQ=255)
                                                                                            acGGTAAAACCATCGTCGTTTACGGCACTCGCTAATCCCGTCGGAGTGGCGCGTTAc                                                                                                                                                  <  4:281015/57‑1 (MQ=255)
                                                                                            acGGTAAAACCATCGTCGTTTACGGCACTCGCTAATCCCGTCGGAGTGGCGCGTTAc                                                                                                                                                  >  3:281015/1‑57 (MQ=255)
                                                                                                tAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGttt                                                  >  1:848806/1‑149 (MQ=255)
                                                                                                 aaaaCCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGttt                                                  >  3:1041743/1‑148 (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                                                        <  2:302452/137‑1 (MQ=255)
                                                                                                      cATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCt                                                        >  1:302452/1‑137 (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                                                             <  4:370802/118‑1 (MQ=255)
                                                                                                                    gCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTg                                                             >  3:370802/1‑118 (MQ=255)
                                                                                                                     cACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGcc                                               >  2:283096/1‑131 (MQ=255)
                                                                                                                     cACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGcc                                               <  1:283096/131‑1 (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                                                                               <  1:760833/93‑1 (MQ=255)
                                                                                                                           ttAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTAt                                                                               >  2:760833/1‑93 (MQ=255)
                                                                                                                           ttAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGTTAATGCTTg                        <  3:867925/148‑1 (MQ=255)
                                                                                                                             aaTCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCTGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTg                                           >  1:290989/1‑127 (MQ=255)
                                                                                                                             aaTCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTg                                           <  2:290989/127‑1 (MQ=255)
                                                                                                                                cccGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGTTAATGCTTGGCGaa                   >  4:252197/1‑148 (MQ=255)
                                                                                                                                  cGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGTTAATGCTTGGCGAAtt                 <  4:854415/148‑1 (MQ=255)
                                                                                                                                   gTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGTTAATGCTTGGCGAATTa                >  3:163618/1‑148 (MQ=255)
                                                                                                                                   gTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGTTAATGCTTGGCGAATTa                <  4:329906/148‑1 (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)
                                                                                                                                                 |                                                                                                                                                   
CCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGCGCAGACTGATTTATCACCCCGATATCAACTATGCACTTCGACAAACGCTGGTGCTATGTTTGCCCGTGGCCGTTGGGTTAATGCTTGGCGAATTACGATTCGGTCTGCT  >  NC_000913/3486626‑3486918

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: