Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,993,111 (A)5→4 100% coding (1018/1377 nt) ygeH → putative transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,993,1070A.75.5% 428.7 / 219.2 200coding (1014/1377 nt)ygeHputative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (19/30);  new base . (70/81);  total (89/111)
Fisher's exact test for biased strand distribution p-value = 4.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 211 total aligned reads displayed.

TATGTCTTTGGCACTTCATGGAAAAAGTGAACTTGAACTTGCTGCTCAAAAAGCATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGGAAAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAGTATCTCATATCTTATTTGAACAGGCTAAGATACACTCAACTGATATA  >  NC_000913/2993008‑2993206
                                                                                                   |                                                                                                   
  tGTCTTTGGCACTTCATGGAAAAAGTGAACTTGAACTTGCTGCTCAAAAAGCATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGGaaa                                                                                                   >  2:471213/1‑100 (MQ=255)
  tGTCTTTGGCACTTCATGGAAAAAGTGAACTTGAACTTGCTGCTCAAAAAGCATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGGaaa                                                                                                   <  1:4057666/100‑1 (MQ=255)
    tCTTTGGCACTTCATGGAAAAAGTGAACTTGAACTTGCTGCTCAAAAAGCATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGGaaaat                                                                                                 >  1:1970285/1‑99 (MQ=255)
      cgatctaCTTCATGGAAAAAGTGAACTTGAACTTGCTGCTCAAAAAGCATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGGAAAAttt                                                                                               <  1:1069431/94‑1 (MQ=255)
        tGGCACTTCATGGAAAAAGTGAACTTGAACTTGCTGCTCAAAAAGCATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGGAAAAtttta                                                                                             <  1:2214640/100‑2 (MQ=255)
        tGGCACTTCATGGAAAAAGTGAACTTGAACTTGCTGCTCAAAAAGCATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGGAAAAtttta                                                                                             <  2:1211584/100‑2 (MQ=255)
        tGGCACTTCATGGAAAAAGTGAACTTGAACTTGCTGCTCAAAAAGCATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGGAAAAtttta                                                                                             <  1:994897/100‑2 (MQ=255)
         ggCACTTCATGGAAAAAGTGACCTTGCACTTGCTGCTCCAAACGCATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGG‑AAAATTTTAg                                                                                           <  1:1230604/100‑1 (MQ=255)
           cACTTCATGGAAAAAGTGAACTTGAACTTGCTGCTCAAAAAGCATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGG‑AAAATTTTAg                                                                                           <  2:2702560/98‑1 (MQ=255)
            aCTTCATGGAAAAAGTGAACTTGAACTTGCTGCTCAAAAAGCATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGGaaaatttagatcg                                                                                         >  2:1069427/1‑91 (MQ=255)
             cTTCATGGAAAAAGTGAACTTGAACTTGCTGCTCAAAAAGCATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGG‑AAAATTTTAGCtat                                                                                       >  1:2137848/1‑100 (MQ=255)
                cATGGAAAAAGTGAACTTGAACTTGCTGCTCAAAAAGCATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGG‑AAAATTTTAGCtattat                                                                                    >  1:3558735/1‑100 (MQ=255)
                 aTGGAAAAAGTGAACTTGAACTTGCTGCTCAAAAAGCATTAGAGCTTTTAGATTATGTATCAGCCATAACCACTGTCGATGG‑AAAATTTTAGCTATTATg                                                                                   >  1:668656/1‑100 (MQ=255)
                 aTGGAAAAAGTGAACTTGAACTTGCTGCTCAAAAAGCATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATg                                                                                   <  1:1792761/100‑1 (MQ=255)
                   ggAAAAAGTGAACTTGAACTTGCTGCTCAAAAAGCATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATggg                                                                                 <  1:156799/100‑1 (MQ=255)
                    gAAAAAGTGAACTTGAACTTGCTGCTCAAAAAGCATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGa                                                                                >  2:280667/1‑100 (MQ=255)
                     aaaaaGTGAACTTGAAATAGCAGCTCAAAAAGCATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGAc                                                                               <  2:3567163/100‑1 (MQ=255)
                      aaaaGTGAACTTGAACTTGCTGCTCAAAAAGCATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACt                                                                              <  1:3559692/100‑1 (MQ=255)
                      aaaaGTGAACTTGAACTTGCTGCTCAAAAAGCATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACt                                                                              <  1:2408084/100‑1 (MQ=255)
                       aaaGTGAACTTGAACTTGCTGCTCAAAAAGCATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTg                                                                             <  1:4014340/100‑1 (MQ=255)
                        aaGTGAACTTGAACTTGCTGCTCAAAAAGCATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGa                                                                            <  2:3883071/100‑1 (MQ=255)
                         aGTGAACTTGAACTTGCTGCTCAAAAAGCATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGAt                                                                           <  2:1323250/100‑1 (MQ=255)
                          gTGAACTTGAACTTGCTGCTCAAAAAGCATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGG‑CAAATTTTAGCTATTATGGGACTGATa                                                                          <  2:3768789/100‑1 (MQ=255)
                          gTGAACTTGAACTGGCTGCTCAAAAAGCATTAGAGATTTTAGATTATGTATCAGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATa                                                                          <  2:253250/100‑1 (MQ=255)
                            gaacttgaacttGCTGCTCAAAAAGCATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAAc                                                                        <  2:458277/100‑1 (MQ=255)
                            gaacttgaacttGCTGCTCAAAAAGCATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAAc                                                                        <  2:3848510/100‑1 (MQ=255)
                               cttgaacttgCTGCTCAAAAAGGATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTgg                                                                     >  2:2422136/1‑100 (MQ=255)
                               cttgaacttgCTGCTCAAAAAGCATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTgg                                                                     <  1:226159/100‑1 (MQ=255)
                               cttgaacttgCTGCTCAAAAAGCATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTgg                                                                     <  2:2212746/100‑1 (MQ=255)
                               cttgaacttgCTGCTCAAAAAGCATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTgg                                                                     <  2:2798328/100‑1 (MQ=255)
                                 tgaacttgCTGCTCAAAAAGCATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGGAAAAATTTTAGCTATTATGGGACTGATAACTGgt                                                                    >  2:1853530/1‑100 (MQ=255)
                                 tgaacttgCTGCTCAAAAAGCATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGgtc                                                                   <  1:3010520/100‑1 (MQ=255)
                                    acttgCTGCTCAAAAAGCATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGgtctgt                                                                <  2:3160667/100‑1 (MQ=255)
                                      ttgCTGCTCAAAAAGCATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGgtctgtct                                                              <  2:1592883/100‑1 (MQ=255)
                                      ttgCTGCTCAAAAAGCATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGgtctgtct                                                              <  1:1123670/100‑1 (MQ=255)
                                       tgctgcTCAAAAAGCATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGGAAAAATTTTAGCTATTATGGGACTGATAACTGgtctgtct                                                              >  2:3868878/1‑100 (MQ=255)
                                       tgctgcTCAAAAAGCATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGGtctgtctg                                                             >  2:3384235/1‑100 (MQ=255)
                                        gctgctCAAAAAGCATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGG‑TAAATTTTAGCTATTATGGGACTGATAACTGGtctgtctg                                                             >  2:1993468/1‑99 (MQ=255)
                                        gctgctCAAAAAGCATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTgg                                                            <  2:2283052/100‑1 (MQ=255)
                                        gctgctCAAAAAGCATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTgg                                                            <  1:137412/100‑1 (MQ=255)
                                           gctCAAAAAGCATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACa                                                         >  1:924489/1‑100 (MQ=255)
                                           gctCAAAAAGCATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACa                                                         >  1:577753/1‑100 (MQ=255)
                                           gctCAAAAAGCATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACa                                                         >  2:2466206/1‑100 (MQ=255)
                                            ctCAAAAAGCATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACaa                                                        >  2:677973/1‑100 (MQ=255)
                                              cAAAAAGCATTAGAGCTTTTAGATTATGTATCTGTTTTAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGc                                                      <  1:150264/100‑1 (MQ=255)
                                              cAAAAAGCATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGc                                                      <  1:1047238/100‑1 (MQ=255)
                                              cAAAAAGCATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGc                                                      >  2:577528/1‑100 (MQ=255)
                                              cAAAAAGCATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGc                                                      >  2:2489547/1‑100 (MQ=255)
                                               aaaaaGCATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCa                                                     >  1:3511530/1‑100 (MQ=255)
                                                 aaaGCATTAGAGCTTTTAGATTATGTATCAGACATAAAAAAAAAAAAAAAAAAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCaa                                                    <  2:1761597/100‑1 (MQ=255)
                                                   aGCATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGGAAAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCaaaa                                                  >  1:2379655/1‑100 (MQ=255)
                                                    gCATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGGAAAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAg                                                 <  1:3221205/100‑1 (MQ=255)
                                                     cATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAGTa                                               <  2:3180926/100‑1 (MQ=255)
                                                      aTTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGGAAAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAGTa                                               <  1:1722823/100‑1 (MQ=255)
                                                       ttAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAGTAtc                                             >  1:2806050/1‑100 (MQ=255)
                                                       ttAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAGTAtc                                             <  1:3415855/100‑1 (MQ=255)
                                                        tAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAGTAtct                                            >  1:1669118/1‑100 (MQ=255)
                                                          gagCTTTTAGATTATGTATCAGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAGTATCTCa                                          <  2:4103777/100‑1 (MQ=255)
                                                             cTTTTAGATTATGTATCAGACATAACCACTGTCGATGGAAAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAGTATCTCata                                        >  2:2688968/1‑100 (MQ=255)
                                                              ttttAGATTATGTATCAGACATAACCACTGTCGATGGAAAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAGTATCTCatat                                       <  2:2305643/100‑1 (MQ=255)
                                                                 tAGATTATGTATCAGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAGTATCTCATATCtt                                    >  1:3918999/1‑99 (MQ=255)
                                                                 tAGATTATGTATCAGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAGTATCTCATATCTTa                                   <  1:756850/100‑1 (MQ=255)
                                                                 tAGATTATGTATCAGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAGTATCTCATATCTTa                                   >  1:3442281/1‑100 (MQ=255)
                                                                 tAGATTATGTATCAGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAGTATCTCATATCTTa                                   >  1:3539851/1‑100 (MQ=255)
                                                                  aGATTATGTATCAGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAGTATCTCATATCTTAt                                  >  2:2970298/1‑100 (MQ=255)
                                                                   gATTATGTATCAGACATAACCACTGTCGATGGAAAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAGTATCTCATATCTTAt                                  >  1:3106967/1‑100 (MQ=255)
                                                                    aTTATGTATCAGACATAACCACTGTCGATGGAAAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAGTATCTCATATCTTAtt                                 >  1:1202662/1‑100 (MQ=255)
                                                                    aTTATGTATCAGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAGTATCTCATATCTTAttt                                <  2:14738/100‑1 (MQ=255)
                                                                     ttATGTATCAGACATAACCACTGTCGATTGAAAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAGTATCTCATATCTTAttt                                <  2:844629/100‑1 (MQ=255)
                                                                     ttATGTATCAGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAGTATCTCATATCTTATTTg                               <  2:3539824/100‑1 (MQ=255)
                                                                      tATGTATCAGACATAACCACTGTCGATGGAAAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAGTATCTCATATCTTATTTg                               <  2:552947/100‑1 (MQ=255)
                                                                        tGTATCAGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAGTATCTCATATCTTATTTGAAc                            <  2:1837373/100‑1 (MQ=255)
                                                                         gTATCAGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAGTATCTCATATCTTATTTGAACa                           >  1:3034326/1‑100 (MQ=255)
                                                                           aTCAGACATAACCACTGTCGATGG‑ACAATTTTAGCTATTATGGGACAGATAACTGGTAAGTCTGGACACGCAAGAGGACCTAATAGCTTATTTGAAc                            >  2:3592676/1‑97 (MQ=255)
                                                                           aTCAGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAGTATCTCATATCTTATTTGAACAgg                         >  2:3337951/1‑100 (MQ=255)
                                                                            tCAGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAGTATCTCATATCTTATTTGAACAGGc                        <  2:2578465/100‑1 (MQ=255)
                                                                             cAGACATAACCACTGTCGATGGAAAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAGTATCTCATATCTTATTTGAACAGGc                        <  1:255870/100‑1 (MQ=255)
                                                                             cAGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAGTATCTCATATCTTATTTGAACAGGCt                       <  1:2857629/100‑1 (MQ=255)
                                                                              aGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAGTATCTCATATCTTATTTGAACAGGCTa                      <  1:2520671/100‑1 (MQ=255)
                                                                              aGACATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAGTATCTCATATCATATTTGAACAGGCTa                      >  1:4113907/1‑100 (MQ=255)
                                                                               gACATAACCACTGTCGATGGAAAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAGTATCTCATATCTTATTTGAACAGGCTa                      >  2:755800/1‑100 (MQ=255)
                                                                                aCATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAGTATCTCATATCTTATTTGAACAGGCTAAg                    <  2:2719359/100‑1 (MQ=255)
                                                                                aCATAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAGTATCTCATATCTTATTTGAACAGGCTAAg                    <  1:3695751/100‑1 (MQ=255)
                                                                                  aTAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAGTATCTCATATCTTATTTGAACAGGCTAAGAt                  <  1:1082214/100‑1 (MQ=255)
                                                                                   tAACCACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAGTATCTCATATCTTATTTGAACAGGCTAAGATa                 >  2:508427/1‑100 (MQ=255)
                                                                                      ccACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAGTATCTCATATCTTATTTGAACAGGCTAAGATacac              >  1:2496993/1‑100 (MQ=255)
                                                                                      ccACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAGTATCTCATATCTTATTTGAACAGGCTAAGATacac              <  1:1687658/100‑1 (MQ=255)
                                                                                      ccACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAGTATCTCATATCTTATTTGAACAGGCTAAGATacac              >  1:435046/1‑100 (MQ=255)
                                                                                      ccACTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAGTATCTCATATCTTATTTGAACAGGCTAAGATacac              <  1:1021515/100‑1 (MQ=255)
                                                                                        aCTGTCGATGGAAAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAGTATCTCATATCTTATTTGAACAGGCTAAGATACACt             <  2:295054/100‑1 (MQ=255)
                                                                                         cTGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAGTATCTCATATCTTATTTGAACAGGCTAAGATACACTCa           <  1:2692165/100‑1 (MQ=255)
                                                                                          tGTCGATGGAAAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAGTATCTCATATCTTATTTGAACAGGCTAAGATACACTCa           >  2:3229581/1‑100 (MQ=255)
                                                                                          tGTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAGTATCTCATATCTTATTTGAACAGGCTAAGATACACTCaa          >  1:825950/1‑100 (MQ=255)
                                                                                           gTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAGTATCTCATATCTTATTTGAACAGGCTAAGATACACTCAAc         <  1:1963416/100‑1 (MQ=255)
                                                                                           gTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAGTATCTCATATCTTATTTGAACAGGCTAAGATACACTCAAc         >  1:1265608/1‑100 (MQ=255)
                                                                                           gTCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAGTATCTCATATCTTATTTGAACAGGCTAAGATACACTCAAc         >  1:4036648/1‑100 (MQ=255)
                                                                                            tCGATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAGTATCTCATATCTTATTTGAACAGGCTAAGATACACTCAACt        >  1:4077611/1‑100 (MQ=255)
                                                                                             cGATGGAAAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAGTATCTCATATCTTATTTGAACAGGCTAAGATACACTCAACt        <  1:3168013/100‑1 (MQ=255)
                                                                                              gATGG‑AAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAGTATCTCATATCTTATTTGAACAGGCTAAGATACACTCAACTGa      <  2:1695506/100‑1 (MQ=255)
                                                                                                 ggAAAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAGTATCTCATATCTTATTTGAACAGGCTAAGATACACTCAACTGata    <  1:2475835/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
TATGTCTTTGGCACTTCATGGAAAAAGTGAACTTGAACTTGCTGCTCAAAAAGCATTAGAGCTTTTAGATTATGTATCAGACATAACCACTGTCGATGGAAAAATTTTAGCTATTATGGGACTGATAACTGGTCTGTCTGGACAAGCAAAAGTATCTCATATCTTATTTGAACAGGCTAAGATACACTCAACTGATATA  >  NC_000913/2993008‑2993206

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