Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,580,128 A→G 21.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,1280AG21.5% 210.1 / 55.4 111F221F (TTT→TTCydeMputative YdeN‑specific sulfatase‑maturating enzyme
Reads supporting (aligned to +/- strand):  ref base A (40/47);  new base G (12/12);  total (52/59)
Fisher's exact test for biased strand distribution p-value = 8.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 113 total aligned reads displayed.

TGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGC  >  NC_000913/1580030‑1580227
                                                                                                  |                                                                                                   
tGCCCATTCCCCACAAAACGGCTGACAAAGCTTTCACACTGCCGGATGAAAATTTCACCCACGTTGTTTTTAACCCATTGCATAAACATGGTTGACATaa                                                                                                    >  2:3818138/1‑100 (MQ=255)
  cccATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAAc                                                                                                  >  1:2729861/1‑100 (MQ=255)
   ccATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACt                                                                                                 >  1:715835/1‑100 (MQ=255)
   ccATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACt                                                                                                 <  2:2283407/100‑1 (MQ=255)
   ccATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTCACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACt                                                                                                 <  2:2805150/100‑1 (MQ=255)
    cATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACtt                                                                                                >  1:3565130/1‑100 (MQ=255)
    cATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACtt                                                                                                <  1:2368642/100‑1 (MQ=255)
     aTTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTg                                                                                               <  1:1242278/100‑1 (MQ=255)
     aTTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTg                                                                                               <  1:2672870/100‑1 (MQ=255)
     aTTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTg                                                                                               <  2:3934157/100‑1 (MQ=255)
      ttCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGc                                                                                              <  2:1770429/100‑1 (MQ=255)
       tCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGcc                                                                                             <  1:3694681/100‑1 (MQ=255)
       tCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGcc                                                                                             <  2:1668141/100‑1 (MQ=255)
        ccccAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCa                                                                                            >  2:2509626/1‑100 (MQ=255)
        ccccAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCa                                                                                            >  2:188342/1‑100 (MQ=255)
         cccAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATACAAATGGTTGACATAAAATTGCCAt                                                                                           >  1:3078173/1‑100 (MQ=255)
         cccAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCAt                                                                                           >  1:3339039/1‑100 (MQ=255)
          ccAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATa                                                                                          >  1:2465664/1‑100 (MQ=255)
          ccAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATa                                                                                          >  1:471864/1‑100 (MQ=255)
          ccAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATa                                                                                          >  1:135481/1‑100 (MQ=255)
           cAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCaaga                                                                                         >  1:4097861/1‑97 (MQ=255)
           cAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATaa                                                                                         >  2:1000452/1‑100 (MQ=255)
             aaaaaCGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGc                                                                                       <  2:2649407/100‑1 (MQ=255)
             aaaaaCGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGCTGACATAAACTTGCCATAAGc                                                                                       >  2:2847688/1‑100 (MQ=255)
                  cGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTgg                                                                                  <  1:128478/100‑1 (MQ=255)
                   ggCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTggg                                                                                 >  1:1289417/1‑100 (MQ=255)
                    gCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGa                                                                                <  1:33866/100‑1 (MQ=255)
                     cTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAg                                                                               <  1:2188909/100‑1 (MQ=255)
                      tGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAgg                                                                              <  1:3591353/100‑1 (MQ=255)
                      tGACAAAGCTTTCAAACTGACGGAGGAAAATTTCACCCACATTGGTCTTAACCCATTGCGTAAAAATGGTTTACATAAACTTGCCATAGGCCGTGGGAgg                                                                              <  2:1216360/100‑1 (MQ=255)
                       gACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGc                                                                             <  1:210167/100‑1 (MQ=255)
                       gACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGc                                                                             <  1:2509641/100‑1 (MQ=255)
                       gACAAAGCTTTCAAACTGACGGATGAAAAATTCCCACACATTGTTTTTAACCCATTGCTTAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGc                                                                             <  2:877321/100‑1 (MQ=255)
                        aCAAAGCTTTCAAACTGACGGATGAAAATTTTATCCACATTGTTTTTAAAAAAAAAAAAAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGc                                                                             <  2:2319353/99‑1 (MQ=255)
                        aCAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGca                                                                            <  2:1022756/100‑1 (MQ=255)
                         cAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGcac                                                                           >  1:419443/1‑100 (MQ=255)
                         cAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGcac                                                                           >  1:962810/1‑100 (MQ=255)
                           aaGCTTTCAAACTGACGGATGAAAATTTCACTCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAg                                                                         >  1:2813713/1‑100 (MQ=255)
                           aaGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAg                                                                         <  2:2041895/100‑1 (MQ=255)
                           aaGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAg                                                                         >  1:2107611/1‑100 (MQ=255)
                           aaGCTTTCAAACTAACGGATGAAAATTTCACGCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAg                                                                         >  1:2813667/1‑100 (MQ=255)
                              cTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGaaa                                                                      >  1:2255556/1‑100 (MQ=255)
                              cTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGaaa                                                                      >  2:457966/1‑100 (MQ=255)
                                   aaaCTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCa                                                                 <  2:1599438/100‑1 (MQ=255)
                                     aCTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAAt                                                               <  1:2003431/100‑1 (MQ=255)
                                     aCTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAAt                                                               <  1:367654/100‑1 (MQ=255)
                                      cTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATg                                                              <  2:2936183/100‑1 (MQ=255)
                                           aaataaaaaTTTCTCGTACATTTTGTTTAACCCATTGCATAAAAATGGCTAACATAAACTTGCCATAAGCCGTGAGTTTCACAGAAAAATCAATGacac                                                          <  2:4085912/94‑1 (MQ=255)
                                              tGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGaa                                                      >  1:2379546/1‑100 (MQ=255)
                                               gAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGaa                                                      >  2:3463612/1‑99 (MQ=255)
                                                aaaaTTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGAGAGGCACAGAAAAATCAATGATACGGAAtg                                                    >  2:3108142/1‑100 (MQ=255)
                                                     ttCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTtc                                               <  2:1130236/100‑1 (MQ=255)
                                                     atCACCAACATTCTTTTTAACAAAATGAATAAAAACGGATGACATAAACATGCAATCAGACGTGAAAAGCACAGAAAAATCAATGATACGGaaagaaaac                                               <  2:3749303/99‑8 (MQ=255)
                                                        aCCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCa                                            <  2:1045542/100‑1 (MQ=255)
                                                        aCCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCa                                            <  2:1045760/100‑1 (MQ=255)
                                                        aCCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCa                                            >  2:2121009/1‑100 (MQ=255)
                                                        aCCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCa                                            <  1:1957417/100‑1 (MQ=255)
                                                         cccACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCAc                                           >  1:783226/1‑100 (MQ=255)
                                                           cacaTTGTTTTTAACCCATTGCATACAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATCCGGAATGTGTTCTCACTa                                         >  2:3869361/1‑100 (MQ=255)
                                                           cacaTTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTa                                         <  2:155528/100‑1 (MQ=255)
                                                           cacaTTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTa                                         >  1:473578/1‑100 (MQ=255)
                                                           cacaTTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTa                                         <  1:1979905/100‑1 (MQ=255)
                                                            acaTTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTAt                                        >  1:2485602/1‑100 (MQ=255)
                                                             cattgTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATg                                       <  2:2729845/96‑1 (MQ=255)
                                                              attgTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGa                                      >  1:1764353/4‑100 (MQ=255)
                                                              attgTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGa                                      >  2:3757729/4‑100 (MQ=255)
                                                               ttgTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGAc                                     <  1:188345/98‑1 (MQ=255)
                                                                tgTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGAcc                                    <  2:1542776/99‑1 (MQ=255)
                                                                tgTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGAcc                                    <  1:3087683/99‑1 (MQ=255)
                                                                tgTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGAcc                                    <  2:3524826/99‑1 (MQ=255)
                                                                  tttttAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCAc                                  >  2:774632/1‑100 (MQ=255)
                                                                   ttttAACCCATTGCATAAAAATGGTTGACAAGAACTGGCCATAAGCGGGGGGAGGCACAGAAAAATCAATGATACGGCAGGTGTGCTCACTATGACCACt                                 >  2:1385209/1‑100 (MQ=255)
                                                                    tttAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTg                                <  1:2487867/100‑1 (MQ=255)
                                                                    tttAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTg                                <  1:3867156/100‑1 (MQ=255)
                                                                        aCCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAAt                            <  1:3164045/100‑1 (MQ=255)
                                                                         cccATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATc                           >  1:3766648/1‑100 (MQ=255)
                                                                         cccATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATc                           <  1:1291001/100‑1 (MQ=255)
                                                                           cATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCaa                         <  2:1994114/100‑1 (MQ=255)
                                                                           cATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCaa                         <  2:680322/100‑1 (MQ=255)
                                                                               gCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATAtt                     >  1:2953528/1‑100 (MQ=255)
                                                                                cATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTc                    >  1:2994124/1‑100 (MQ=255)
                                                                                  tAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCgg                  >  1:2756799/1‑100 (MQ=255)
                                                                                   aaaaaTGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGc                 <  2:1151059/100‑1 (MQ=255)
                                                                                    aaaaTGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCg                >  1:173253/1‑100 (MQ=255)
                                                                                    aaaaTGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCg                >  1:524085/1‑100 (MQ=255)
                                                                                    aaaaTGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCg                >  1:2786382/1‑100 (MQ=255)
                                                                                     aaaTGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGc                 >  2:553927/1‑98 (MQ=255)
                                                                                     aaaTGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGt               >  2:2651583/1‑100 (MQ=255)
                                                                                        tGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTccc            >  1:145398/1‑100 (MQ=255)
                                                                                         ggTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCg           >  1:2234278/1‑100 (MQ=255)
                                                                                          gTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGAGTCCCgg          >  1:849071/1‑100 (MQ=255)
                                                                                          gTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCgg          <  1:1473998/100‑1 (MQ=255)
                                                                                          gTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCgg          >  2:1870455/1‑100 (MQ=255)
                                                                                          gTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCgg          <  1:585154/100‑1 (MQ=255)
                                                                                            tGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGtt        <  2:1395382/100‑1 (MQ=255)
                                                                                             gACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGttt       <  1:110779/100‑1 (MQ=255)
                                                                                              aCATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTc      <  1:2988081/100‑1 (MQ=255)
                                                                                              aCATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTc      <  2:595237/100‑1 (MQ=255)
                                                                                                 tAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAg   <  2:133004/100‑1 (MQ=255)
                                                                                                  aaaCTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGc  <  2:456012/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
TGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGC  >  NC_000913/1580030‑1580227

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