Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,220,971 G→A pseudogene (1371/2648 nt) ycgH → putative transporter component, N‑terminal fragment

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,220,9710GA100.0% 249.7 / NA 69pseudogene (1371/2648 nt)ycgHputative transporter component, N‑terminal fragment
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (36/33);  total (36/33)

TTCAGACAGTATTATCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCGGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGAATGGCGATGAATCTCCAGCAG  >  NC_000913/1220836‑1221107
                                                                                                                                       |                                                                                                                                         
ttCAGACAGTATTATCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTa                                                                                                                                       <  2:1041573/140‑1 (MQ=255)
   aGACAGTATTATCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACgg                                                                                                                                   <  2:97215/140‑1 (MQ=255)
     aCAGTATTATCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGAc                                                                                                                                 <  2:705314/140‑1 (MQ=255)
     aCAGTATTATCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGAc                                                                                                                                 <  2:714333/140‑1 (MQ=255)
         tattatCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAATAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGAc                                                                                                                             >  2:425971/1‑140 (MQ=255)
         tattatCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGAc                                                                                                                             >  1:501822/1‑140 (MQ=255)
         tattatCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGAc                                                                                                                             >  1:581486/1‑140 (MQ=255)
           ttatCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACtt                                                                                                                           >  1:868447/1‑140 (MQ=255)
            tatCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTa                                                                                                                          <  2:132213/140‑1 (MQ=255)
                 aTCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTAttgtt                                                                                                                     >  2:852137/1‑140 (MQ=255)
                 aTCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTAttgtt                                                                                                                     <  1:227424/140‑1 (MQ=255)
                 aTCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTAttgtt                                                                                                                     <  1:339710/140‑1 (MQ=255)
                   cTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAGTGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGa                                                                                                                   >  2:683217/1‑140 (MQ=255)
                   cTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAGTCAGCTACGGAc                                                                                                                                 <  2:93566/126‑1 (MQ=255)
                   cTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGAc                                                                                                                                 >  1:93566/1‑126 (MQ=255)
                   cTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGa                                                                                                                   <  1:172569/140‑1 (MQ=255)
                   cTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGa                                                                                                                   >  2:637417/1‑140 (MQ=255)
                     ttCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAgg                                                                                                                 >  2:533033/1‑140 (MQ=255)
                            gatgatgGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATAc                                                                                                          >  1:261937/1‑140 (MQ=255)
                                atgGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCg                                                                                                      <  2:162806/140‑1 (MQ=255)
                                atgGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCg                                                                                                      >  2:109861/1‑140 (MQ=255)
                                atgGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCg                                                                                                      >  1:843020/1‑140 (MQ=255)
                                 tgGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCgg                                                                                                     <  2:868447/140‑1 (MQ=255)
                                  ggTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCggg                                                                                                    <  2:665479/140‑1 (MQ=255)
                                  ggTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCggg                                                                                                    <  2:969894/140‑1 (MQ=255)
                                    tCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGtt                                                                                                  <  1:1039771/140‑1 (MQ=255)
                                    tCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGtt                                                                                                  <  2:738529/140‑1 (MQ=255)
                                       aaCATTTACACAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCaa                                                                                               <  2:640759/140‑1 (MQ=255)
                                             tACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCg                                                                                         >  2:548923/1‑140 (MQ=255)
                                              aCCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGt                                                                                        <  2:261937/140‑1 (MQ=255)
                                              aCCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGt                                                                                        <  1:109861/140‑1 (MQ=255)
                                                               aaCCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTg                                                                       >  1:216789/1‑140 (MQ=255)
                                                                          gTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAgcggc                                                            <  1:533033/140‑1 (MQ=255)
                                                                            aaTTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAgcggcgc                                                          >  1:630296/1‑140 (MQ=255)
                                                                              ttATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAgcggcgcgg                                                        >  1:98927/1‑140 (MQ=255)
                                                                              ttATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAgcggcgcgg                                                        <  1:637417/140‑1 (MQ=255)
                                                                               tATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAgcggcgcggc                                                       <  1:800657/140‑1 (MQ=255)
                                                                                   tCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACg                                                   >  1:95390/1‑140 (MQ=255)
                                                                                   tCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACg                                                   >  1:669661/1‑140 (MQ=255)
                                                                                   tCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACg                                                   <  2:581486/140‑1 (MQ=255)
                                                                                    cGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGc                                                  >  1:810261/1‑140 (MQ=255)
                                                                                               gTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATc                                       <  2:630296/140‑1 (MQ=255)
                                                                                               gTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATc                                       <  1:852137/140‑1 (MQ=255)
                                                                                               gTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATc                                       <  2:216789/140‑1 (MQ=255)
                                                                                                 cAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGa                                     >  1:316706/1‑140 (MQ=255)
                                                                                                  aGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGaa                                    >  1:1004207/1‑140 (MQ=255)
                                                                                                    cTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACt                                  <  1:105199/140‑1 (MQ=255)
                                                                                                    cTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACt                                  <  2:316706/140‑1 (MQ=255)
                                                                                                    cTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACt                                  >  1:954116/1‑140 (MQ=255)
                                                                                                    cTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCAAGCTAAACGGGATCGAACt                                  <  2:98927/140‑1 (MQ=255)
                                                                                                      tAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCa                                <  2:458011/140‑1 (MQ=255)
                                                                                                            ccGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGCGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACtg                          >  1:941278/1‑140 (MQ=255)
                                                                                                               cACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGa                       >  1:51246/1‑140 (MQ=255)
                                                                                                                  cGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGAATg                    <  2:954116/140‑1 (MQ=255)
                                                                                                                    tATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGAATGGc                  <  2:117768/140‑1 (MQ=255)
                                                                                                                         ggTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGAATGGCGATGa             <  2:810261/140‑1 (MQ=255)
                                                                                                                          gTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGAATGGCGATGaa            <  2:95390/140‑1 (MQ=255)
                                                                                                                            gatgatAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGAATGGCGATGAAtc          >  2:231609/1‑140 (MQ=255)
                                                                                                                            gatgatAAATCAGCCACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGAATGGCGATGAAtc          >  2:733297/1‑140 (MQ=255)
                                                                                                                             atgatAAATCAGCTACGGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGAATGGCGATGAAtc          >  2:990472/1‑140 (MQ=255)
                                                                                                                             atgatAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGAATGGCGATGAAtct         >  2:1006624/1‑140 (MQ=255)
                                                                                                                              tgatAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGAATGGCGATGa             >  2:690663/1‑135 (MQ=255)
                                                                                                                              tgatAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGAATGGCGATGa             <  1:690663/135‑1 (MQ=255)
                                                                                                                              tgatAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGAATGGCGATGAAtctc        >  2:707458/1‑140 (MQ=255)
                                                                                                                              tgatAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGAATGGCGATGAAtctc        >  2:492182/1‑140 (MQ=255)
                                                                                                                              tgatAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGAATGGCGATGAAtctc        >  2:407502/1‑140 (MQ=255)
                                                                                                                              tgatAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGAATGGCGATGAAtctc        >  1:659382/1‑140 (MQ=255)
                                                                                                                              tgatAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGAATGGCGATGAAtctc        >  1:350462/1‑140 (MQ=255)
                                                                                                                              tgatAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGAATGGCGATGAAtctc        >  1:248627/1‑140 (MQ=255)
                                                                                                                                atAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGAATGGCGATGAATCTCca      >  1:1065951/1‑140 (MQ=255)
                                                                                                                                    aTCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGAATGGCGATGAATCTCcagcag  <  2:53385/140‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                         
TTCAGACAGTATTATCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCGGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGAATGGCGATGAATCTCCAGCAG  >  NC_000913/1220836‑1221107

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