Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,067,320 A→T 7.7% intergenic (‑147/+20) ygfI ← / ← yggE putative LysR‑type DNA‑binding transcriptional regulator YgfI/uncharacterized protein YggE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,067,3200AT 7.7% 179.2 / 6.6 66intergenic (‑147/+20)ygfI/yggEputative LysR‑type DNA‑binding transcriptional regulator YgfI/uncharacterized protein YggE
Reads supporting (aligned to +/- strand):  ref base A (26/34);  new base T (3/2);  total (29/37)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.16e-01

CATCCTTGTCCGCTTAACCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGG  >  NC_000913/3067186‑3067459
                                                                                                                                       |                                                                                                                                           
cATCCTTGTCCGCTTAACCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACATCGTAATTGTTCAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCg                                                                                                                                     <  1:989804/143‑1 (MQ=255)
  tCCTTGTCCGCTTAACCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAAGATCGGGTAAAGGCATTAGCCCGACGTATTCGTCggg                                                                                                                                   <  1:868251/143‑1 (MQ=255)
   ccTTGTCCGCTTAACCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCAGGTAAAGCCATTAGCCCGACGAATACGTCAGGc                                                                                                                                  >  2:74931/1‑143 (MQ=255)
         ccGCATAACCGTTTCGTGTCGTCTATCTTTGGTCAAACCTCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCGGGCTaatgc                                                                                                                            <  2:477254/143‑5 (MQ=255)
          cGCTTAACCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTGTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCAGGCTACACat                                                                                                                           >  1:913396/1‑143 (MQ=255)
                gccGTTTCAGTCGGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATattatt                                                                                                                     <  2:1095151/142‑1 (MQ=255)
                aCCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGACGCAATATCGGGTAAAGGCATTAGCCCGAAGAATACGCCGGGCTCCAAAaagtatc                                                                                                                     >  2:174877/1‑142 (MQ=255)
                aCCGTGTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATattatt                                                                                                                     <  1:301732/143‑1 (MQ=255)
                  cGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAGAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATtgt                                                                                                                   >  2:264846/1‑143 (MQ=255)
                  cGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGCCTGTCTACAAACATTATtgt                                                                                                                   >  2:822955/1‑143 (MQ=255)
                   gTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATtgtg                                                                                                                  <  1:234120/143‑1 (MQ=255)
                    tttCAGTCAGCCTATCCTTGCTGAAACCGCGAGTAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGtgc                                                                                                                 <  1:855901/143‑1 (MQ=255)
                    tttCAGTAAACCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGtgc                                                                                                                 >  2:191763/1‑143 (MQ=255)
                     ttCAGTCAGCCTAGCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCCTGGTTTTACAGCAAATTACAGGAAGAGGAGGCAATAGCGGGTTAAGTCATTAGCCCGACGAATACGTCGGGCTACAAGTGATAGTGtgct                                                                                                                >  2:560134/1‑143 (MQ=255)
                           tcagCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGACAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAgg                                                                                                           >  2:355531/1‑143 (MQ=255)
                                   tctttGATGAAAGCGAGAGGAGAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAg                                                                                                   <  1:1011649/140‑1 (MQ=255)
                                    ccTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGc                                                                                                  <  1:328895/143‑1 (MQ=255)
                                     cTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCg                                                                                                 >  1:568501/1‑143 (MQ=255)
                                     cTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGGCTTTTATCg                                                                                                 >  1:544810/1‑143 (MQ=255)
                                      ttGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCgg                                                                                                <  1:565855/143‑1 (MQ=255)
                                          tgaAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGgttg                                                                                            <  2:827283/143‑1 (MQ=255)
                                           gaAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAGGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGgttgt                                                                                           <  1:327527/143‑1 (MQ=255)
                                           gaAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGgttgt                                                                                           >  1:1040233/1‑143 (MQ=255)
                                           gaAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGGAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGgttgt                                                                                           <  1:1080554/143‑1 (MQ=255)
                                            aaaCCGCGAGCAAGGAGAGGTGATTACGTCATGGTGTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGgttgtt                                                                                          <  2:287078/143‑1 (MQ=255)
                                              aCCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGGGGGTTGTTGa                                                                                        >  1:103450/1‑143 (MQ=255)
                                              aCCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAAGATCGGGTAAAGGCATAAGCAAGACGAATACGTCGGGCTACAAATAGTATTGTGGTGCAGGTGTATTAGCGGGTTGTTGa                                                                                        >  2:957445/1‑143 (MQ=255)
                                                cgcgAGCAAAGATAGGGGATTACGTCGTGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATCAGCCTGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATc                                                                                      <  1:399671/143‑1 (MQ=255)
                                                 gcgAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACATATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATcc                                                                                     >  1:399893/1‑143 (MQ=255)
                                                   caGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCcac                                                                                   <  1:30332/142‑1 (MQ=255)
                                                     gccaataTAGGCGCTTGCGTCATGGTTTTTCAGAGAATTACCTAAAAAGGGTGCTATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAg                                                                                 <  2:828630/137‑1 (MQ=255)
                                                      cAAAGATAGGTGATTACGTCATGGTATTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAgg                                                                                <  1:138913/143‑1 (MQ=255)
                                                             aGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAAc                                                                         <  1:319615/143‑1 (MQ=255)
                                                              ggttgttACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACt                                                                        <  2:514698/138‑1 (MQ=255)
                                                               gTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTg                                                                       >  1:600213/1‑143 (MQ=255)
                                                               gTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATAATTGTGCTGCAGGGGTTTTAGCGGGTTGTTGATCCACAGGGTCTGAcgg                                                                       >  2:753945/1‑141 (MQ=255)
                                                                   ttACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAg                                                                   >  1:301583/1‑143 (MQ=255)
                                                                    tACTTCATTGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGa                                                                  <  1:587029/143‑1 (MQ=255)
                                                                       gTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACca                                                               >  1:850212/1‑143 (MQ=255)
                                                                       gTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTCTGGGTGTTGTTGATCCAGAGGTTCTAACTGGAAggcga                                                               >  1:1041881/1‑139 (MQ=255)
                                                                        tCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTCTTGCGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACcac                                                              <  1:636930/143‑1 (MQ=255)
                                                                        tCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACcac                                                              >  1:413618/1‑143 (MQ=255)
                                                                         cATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACcaca                                                             >  1:595653/1‑143 (MQ=255)
                                                                         cATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACcaca                                                             >  2:63261/1‑143 (MQ=255)
                                                                          aTGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTAGTTTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCGCAGGGTCTCACTGGACGACCACAt                                                            >  2:664635/1‑143 (MQ=255)
                                                                            ggTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCg                                                          <  1:934721/143‑1 (MQ=255)
                                                                              ttttGCAGAAAATTACAGAAAAAGGAGGCAATACAGGGTCATCTCATGAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGAc                                                        <  2:1072250/143‑1 (MQ=255)
                                                                              tttcACAGAAAATTGCAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGAc                                                        <  1:938482/143‑1 (MQ=255)
                                                                                  acagaaaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGa                                                    <  1:472846/143‑1 (MQ=255)
                                                                                  acagaaaattacagaaaaAGGAGGCAATATCGGGACAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCCGCAAGTGTGTTAGCGGGTTGTTGATCCACAGATTCTAACTTGAAGACAACATCGACCGGa                                                    >  2:602060/1‑143 (MQ=255)
                                                                                   cagaaaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGat                                                   <  1:854754/143‑1 (MQ=255)
                                                                                    agaaagttagagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAAGATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGatc                                                  <  2:1061662/143‑1 (MQ=255)
                                                                                          ttacagcaaaAGGAGGCCACCTCGTGTAAAGGCATTCGCCCGACGCATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCaaa                                            <  1:152804/143‑1 (MQ=255)
                                                                                          ttacagaaaaAGGAGGCAATATCTGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCaaa                                            <  2:825074/143‑1 (MQ=255)
                                                                                              agaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAAAACGTCGGGCTACAAATATTATTGTACTGCAGGTGTTTTAGCGGGGTGTTGATCCACAGGTTCTAACTGGAAGACCACATCCACCTGATCATCAAACTGa                                        >  2:282028/1‑143 (MQ=255)
                                                                                                   aaGGAGGCAATAGCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGc                                   <  1:103532/143‑1 (MQ=255)
                                                                                                      taggCAAGACCGCGTAAAGGCAGTAGCCCGACGAATACGTCGGGCTACAACTCTTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGc                                <  1:449193/142‑1 (MQ=255)
                                                                                                       aggCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCCGGTGTTTCAGCGGGTTGTTTATCGACAGGTTCTAACTGGAAGACCACATCGACGTGATCATCAAACTGAATGACGGcg                               >  2:844777/1‑142 (MQ=255)
                                                                                                            atatCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTc                          <  2:301583/143‑1 (MQ=255)
                                                                                                                cGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTaa                      <  2:568501/143‑1 (MQ=255)
                                                                                                                         gCATTAGCCCGCCTAATGCGTCGGGCTACAAAGATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTgg             <  2:656784/143‑1 (MQ=255)
                                                                                                                          cATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTggg            <  1:63261/143‑1 (MQ=255)
                                                                                                                             tAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGTTCATCAAACTGTATAGCGGCCTGCTCGTAAGTTTCCTGGGCgg         <  2:608100/143‑1 (MQ=255)
                                                                                                                               gCCCGACGTATTCGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGGTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGac       >  1:1071762/1‑143 (MQ=255)
                                                                                                                                cccGACGTATTCGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGaca      >  2:801014/1‑143 (MQ=255)
                                                                                                                                cccGACGTATTCGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGCAAGACCCCATCGACCTGATCATCCCACTGAATAGCGGACTCACCGTACGCTTCCTCGGCGGaca      >  2:891817/1‑143 (MQ=255)
                                                                                                                                    acgaATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCgg  >  1:791283/1‑143 (MQ=255)
                                                                                                                                       |                                                                                                                                           
CATCCTTGTCCGCTTAACCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGG  >  NC_000913/3067186‑3067459

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