Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,580,128 A→G 56.5% F221F (TTT→TTC ydeM ← putative YdeN‑specific sulfatase‑maturating enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,580,1280AG56.5% 43.9 / 130.1 111F221F (TTT→TTCydeMputative YdeN‑specific sulfatase‑maturating enzyme
Reads supporting (aligned to +/- strand):  ref base A (25/23);  new base G (26/36);  total (51/60)
Fisher's exact test for biased strand distribution p-value = 3.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.53e-01

Only 100 of 114 total aligned reads displayed.

TGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAG  >  NC_000913/1580030‑1580226
                                                                                                  |                                                                                                  
tGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATaa                                                                                                   <  2:3032535/100‑1 (MQ=255)
 gCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATaaa                                                                                                  >  2:2187648/1‑100 (MQ=255)
  cccATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAAc                                                                                                 >  2:1396535/1‑100 (MQ=255)
     aTTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTg                                                                                              <  1:1375552/100‑1 (MQ=255)
      ttCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGc                                                                                             <  1:1496104/100‑1 (MQ=255)
       tCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTggc                                                                                            <  2:3961402/100‑3 (MQ=255)
       tCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGcc                                                                                            <  1:746284/100‑1 (MQ=255)
         cccAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCAt                                                                                          >  1:1084997/1‑100 (MQ=255)
          ccAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATa                                                                                         >  1:253203/1‑100 (MQ=255)
           cAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATaa                                                                                        >  2:2125537/1‑100 (MQ=255)
             aaaaaCGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGc                                                                                      <  2:3791031/100‑1 (MQ=255)
              aaaaCGGCTGACAAAACTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTCGACATAAACTTGCCATAAGcc                                                                                     <  2:3651560/100‑1 (MQ=255)
               aaaCGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCg                                                                                    >  1:2465698/1‑100 (MQ=255)
                  cGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTg                                                                                  >  1:2941225/1‑99 (MQ=255)
                   ggCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTggg                                                                                >  2:3293017/1‑100 (MQ=255)
                   ggCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTggg                                                                                >  1:1694498/1‑100 (MQ=255)
                   ggCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTggg                                                                                <  1:1467730/100‑1 (MQ=255)
                    gCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGa                                                                               >  1:3962486/1‑100 (MQ=255)
                    gCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGa                                                                               <  2:3963128/100‑1 (MQ=255)
                       gACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGc                                                                            >  1:1241205/1‑100 (MQ=255)
                         cAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGcac                                                                          <  2:2756559/100‑1 (MQ=255)
                          aaaGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGcaca                                                                         <  2:3210953/100‑1 (MQ=255)
                           aaGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAg                                                                        >  1:1915450/1‑100 (MQ=255)
                           aaGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAg                                                                        >  2:3831818/1‑100 (MQ=255)
                            aGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGa                                                                       >  2:1865563/1‑100 (MQ=255)
                              cTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGaaa                                                                     >  2:128390/1‑100 (MQ=255)
                               tttCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGaaaa                                                                    >  2:3121048/1‑100 (MQ=255)
                                 tCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAAt                                                                  <  1:2462174/100‑1 (MQ=255)
                                  cAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAAtn                                                                 <  2:2975449/100‑2 (MQ=255)
                                  cAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATc                                                                 <  2:4175750/100‑1 (MQ=255)
                                   aaaCTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCa                                                                <  1:1989375/100‑1 (MQ=255)
                                   aaaCTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCa                                                                <  1:2818576/100‑1 (MQ=255)
                                   aaaCTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCa                                                                <  1:18168/100‑1 (MQ=255)
                                   aaaCTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCa                                                                <  2:1451078/100‑1 (MQ=255)
                                   aaaCTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCa                                                                <  2:571047/100‑1 (MQ=255)
                                   aaaCTGACGGATGAAAATTTCACCCACATTGCTTTTAAACCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCa                                                                >  1:1005186/1‑100 (MQ=255)
                                     aCTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAAt                                                              <  1:622790/100‑1 (MQ=255)
                                      cTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATg                                                             <  2:2941237/100‑1 (MQ=255)
                                       tGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGa                                                            >  1:165960/1‑100 (MQ=255)
                                        gACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGa                                                            >  1:3692959/1‑99 (MQ=255)
                                        gACAGATGAAAATTTCAGCCACATTGTTTTAAACCCATGGTATAAGATGGGTTGACATGACCTTGCCGTAAGGCAGTGAAGGCACAGAAAAATCAAt                                                              <  1:395604/97‑1 (MQ=255)
                                         aCGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATg                                                             <  1:1528269/97‑1 (MQ=255)
                                            gATGAAAATTTCACCCACATTGTTTTGAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACgg                                                       <  1:2991416/100‑1 (MQ=255)
                                             aTGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGa                                                      >  2:3654535/1‑100 (MQ=255)
                                              tGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGaa                                                     >  1:2865233/1‑100 (MQ=255)
                                              tGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGaa                                                     >  1:3299221/1‑100 (MQ=255)
                                               gAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAAt                                                    <  1:3111239/100‑1 (MQ=255)
                                               gAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAAt                                                    >  2:2883900/1‑100 (MQ=255)
                                               gAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAAt                                                    >  2:2796904/1‑100 (MQ=255)
                                               gAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAAt                                                    <  1:4228345/100‑1 (MQ=255)
                                                aaaaTTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAAtg                                                   <  2:3896075/100‑1 (MQ=255)
                                                 aaaTTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAAtgt                                                  >  2:3216512/1‑100 (MQ=255)
                                                 aaaTTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAAtgt                                                  >  2:3236127/1‑100 (MQ=255)
                                                  aaTTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAAtgtg                                                 >  2:583302/1‑100 (MQ=255)
                                                    tttCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGtt                                               >  2:1386429/1‑100 (MQ=255)
                                                     ttCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTtc                                              >  1:446300/1‑100 (MQ=255)
                                                     ttCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTtc                                              <  1:1996265/100‑1 (MQ=255)
                                                       cACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTtctc                                            <  1:129083/100‑1 (MQ=255)
                                                        aCCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCa                                           >  1:995383/1‑100 (MQ=255)
                                                          ccACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACt                                         <  2:3930693/100‑1 (MQ=255)
                                                          ccACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGAAGAAACAAAAAAATCAATGATACGGAATGTGTTCTCACt                                         >  1:2948534/1‑100 (MQ=255)
                                                           cacaTTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTa                                        <  1:3831801/100‑1 (MQ=255)
                                                            acaTTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTAt                                       <  1:750231/100‑1 (MQ=255)
                                                            acaTTGTTTTTAACCCATTGCATAAAAAAGGTTGACATGAACTTGCCATAAGCCGTGAGAGAAAAAAAAAAATCAATGATACGGAATGTGTTCTCACTa                                        <  1:756921/99‑1 (MQ=255)
                                                             cattgTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATg                                      <  1:3874639/96‑1 (MQ=255)
                                                             cattgTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATg                                      <  2:481224/96‑1 (MQ=255)
                                                             caTTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATg                                      <  1:3029157/100‑1 (MQ=255)
                                                               ttgTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGAc                                    >  1:2257466/3‑100 (MQ=255)
                                                               ttgTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGAc                                    <  1:3010512/98‑1 (MQ=255)
                                                                 gTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCa                                  >  1:1110507/1‑100 (MQ=255)
                                                                  tttttAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCAc                                 <  1:867143/100‑1 (MQ=255)
                                                                   ttgtAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGAAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACt                                <  1:1328292/97‑1 (MQ=255)
                                                                    tttAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTg                               >  2:2062922/1‑100 (MQ=255)
                                                                     ttAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGa                              <  1:3445387/100‑1 (MQ=255)
                                                                        aCCCATTGCATAAAAATGGTTGACATTAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGaaa                            <  1:3764313/99‑1 (MQ=255)
                                                                        aCCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAAt                           <  1:4053518/100‑1 (MQ=255)
                                                                         cccATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATc                          <  2:691305/100‑1 (MQ=255)
                                                                         cccATTGCATAAAAACGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATc                          >  1:1568549/1‑100 (MQ=255)
                                                                          ccATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAAtga                         >  1:2586623/1‑98 (MQ=255)
                                                                           cATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCaa                        <  2:3642114/100‑1 (MQ=255)
                                                                           cATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCaa                        <  1:583300/100‑1 (MQ=255)
                                                                           cATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCaa                        <  1:1036171/100‑1 (MQ=255)
                                                                           cATTGCATAAAAATAGTTAACATGAACTTGCCATAAGCCGTGAAAAACAAAAAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCa                         <  1:4182450/99‑1 (MQ=255)
                                                                            aTTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAat                       <  2:382344/100‑1 (MQ=255)
                                                                              tGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATGCGGAATGTGTTCTCACTATGACCACTGAAATCAata                      <  2:3011442/99‑1 (MQ=255)
                                                                               gCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATAtt                    >  2:4178314/1‑100 (MQ=255)
                                                                               gCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATAtt                    >  2:1903774/1‑100 (MQ=255)
                                                                                cATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTc                   >  2:3305025/1‑100 (MQ=255)
                                                                                 aTAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCg                  >  2:2218415/1‑100 (MQ=255)
                                                                                 aTAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCg                  <  1:2074975/100‑1 (MQ=255)
                                                                                  tAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCgg                 >  2:549472/1‑100 (MQ=255)
                                                                                    aaaaTGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCg               <  1:100324/100‑1 (MQ=255)
                                                                                     aaaTGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTGTCACTATGACCACTGAAATCAATATTCGGc                >  1:313096/1‑98 (MQ=255)
                                                                                         ggTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCg          >  1:2381906/1‑100 (MQ=255)
                                                                                          gTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCgg         >  2:3485504/1‑100 (MQ=255)
                                                                                          gTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCgg         >  2:3540895/1‑100 (MQ=255)
                                                                                            tGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGtt       <  1:2445420/100‑1 (MQ=255)
                                                                                              aCATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTc     >  2:2724805/1‑100 (MQ=255)
                                                                                              aCATAAACATGACATAATCCGAGGGAGGCACATAAAAAACAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAAGATTCTGCGTCTCGGttt      <  1:196659/99‑1 (MQ=255)
                                                                                                  aaaCTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAg  <  1:1791331/99‑1 (MQ=255)
                                                                                                  |                                                                                                  
TGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAG  >  NC_000913/1580030‑1580226

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