Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,221,641 T→A 7.4% pseudogene (2041/2574 nt) ycgH → pseudogene; putative ATP‑binding component of a transport system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,221,6410TA 7.4% 165.2 / 4.2 54pseudogene (2041/2574 nt)ycgHpseudogene; putative ATP‑binding component of a transport system
Reads supporting (aligned to +/- strand):  ref base T (29/21);  new base A (2/2);  total (31/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.53e-01

GGCGATCTGTTAACCGGTGGTTTTACCGATAGTGACAGTTGGCGTTTGGGAGTGATGGCTGGTTATGCCCGCGACTACAACTTAACTCATTCCAGCGTGTCGGATTATCGTTCGAAAGGGAGTGTCAGAGGCTATAGCGCAGGGCTGTATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGGCGCATACATTGACTCCTGGGCGCAATATAGCTGGTTTAAAAACTCGGTGAAAGGGGATGAATTAGCCTATGAATCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAG  >  NC_000913/1221494‑1221788
                                                                                                                                                   |                                                                                                                                                   
gtcGATATGTTTTCGGGTGGTTTTACGGATAGTGACAGTTGTCGTTTGGGAGTGATGGCTGGTTTTGCCAGCGACTACAACTTAACTCATTCGAGCGTGTCGGATTATCGTTCGATAGTGATTGTCAGAGGCTATAGCGCAGGGCTGt                                                                                                                                                     <  7:443392/146‑1 (MQ=255)
ggCGATCTGTTACCCGGTGGGTTTACCGATAGTGACAGTTGGCGTTTGGGAGTGATGGCTGGTTATGCCCGCGACTACAACTTAACTCATTCCAGCGTGTCGGATTATCGTTCGAAAGGGAGTGTCAGAGGCTATAGCGCAGGGCTGt                                                                                                                                                     <  7:358961/148‑1 (MQ=255)
ggCGATCTGTTAACCGGTGGTTTTACCGATAGTGACAGTTGGCGTTTGGGAGTGATGGCTGGTTATGCCCGCGACTACAACTTAACTCATTCCAGCGTGTCGGATTATCGTTCGAAAGGGAGTGTCAGAGGCTATAGCGCAGGGCTGt                                                                                                                                                     <  2:77291/148‑1 (MQ=255)
ggCGATCTGTTAACCGGTGGTTTTACCGATAGTGACAGTTGGCGTTTGGGAGTGATGGCTGGTTATGCCCGCGACTACAACTTAACTCATTCCAGCGTGTCGGATTATCGTTCGAAAGGGAGTGTCAGAGGCTATAGCGCAGGGCTGt                                                                                                                                                     <  7:177457/148‑1 (MQ=255)
  cGATCTGTTAACCGGTGGTTTTACCGATAGTGACAGTTGGCGTTGGGGAGTGATGGCTGGTTATGCCCGCGACAACAACTTAACTCATTCCAGCGTGTCGGATTATCGTTCGAAAGGGAGTGTCAGAGGCTATAGCGCAGGGCTGTAt                                                                                                                                                   <  2:280966/148‑1 (MQ=255)
         ttAACCGGTGGTTTTACCGATAGTGACAGTTGGCGTTGGGGAGTGATGGCTGGTTATGCCCGCGACTACAACTTAACTCATTCCAGCGTGTCGGATTATCGTTCGAAAGGGAGTGTCAGAGGCTATAGCGCAGGGCGGTATGcc                                                                                                                                                <  5:294150/144‑1 (MQ=255)
          tAACCGGTGGTTTTACCGATAGTGCCAGTTGGCGTTTGGGAGTGATGGCTGGTTATGCACGCGACTACAACTTAACTCATTCCAGCGTGTCGGATTATCGTTCGAAAGGGAGTGTCAGAGGCTATAGCGCAGGGCTGTATGCCACCTg                                                                                                                                           >  8:400539/1‑148 (MQ=255)
               ggtggtTGTACCGATAGTGACAGTTGGCGTTTGGGAGTGATGGCTGGTTATGCCCGCGACTACAACTTAACTCATTCCAGCGTGTCGGATTATCGTTCGAAAGGGAGTGTCAGAGGCTATAGCGCAGGGCTGt                                                                                                                                                     <  5:391823/133‑1 (MQ=255)
                gtggtTTTACCGATAGTGACAGTTGGCGTTTGGGAGTGATGGCTGGTTATGCCCGCGACTACAACTTAACTCATTCCAGCGTGTCGGAGTATCGTTCGAAAGGGAGTGTCAGAGGCTAAAGCGCAGGGGTGTa                                                                                                                                                    >  6:391823/1‑133 (MQ=255)
                   gtTTTACCGATAGTGACAGTTGGCGTTTGGGAGTGATGGCTGGTTATGCCCGCGACTACAACTTAACTCATTCCAGCGTGTCGGATTATCGTTCGAAAGGGAGTGTCAGAGGCTATAGCGCAGGGCTGTATGCCACCTGGTTTGCCGa                                                                                                                                  <  3:198757/148‑1 (MQ=255)
                                 gACAGTTGGCGTTTGGGAGTGATGGCTGGTTATGCCCGCGACTACAACTTAACTCATTCCAGCGTGTCGGATTATCGTTCGAAAGGGACTGTCAGAGGCTATAGCGCAGGGCTGTATGCCACCTGGTTTGCCGATGACATCAGTaaa                                                                                                                     >  3:88405/1‑147 (MQ=255)
                                      ttGGCGTTTGGGAGTGATGGCTGGTTATGCCCGCGACTACAACTTAACTCATTCCAGCGTGTCGGATTATCGTTCGAAAGGGAGTGTCAGAGGCTATAGCGCAGGGCTGt                                                                                                                                                     >  2:121743/1‑110 (MQ=255)
                                      ttGGCGTTTGGGAGTGATGGCTGGTTATGCCCGCGACTACAACTTAACTCATTCCAGCGTGTCGGATTATCGTTCGAAAGGGAGTGTCAGAGGCTATAGCGCAGGGCTGt                                                                                                                                                     <  1:121743/110‑1 (MQ=255)
                                           gTTTGGGAGTGATGGCGGGTTATGCCCGCGACTACAACTTAACTCATTCCAGCGTGTCGGATTATCGTTCGAAAGGGAGTGTCAGAGGCTATAGCGCAGGGCTGTATGCCACCTGGTTTACCGATGACATCAGTAAAAAAGGCGCATa                                                                                                          <  5:14206/148‑1 (MQ=255)
                                             ttGGGAGTGATGGCTGGTTATGCCCGCGACTACAACTTAACTCATTCCAGCGTGTCGGATTATCGTTCGAAAGGGAGTGTCAGAGGCTATAGCGCAGGGCTGt                                                                                                                                                     >  4:155658/1‑103 (MQ=255)
                                             ttGGGAGTGATGGCTGGTTATGCCCGCGACTACAACTTAACTCATTCCAGCGTGTCGGATTATCGTTCGAAAGGGAGTGTCAGAGGCTATAGCGCAGGGCTGt                                                                                                                                                     <  3:155658/103‑1 (MQ=255)
                                                ggAGTGATGGCTGGTTATGCCCGCGACTACAACTTAACTCATTCCAGCGTGTCGGATTATCTTTCGAAAGGGAGTGTCAGAGGCTATAGCGCAGGGCTGTATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGGCGCATACATTg                                                                                                     <  1:244379/148‑1 (MQ=255)
                                                  aGTGATGGCTGGTTATGCCCGCGACTACAACTTAACTCATTCCAGCGTGTCGGATTATCGTTCGAAAGGGAGTGTCAGAGGCTATAGCGCAGGGCTGTATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGGCGCATACATTGAc                                                                                                   >  7:344481/1‑148 (MQ=255)
                                                     gATGGCTGGTTATGCCCGCGACTACAACTTAACTCATTCCAGCGTGTCGGATTATCGTTCGAAAGGGAGTGTCAGAGGCTATAGCGCAGGGCTGTATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGGCGCATACATTg                                                                                                     <  1:223539/143‑1 (MQ=255)
                                                     gATGGCTGGTTATGCCCGCGACTACAACTTAACTCATTCCAGCGTGTCGGATTATCGTTCGAAAGGGAGTGTCAGAGGCTATAGCGCAGGGCTGTATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGGCGCATACATTg                                                                                                     >  2:223539/1‑143 (MQ=255)
                                                     gATGGCTGGTTATGCCCGCGACTACAACTTAACTCATTCCAGCGTGTCGGATTATCGTTCGAAAGGGAGTGTCAGAGGCTATAGCGCAGGGCTGTATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGGCGCATACATTGACTcc                                                                                                >  5:315243/1‑148 (MQ=255)
                                                       tGGCTGGTTATGCCCGCGACTACAACTTAACTCATTCCAGCGTGTCGGATTATCGTTCGAAAGGGAGTGTCAGAGGCTATAGCGCAGGGCTGTATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGGCGCATACATTGACTccgg                                                                                              >  5:51720/1‑146 (MQ=255)
                                                       tGGCTGGTTATGCCCGCGACTACAACTTAACTCATTCCAGCGTGTCGGATTATCGTTCGAAAGGGAGTGTCAGAGGCTATAGCGCAGGGCTGTATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGGCGCATACATTGACTCCt                                                                                               <  3:48388/147‑1 (MQ=255)
                                                            ggTTATGCCCGCGACTACAACTTAACTCATTCCAGCGTGTCGGATTATCGTTCGAAAGGGAGTGTCAGAGGCTATAGCGCAGGGCTGTATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAgg                                                                                                                <  1:214275/125‑1 (MQ=255)
                                                            ggTTATGCCCGCGACTACAACTTAACTCATTCCAGCGTGTCGGATTATCGTTCGAAAGGGAGTGTCAGAGGCTATAGCGCAGGGCTGTATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAgg                                                                                                                >  2:214275/1‑125 (MQ=255)
                                                                 tGCCCGCGACTACAACTTAACTCATTCCAGCGTGTCGGATTATCGTTCGAAAGGGAGTGTCAGAGGCTATAGCGCAGGGCTGTATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGGCGCATACATTGACTCCTGGGCGCAata                                                                                     <  4:88405/147‑1 (MQ=255)
                                                                   cccGCGACTACAACTTAACTCATTCCAGCGTGTCGGATTATCGTTCGAAAGGGAGTGTCAGAGGCTATAGCGCAGGGCTGTATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGGCGCATACATTGACTCCTGGGCGCAATATAg                                                                                  >  2:3980/1‑148 (MQ=255)
                                                                           tACAACTTAACTCATTCCAGCGTGTCGGATTATCGTTCGAAAGGGAGTGTCAGAGGCTATAGCGCAGGGCTGTATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGGCGCATACATTGACTCCTGGGCGCAATATAGCTGGttt                                                                           >  7:279717/1‑147 (MQ=255)
                                                                            aCAACTTAACTCATTCCAGCGTGTCGGATTATCGTTCGAAAGGGAGTGTCAGAGGCTATAGCGCAGGGCTGt                                                                                                                                                     <  1:204158/72‑1 (MQ=255)
                                                                            aCAACTTAACTCATTCCAGCGTGTCGGATTATCGTTCGAAAGGGAGTGTCAGAGGCTATAGCGCAGGGCTGt                                                                                                                                                     >  2:204158/1‑72 (MQ=255)
                                                                            aCAACTTAACTCATTCCAGCGTGTCGGATTATCGTTCGAAAGGGAGTGTCAGAGGCTATAGCGCAGGGCTGTATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGGCGCATACATTGACTCCTGGGCGCAATATAGCTGGTTTa                                                                          >  5:41276/1‑147 (MQ=255)
                                                                               aCTTAACTCATTCCAGCGTGTCGGATTATCGTTCGAAAGGGAGTGTCAGAGGCTATAGCGCAGGGCTGTATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGGCGCATACATTGACTCCTGGGCGCAATATAGCTGGTTTaaaaa                                                                      >  1:116770/1‑148 (MQ=255)
                                                                               aCTTAACTCATTCCAGCGTGTCGGATTATCGTTCGAAAGGGAGTGTCAGAGGCTATAGCGCAGGGCTGTATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGGCGCATACATTGACTCCTGGGCGCAATATAGCTGGTTTaaaa                                                                       >  2:153498/1‑147 (MQ=255)
                                                                               aCTTAACTCATTCCAGCGTGTCGGATTATCGTTCGAAAGGGAGTGTCAGAGGCTATAGCGCAGGGCTGTATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGGCGCATACATTGACTCCTGGGCGCAATATAGCTGGTTTaaaa                                                                       <  1:153498/147‑1 (MQ=255)
                                                                                   aaCTCATTCCAGCGTGTCGGATTATCGTTCGAAAGGGAGTGTCAGAGGCTATAGCGCAGGGCTGTATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGGCGCATACATTGACTCCTGGGCGCAATATAGCTGGTTTAAAAACTc                                                                   <  3:248124/147‑1 (MQ=255)
                                                                                    aCTCATTCCAGCTTGGCGGATTATCGTTCGAAAGGGAGTGTCAGAGGCTATAGCGCAGGTCTGTATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGGCGCATACATTGACTCCTGGGCGCAATATAGCTGGTTTAAAAACTCgg                                                                 <  2:320527/148‑1 (MQ=255)
                                                                                    aCTCATTCCAGCGTGTCGGATTATCGTTCGAAAGGGAGTGTCAGAGGCTATAGCGCAGGGCTGTATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGGCGCATACATTGACTCCTGGGCGCAATATAGCTGGTTTAAAAACTCg                                                                  <  6:315243/147‑1 (MQ=255)
                                                                                      tCATTCCAGCGTGTCGGATTATCGTTCGAAAGGGAGTGTCAGAGGCTATAGCGCAGGGCTGTATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGGCGCATACATTGACTCCTGGGCGCAATATAGCTGGTTTAAAAACTCGGt                                                                >  1:297005/1‑147 (MQ=255)
                                                                                       cATTCCAGCGTGTCGGATTATCGTTCGAAAGGGAGTGTCAGAGGCTATAGCGCAGGGCTGTATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGGCGCATACATTGACTCCTGGGCGCAATATAGCTGGTTTAAAAACTCGGTg                                                               >  8:57796/1‑147 (MQ=255)
                                                                                       cATTCCAGCGTGTCGGATTATCGTTCGAAAGGGAGTGTCAGAGGCTATAGCGCAGGGCTGTATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGGCGCATACATTGACTCCTGGGCGCAATATAGCTGGTTTAAAAACTCGGTg                                                               >  4:55266/1‑147 (MQ=255)
                                                                                        aTTCCAGCTGGTCGGATTATCGTTCGAAAGGGAGTGTCAGAGGCTATAGCGCAGGGCTGTATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGGCGCATACATTGACTCCTTGGCGCAATATAGCTGGTTTAAAAACTCGGTGaa                                                             <  1:52066/148‑1 (MQ=255)
                                                                                         ttCCAGCGTGTCGGATTATCGTTCGAAAGGGAGTGTCAGAGGCTATAGCGCAGGGCTGAATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGGCGCAt                                                                                                           >  8:126571/1‑101 (MQ=255)
                                                                                         ttCCAGCGTGTCGGATTATCGTTCGAAAGGGAGTGTCAGAGGCTATAGCGCAGGGCTGAATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGGCGCAt                                                                                                           >  8:132878/1‑101 (MQ=255)
                                                                                         ttCCAGCGTGTCGGATTATCGTTCGAAAGGGAGTGTCAGAGGCTATAGCGCAGGGCTGAATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGGCGCAt                                                                                                           <  7:132878/101‑1 (MQ=255)
                                                                                         ttCCAGCGTGTCGGATTATCGTTCGAAAGGGAGTGTCAGAGGCTATAGCGCAGGGCTGAATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGGCGCAt                                                                                                           <  7:126571/101‑1 (MQ=255)
                                                                                          tCCAGCGTGTCGGATTATCGTTCGAAAGGGAGTGTCAGAGGCTATAGCGCAGGGCTGTATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGGCGCATACATTGACTCCTGGGCGCAATATAGCTGGTTTAAAAACTCGGTGAAAg                                                           >  7:368001/1‑148 (MQ=255)
                                                                                          tCCAGCGTGTCGGATTATCGTTCGAAAGGGAGTGTCAGAGGCTATAGCGCAGCTCTGTATTCCCCCTGGTTTGCCGATGACATCAGTAAAAAAGGAGCATACATTGACTCCTTGGCGCAATAAAGAGGGTTTAAAAACACCGTGAAAg                                                           >  5:140013/1‑148 (MQ=255)
                                                                                           ccAGCGTGTCGGATTATCGTTCGAAAGGGAGTGTCAGAGGCTATAGCGCAGGGCTGTATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGGCGCATACATTGACTCCTGGGCGCAATATAGCTGGTTTAAAAACTCGGTGAAAg                                                           <  2:297005/147‑1 (MQ=255)
                                                                                                 tgtCGGATTATCGTTCGAAAGGGAGTGTCAGAGGCTATAGCGCAGGGCTGTATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGGCGCATACATTGACTCCTGGGCGCAATATAGCTGGTTTAAAAACTCGGTGAAAGGGGATGa                                                    <  7:400539/148‑1 (MQ=255)
                                                                                                                cGAAAGGGAGTGTCAGAGGCTATAGCGCAGGGCTGt                                                                                                                                                     <  1:115899/36‑1 (MQ=255)
                                                                                                                cGAAAGGGAGTGTCAGAGGCTATAGCGCAGGGCTGt                                                                                                                                                     >  2:115899/1‑36 (MQ=255)
                                                                                                                  aaaGGGAGTGTCAGAGGCTATAGCGCAGGGCTGTATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGGCGCATACATTGACTCCTGGGCGCAATATAGCTGGTTTAAAAACTCGGTGAAAGGGGATGAATTAGCCTATGAATCCt                                   <  2:116770/148‑1 (MQ=255)
                                                                                                                    aGGGAGTGTCAGAGGCTATAGCGCAGGGCTGTATGCCACCTGGTTGCCGATGACATCAGTAAAAAAGgcgc                                                                                                             >  3:371348/1‑71 (MQ=255)
                                                                                                                    aGGGAGTGTCAGAGGCTATAGCGCAGGGCTGTATGCCACCTGGTTGCCGATGACATCAGTAAAAAAGgcgc                                                                                                             <  4:371348/71‑1 (MQ=255)
                                                                                                                           gtCAGAGGCTATAGCGCAGGGCTGTATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGgcgc                                                                                                             >  1:175072/1‑65 (MQ=255)
                                                                                                                           gtCAGAGGCTATAGCGCAGGGCTGTATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGgcgc                                                                                                             <  2:175072/65‑1 (MQ=255)
                                                                                                                           gtCAGACGCTATAGCGCAGGGCTGTATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGGCGCATACATTGACTCCTGGGCGCAATATAGCAGGTTTAAAAACTCGGTGAAAGTGGATGAATTAGACTATGAATCCTAAAgcgag                           >  2:89229/1‑145 (MQ=255)
                                                                                                                                     ataGCGCAGGGCTGTATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGGCGCATACATTGACTCCTGGGCGCAAGATAGCTGGTTTAAAAACTCGGTGAAAGGGGATGAATTAGCCTATGAATCCTATAGCGCGAAAGGTGCAAc                >  6:316537/1‑148 (MQ=255)
                                                                                                                                      taGCGCAGGGCTGTATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGGCGCATACATTGACTCCTGGGCGCAAtata                                                                                   >  1:322633/1‑80 (MQ=255)
                                                                                                                                      taGCGCAGGGCTGTATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGGCGCATACATTGACTCCTGGGCGCAAtata                                                                                   <  2:322633/80‑1 (MQ=255)
                                                                                                                                       aGCGCAGGGCTGTATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGGCGCATACATTGACTCCTGGGCGCAATATAGCTGGTTTAAAAACTCGGTGAAAGGGGATGAATTAGCCTATGAATCCTATAGCGCGAAAGGTGCAACCg              >  7:7220/1‑148 (MQ=255)
                                                                                                                                       aGCGCAGGGCTGTATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGGCGCATACATTGACTCCTGGGCGCAATATAGCTGGTTTAAAAACTCGGTGAAAGGGGATGAATTAGCCTATGAATCCTATAGCGCGAAAGGTGCAACCg              >  4:359517/1‑148 (MQ=255)
                                                                                                                                          gcAGGGCTGTATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGGCGCATACATTGACTCCTGGGCGCAATATAGCTGGTTTAAAAACTCGGTGAAAGGGGATGAATTAGCCTATGAATCCTATAGCGCGAAAGGTGCAACCGtct           >  5:14729/1‑148 (MQ=255)
                                                                                                                                            aGGTCTGTATGACACCTGGTTTTCCGATGACATCAGTAAAAAATGCTCATTCATTGACTACTGGGCGCAATATAGCTGGTTTATATACTCGGTGAATGGGGATTAATTATCCTATGAATCCTATAGCGCGAAAGGTGCAACGGTCTCg         <  6:140013/148‑1 (MQ=255)
                                                                                                                                               gCTGTATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGGCGCATACATTGACTCCTGGGCGCAATATAGCTGGTTTAAAAACTCGGTGAAAGGGGATGAATTATCCTATGAATCCTATAGCGCGAAAGGTGCAACCGTCTCGCTg      >  6:163120/1‑148 (MQ=255)
                                                                                                                                                  gTATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGGCGCATACATTGACTCCTGGGCGCAATATAGCTGGTTTAAAAACTCGGTGAAAGGGGATGAATTAGCCTATGAATCCTATAGCGCGAAAGGTGCAACCGTCTCGCTgg     >  6:379751/1‑146 (MQ=255)
                                                                                                                                                  gTATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGGCCCATACATTGACTCCTGGGCGCAATATAGCTGGTTTAATAACTCGGTGAAAGGGGATGAATTAGCCTATGAATCCTATAGCGCGAAAGGTGCAACCGTCTCTCTgg     >  3:57436/1‑146 (MQ=255)
                                                                                                                                                   tATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGGCGCATACATTGACTCCTGGGCGCAATATAGCTGGGTTAAAAACTCGGTGAAAGGGGATGAATTAGCCTATGAATCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAg  <  7:403968/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GGCGATCTGTTAACCGGTGGTTTTACCGATAGTGACAGTTGGCGTTTGGGAGTGATGGCTGGTTATGCCCGCGACTACAACTTAACTCATTCCAGCGTGTCGGATTATCGTTCGAAAGGGAGTGTCAGAGGCTATAGCGCAGGGCTGTATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGGCGCATACATTGACTCCTGGGCGCAATATAGCTGGTTTAAAAACTCGGTGAAAGGGGATGAATTAGCCTATGAATCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAG  >  NC_000913/1221494‑1221788

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