Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 945,863 T→C 5.8% intergenic (‑306/‑8) ycaC ← / → ycaD putative isochorismatase family hydrolase/putative MFS transporter, inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913945,8630TC 5.8% 220.0 / 3.5 71intergenic (‑306/‑8)ycaC/ycaDputative isochorismatase family hydrolase/putative MFS transporter, inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (29/37);  new base C (2/2);  total (32/39)
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.77e-01

GATGAGACTAATCATCTCTTCAGCATTAGCGTTATGAATGAGTAATATAAATCGATAATACCGCCGCTGGGATATCGCGTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCCACATGGCAGGTAGGCGTTGTC  >  NC_000913/945716‑946008
                                                                                                                                                   |                                                                                                                                                 
gATGAGACTAATCATCTCTTCAGCATTAGCGTTATGAATGAGTAATATAAATCGATAATACCGCCGCTGGGATATCGCGTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGtt                                                                                                                                                   <  2:410696/148‑1 (MQ=255)
 aTGAGACTAATCATCTCTTCAGCATTAGCGTTATGAATGAGTAATATAAATCGATAATACCGCCGCTGGGATATCGCGTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGttt                                                                                                                                                  >  2:210061/1‑148 (MQ=255)
    agaCTAATCATCTCTTCAGCATTAGCGTTATGAATGAGTAATATAAATCGATAATACCGCCGCGGGGATATCGCGTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTa                                                                                                                                                <  2:295969/147‑1 (MQ=255)
      aCTAATCATCTCTTCAGCATTAGCGTTATGAATGAGTAATATAAATCGATAATACCGCCGCTGGGATATCGCGTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGGATGAAATTTGAACTTCCTCACGGTTTTAAtt                                                                                                                                             <  7:189954/148‑1 (MQ=255)
      aCTAATCATCTCTTCAGCATTAGCGTTATGAATGAGTAATATAAATCGATAATACCGCCGCTGGGATATCGCGTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAAtt                                                                                                                                             <  7:174558/148‑1 (MQ=255)
       atAATCATCTCTTCAGCATTAGCGTTATGAATTAGTAATATAAATCGATAATACCGCCGCTGGGATATCGCGTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTc                                                                                                                                            <  5:212354/147‑1 (MQ=255)
          atcatcTCTTCAGCATTAGCGTTATGAATGAGTAATATAAATCGATAATACCGCCGCTGGGATATCGCGTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCa                                                                                                                                           <  4:291648/146‑1 (MQ=255)
                   tCAGCATTAGCGTTATGAATGAGTAATATAAATCGATAATACCGCCGCTGGGATATCGCGTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGTCCACGta                                                                                                                                 >  3:36807/1‑147 (MQ=255)
                    cAGCATTAGCGTTATGAATGAGTAATATAAATCGATAATACCGCCGCTGGGATATCGCGTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGTCCACGt                                                                                                                                  <  8:172017/145‑1 (MQ=255)
                         ttAGCGTTATGAATGAGTAATATAAATCGATAATACCGCCGCTGGGATATCGCGTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGTCCACGTATACCCAg                                                                                                                          <  2:431686/148‑1 (MQ=255)
                          tAGCGTTATGAATGAGTAATATAAATCGATAATACCGCCGCTGGGATATCGCGTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGTCCACGTATACCCAGc                                                                                                                         >  7:346663/1‑148 (MQ=255)
                          tAGCGTTATGAATGAGTAATATAAATCGATAATACCGCCGCTGGGATATCGCGTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGTCCACGTATACCCAGc                                                                                                                         >  4:297797/1‑148 (MQ=255)
                          tAGCGTTATGAATGAGTAATATAAATCGATAATACAGCCGCTGGGATATCGCGTATTTTCACCCATTGACAATGTTATTGGCGGTTGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGTCCACGTATACCCAGc                                                                                                                         >  3:343645/1‑148 (MQ=255)
                               ttATGAATGAGTAATATAAATCGATAATACCGCCTCTTGGATATCGCGTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTc                                                                                                                    <  4:343645/148‑1 (MQ=255)
                                      tgaGTAATATAAATCGATAATACCGCCGCTGGGATATCGCGTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGATATTTGAACTTCCTCACGGTTTTAATTCATGTCCTCGTTTACCCAGCCTGTCATGCttt                                                                                                             <  8:273339/148‑1 (MQ=255)
                                      tgaGTAATATAAATCGATAATACCGCCGCTGGGATATCGCGTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGGTGTAATTCATGTCCACGTATACCCAGCCTGTCATGCttt                                                                                                             <  1:317746/148‑1 (MQ=255)
                                       gaGTAATATAAATCGATAATACCGCCGCTGGGATATCGCGTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTc                                                                                                                                            <  2:400434/116‑1 (MQ=255)
                                       gaGTAATATAAATCGATAATACCGCCGCTGGGATATCGCGTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTc                                                                                                                                            >  1:400434/1‑116 (MQ=255)
                                        aGTAATATAAATCGATAATACCGCCGCTGGGATATCGCGTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCtttt                                                                                                            <  7:425472/147‑1 (MQ=255)
                                             tataAATCGATAATACCGCCGCTGGGATATCGCGTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTc                                                                                                      <  4:232072/148‑1 (MQ=255)
                                              ataAATCGATAATACCGCCGCTGGGATATCGCGTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCt                                                                                                     >  7:433832/1‑148 (MQ=255)
                                                    cGATAATACCGCCGCTGGGATATCGCGTATTTTCACCCATTGACAAGGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTg                                                                                               >  4:358242/1‑148 (MQ=255)
                                                      ataataCCGCCGCTGGGATATCGCGTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTATGGACTGCt                                                                                             >  5:349196/1‑148 (MQ=255)
                                                         ataCCGCCGCTGGGATATCGCGTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCtttt                                                                                          <  3:297797/148‑1 (MQ=255)
                                                            ccgccgCTGGGATATCGCGTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTttgtt                                                                                       >  5:327875/1‑148 (MQ=255)
                                                            ccgccgCTGGGATATCGCGTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTttgt                                                                                        <  1:384157/147‑1 (MQ=255)
                                                            ccgccgCTGGGATATCGCGTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTttgt                                                                                        <  4:352374/147‑1 (MQ=255)
                                                             cgccgaGGGGATTTCGGTTATTTTCACCCATTTACAATTTTTTTGGCGGTGGCTTGATGAGCATTAAATTTGACCTTCCTACCTGTTTTAATTAATGTCCAATTATACCACGCCTGTCATGATTTTGCTGTCTGGCCTGCTTttgttg                                                                                      <  8:385937/148‑1 (MQ=255)
                                                                     ggATATCGCGTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTg                                                                                <  8:346663/146‑1 (MQ=255)
                                                                      gATATCGCGTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGc                                                                              >  1:227507/1‑147 (MQ=255)
                                                                       atatCGCGTATTTTCACACATTGACAATGTTTTTGGCGGTGTCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGa                                                                            >  6:300058/1‑148 (MQ=255)
                                                                           cgcgTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGc                                                                        >  5:100096/1‑148 (MQ=255)
                                                                           cgcgTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGAATCTGGCGATTGc                                                                        >  1:331539/1‑148 (MQ=255)
                                                                            gcgTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCg                                                                       >  1:64048/1‑148 (MQ=255)
                                                                             cgTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCgg                                                                      <  8:433832/148‑1 (MQ=255)
                                                                                 ttttCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGgtgt                                                                   <  6:327875/147‑1 (MQ=255)
                                                                                     cACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTaaa                                                               <  3:19460/147‑1 (MQ=255)
                                                                                      aCCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGGTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATa                                                             >  4:88434/1‑148 (MQ=255)
                                                                                       cccATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTaa                                                                                                                                               <  4:26641/65‑1 (MQ=255)
                                                                                       cccATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTaa                                                                                                                                               >  3:26641/1‑65 (MQ=255)
                                                                                           ttGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATACACTc                                                        <  6:100096/148‑1 (MQ=255)
                                                                                            tGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGACCTTCCTCACGGTTTTAATTCATGTCCCCGTATACCCCGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGACGCTTTCGGTGTTTAATACACACg                                                       >  6:100175/1‑148 (MQ=255)
                                                                                                aaTGTTTATGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTATTAATTCATGTCCACGTATACACAGCCTCACATCCTTATGCAGTCTCGCCTGCTTATCATGAATCAGGCGAATGCCCTGTTAAATACACTCCTGcc                                                   >  3:186701/1‑148 (MQ=255)
                                                                                                 ttGTTTTTGGCGGTGGCATGATGCTCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTTTGGCGATTGCGGTGTTAAATACACTCGTGCCg                                                  <  7:131133/147‑1 (MQ=255)
                                                                                                 aTGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTCTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATaca                                                           >  7:276211/1‑139 (MQ=255)
                                                                                                 aTGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTCTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATaca                                                           >  7:269489/1‑139 (MQ=255)
                                                                                                 aTGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTCTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATaca                                                           <  8:269489/139‑1 (MQ=255)
                                                                                                 aTGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTCTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATaca                                                           <  8:276211/139‑1 (MQ=255)
                                                                                                       ttGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATACACTCGTGCCGCTTTg                                             <  2:227507/147‑1 (MQ=255)
                                                                                                        tGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATACACTCGTGCCGCTTTgg                                            <  3:358242/147‑1 (MQ=255)
                                                                                                             gTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATacac                                                          <  2:307468/128‑1 (MQ=255)
                                                                                                             gTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATacac                                                          >  1:307468/1‑128 (MQ=255)
                                                                                                               ggCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGcc                                      >  2:435326/1‑146 (MQ=255)
                                                                                                                    gatgCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGaa                                >  7:91471/1‑147 (MQ=255)
                                                                                                                     atgCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATACACTCGTGCCGCTTTGGc                                           >  7:6531/1‑135 (MQ=255)
                                                                                                                     atgCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATACACTCGTGCCGCTTTGGc                                           <  8:6531/135‑1 (MQ=255)
                                                                                                                     atgCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAac                               <  1:208027/147‑1 (MQ=255)
                                                                                                                     atgCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGAGTTTGTTGACTCTGGCGATTGCGGTGGTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAaca                              >  3:347928/1‑148 (MQ=255)
                                                                                                                      tgCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAacac                             <  8:91471/148‑1 (MQ=255)
                                                                                                                        cgcATGAAATTTGAACTTCCTCACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTaa                                                                <  4:246306/111‑1 (MQ=255)
                                                                                                                        cgcATGAAATTTGAACTTCCTCACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTaa                                                                >  3:246306/1‑111 (MQ=255)
                                                                                                                                      aCTTCCTCACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGACTGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGTAACACATGTCCACATGGCAg              <  5:300058/147‑1 (MQ=255)
                                                                                                                                       cTTCCTCACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGa                                                                                     >  7:313091/1‑75 (MQ=255)
                                                                                                                                       cTTCCTCACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGa                                                                                     <  8:313091/75‑1 (MQ=255)
                                                                                                                                         tCCTCACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCcac                     >  4:180406/1‑137 (MQ=255)
                                                                                                                                         tCCTCACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCcac                     <  3:180406/137‑1 (MQ=255)
                                                                                                                                         tCCTCACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCCACATgg                 >  6:268965/1‑141 (MQ=255)
                                                                                                                                         tCCTCACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCCACATgg                 <  5:268965/141‑1 (MQ=255)
                                                                                                                                            tCACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCCACATGGCAGGTAcgcg       >  4:234537/1‑148 (MQ=255)
                                                                                                                                             cACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCCACATGGCAGGTAGGCGt      >  5:360708/1‑148 (MQ=255)
                                                                                                                                                 gTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGGTTGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGccc                                     >  2:205547/1‑113 (MQ=255)
                                                                                                                                                 gTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGGTTGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGccc                                     <  1:205547/113‑1 (MQ=255)
                                                                                                                                                 gTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCCACATGGCAGGTAGGCGTTGTc  <  4:36807/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GATGAGACTAATCATCTCTTCAGCATTAGCGTTATGAATGAGTAATATAAATCGATAATACCGCCGCTGGGATATCGCGTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCCACATGGCAGGTAGGCGTTGTC  >  NC_000913/945716‑946008

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