Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 16,923 T→G 7.2% K13T (AAG→ACG)  mokC ← regulatory protein MokC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091316,9230TG 7.2% 161.4 / 4.6 70K13T (AAG→ACG) mokCregulatory protein MokC
Reads supporting (aligned to +/- strand):  ref base T (25/38);  new base G (3/4);  total (28/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.00e-01

CCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATA  >  NC_000913/16783‑17065
                                                                                                                                            |                                                                                                                                              
ccTCCGTCTGGCCAGTCCGGATGTGAACCTCACAGAGGTCTTTTCTCGGTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTa                                                                                                                                              <  2:417636/143‑1 (MQ=255)
      cccggCCAGTGCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCa                                                                                                                                        <  1:61692/140‑1 (MQ=255)
         ggCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCGTAAGGTCaagg                                                                                                                                     <  1:445352/143‑4 (MQ=255)
            cAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAaga                                                                                                                                  >  1:951633/1‑143 (MQ=255)
              gTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAgg                                                                                                                                <  2:666974/143‑1 (MQ=255)
              gTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAgg                                                                                                                                <  2:575069/143‑1 (MQ=255)
                tCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTCTGCTGCTTCATTGCTCTCATCTCCTTGACCTTACGGTCACTAAGCGGCa                                                                                                                              >  1:545090/1‑143 (MQ=255)
                 cGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCAc                                                                                                                             <  2:667004/143‑1 (MQ=255)
                  ggATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGCTCAGTAAGAAGCAct                                                                                                                            >  1:938981/1‑143 (MQ=255)
                    aTGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGGTCAGGGCGACAATCATCGCCTTATGCGGCGTCATTGCGCTCCTCTCCTTGGCCTTCCGGTCAGTACGCAGCActct                                                                                                                          >  1:251347/1‑143 (MQ=255)
                     tgtgAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTa                                                                                                                         >  2:108057/1‑143 (MQ=255)
                        gAACCTCACAGGGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACAt                                                                                                                      <  2:291820/143‑1 (MQ=255)
                        gAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACAt                                                                                                                      <  1:654255/143‑1 (MQ=255)
                          aCCGCACAGAGGGCTTGTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACAtgt                                                                                                                    <  1:61076/143‑1 (MQ=255)
                           ccTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGCTGCCGATCAGGGAGGAAGTCGTCGCCCTATGCTGCTTCATTGCTCTCTTCTCCTTCACCTTACGGTCAGTCAGAGGCACTATAAAtgtg                                                                                                                   >  2:133002/1‑143 (MQ=255)
                               acaGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAAGCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCa                                                                                                               <  1:1007418/143‑1 (MQ=255)
                                     gTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGACGCACTCTACATGTGTgcagcacac                                                                                                         >  2:220675/1‑140 (MQ=255)
                                         tttCTCGTGTCCGGCGCCCCCACTTCGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATagga                                                                                                     <  2:380781/143‑1 (MQ=255)
                                                  aCCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGCGTTCAGCATATAGGAGGCCTCggg                                                                                            >  2:709025/1‑143 (MQ=255)
                                                     aGCGCCGCCACTTCTGTGGTGTTACAGATGACGATCACGGCGACAATCATCGCCTTATGCTGCATCATTGCTCTCCTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTg                                                                                         <  1:268275/143‑1 (MQ=255)
                                                     aGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTg                                                                                         >  2:725740/1‑143 (MQ=255)
                                                      gcgcCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCGTCGGGTTGa                                                                                        >  1:86768/1‑143 (MQ=255)
                                                       cgccgcCCCCACGGCGGTGATACAGATGACGACCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTCAGAGGCACCGTACATGTCCTCAGCATATAGGAGGCCTCGGGTTGAt                                                                                       <  1:520473/143‑1 (MQ=255)
                                                       cgccgcCACGTCGGGGGTTATGCAGAGGACGATCAGGGGGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGAt                                                                                       <  2:655413/143‑1 (MQ=255)
                                                              aCTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATACGACGCCTCGGGTCGATTGTaaaa                                                                                >  2:511381/1‑143 (MQ=255)
                                                              aCTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCCGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTaaaa                                                                                <  1:956437/143‑1 (MQ=255)
                                                               cTACGGCGGTGATACCGAGGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAat                                                                               <  2:231832/143‑1 (MQ=255)
                                                               cTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAat                                                                               <  1:995067/143‑1 (MQ=255)
                                                                       gTGATACAGATGACGCTCAGTGCTACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCg                                                                       <  2:433213/143‑1 (MQ=255)
                                                                         gATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCggg                                                                     <  2:293866/143‑1 (MQ=255)
                                                                          atacttgtGACGATCAGTGCGACCATCATCGCCTTCTGCTGCTTCATTGCCCTCTTCTCCTTGACCTTACGGTCCGTAATAGGCACTCTACATGTGTTCGGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCgggg                                                                    <  2:468974/136‑1 (MQ=255)
                                                                          aTACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCCTATAGGAGGCCTCGGGTTGATGGTAACATATCACTCgggg                                                                    >  2:14635/1‑143 (MQ=255)
                                                                                ctGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTtc                                                              <  2:414774/142‑1 (MQ=255)
                                                                                aTGACGATCAGGGCGACAATCATCGCCTTATGCCTCTTCGTTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTtc                                                              >  2:447933/1‑143 (MQ=255)
                                                                                  gACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTtctc                                                            >  2:36019/1‑143 (MQ=255)
                                                                                    cGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTa                                                          <  1:133002/143‑1 (MQ=255)
                                                                                       tCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGCGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCt                                                       >  1:529592/1‑143 (MQ=255)
                                                                                         aGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGc                                                     <  1:587235/143‑1 (MQ=255)
                                                                                           ggCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTCCTCCTTGACCTGACGGTCAGTAAGAGGCACTCTACAGGTGGTCAGCAGATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCg                                                   <  1:1078584/143‑1 (MQ=255)
                                                                                           ggCGACAATCATCGCCTTATGCTGCTTCATTGCGCTCTTCTCCTTGTCCTTTTGTTCGGTCAGCGGCACCCTACATGGGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCg                                                   <  2:491202/143‑1 (MQ=255)
                                                                                            gCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGt                                                  <  2:81124/143‑1 (MQ=255)
                                                                                             cGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCGGGGATAGGCACACCACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTACAATATCACTCGGGGCTTTTCTCTATCTGCCGtt                                                 >  1:466515/1‑143 (MQ=255)
                                                                                              gACAACCATCGCCTTATGCTGCTTCCTTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTc                                                <  1:391811/143‑1 (MQ=255)
                                                                                               aCAATCCTCGCCTTATGCTGTTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGGTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCa                                               <  1:438429/143‑1 (MQ=255)
                                                                                                       cGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGc                                       >  2:501606/1‑143 (MQ=255)
                                                                                                         ccTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTACTCTATCTGCCGTTCAGCTAATGCCt                                     >  2:1011860/1‑143 (MQ=255)
                                                                                                           ttATGCTGGTTCTTTTCTCCCTTCTCCGTGGCCGTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTga                                   <  2:951633/143‑1 (MQ=255)
                                                                                                           ttATGCTGCTTCATTGCTCTCTTCTCCTTTACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTga                                   <  2:1089951/143‑1 (MQ=255)
                                                                                                           ttATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTga                                   <  1:671779/143‑1 (MQ=255)
                                                                                                           ttATGCTGCTTCATTCCCCTCTTCGCCTTGACCTTACGGTCAGTAAGAGGCACTCCACACGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTga                                   <  1:98571/143‑1 (MQ=255)
                                                                                                            tATGCTGCTTCATTGCTCTCTTCTCCTGGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTgag                                  <  1:943422/143‑1 (MQ=255)
                                                                                                             acgctgctTCAGTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTgaga                                 <  1:366879/141‑1 (MQ=255)
                                                                                                               gctgctTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGagaca                               >  2:202581/1‑143 (MQ=255)
                                                                                                                 tgctTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTCCATGTGTTCAGCATATAGGAGACCTCGGGTTGATGGTAAAATATCACTCGGGGCATTTCTCTATCTGCCGTTCAGCTAATGCCTGtgacagg                             >  1:70168/1‑142 (MQ=255)
                                                                                                                 tgctTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCACAGGACTTTTTGAGAATATCACTCGGTGCTTTTCTCTATCTGCCGTTCAGCTAATGCTTGagacagg                             >  1:808063/1‑142 (MQ=255)
                                                                                                                  gctTCATTGCTCGCTTCTCCGTGGCCGTGCGGTCAGTAAGGGGCGCTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGagacagac                            <  2:227057/143‑1 (MQ=255)
                                                                                                                      cATTGCTCTCTTCTCCTTGGCCTTACGGTCAGTAAGAGGCACTCTACACGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGcc                        <  1:95181/143‑1 (MQ=255)
                                                                                                                        ttGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTc                      >  2:399147/1‑143 (MQ=255)
                                                                                                                          cttctctTCTCCTTGACCTGGCGGTCAGTAAGAGGCACTCTACTTGTTTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAGTATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCaa                    <  2:925209/141‑1 (MQ=255)
                                                                                                                           ctctctTCTCCTTGACCTTGCGGTCAGTAGGAGGCACTCTACATGGGGTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCGAg                   <  1:950038/143‑1 (MQ=255)
                                                                                                                           ctctctTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAg                   <  1:23430/143‑1 (MQ=255)
                                                                                                                           ctctctTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAg                   >  1:603084/1‑143 (MQ=255)
                                                                                                                           cgctctTCTCCTTGACCGTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAg                   <  1:952352/141‑1 (MQ=255)
                                                                                                                             ctctTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAat                                                                               <  1:136329/81‑1 (MQ=255)
                                                                                                                             ctctTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAat                                                                               >  2:136329/1‑81 (MQ=255)
                                                                                                                                ttctCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGCGGCCTCGGGATGATGGTACAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCCGCTAATGTCTGAGACAGACAGCCTCCAGCAccc              >  2:787814/1‑143 (MQ=255)
                                                                                                                                ttctCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCAccc              <  1:987223/143‑1 (MQ=255)
                                                                                                                                ctcttactGACCGTAAGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCAccc              >  1:666067/8‑143 (MQ=255)
                                                                                                                                 tcttactGACCGTAAGGTCAGTAAGAGGCACTCTACATGTGTTCAGCGTATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTGTCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACccg             >  1:192780/7‑143 (MQ=38)
                                                                                                                                 tcttactGACCGTAAGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACccg             >  2:874042/7‑143 (MQ=255)
                                                                                                                                            ttACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATtata  >  1:65028/1‑143 (MQ=255)
                                                                                                                                            |                                                                                                                                              
CCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATA  >  NC_000913/16783‑17065

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: ‑