Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,878,344 C→T 5.0% intergenic (+726/+225) iap → / ← ygbF aminopeptidase in alkaline phosphatase isozyme conversion/CRISPR adaptation ssRNA endonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,878,3440CT 5.0% 253.4 / 2.3 80intergenic (+726/+225)iap/ygbFaminopeptidase in alkaline phosphatase isozyme conversion/CRISPR adaptation ssRNA endonuclease
Reads supporting (aligned to +/- strand):  ref base C (40/35);  new base T (2/2);  total (43/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.83e-01

TCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACATTATCAATTACAACCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGCGTGAGCGTATCGCCGCGCGTCTGCGAAAGCGGTTTATCCCCGCTGGCGCGGGGAACTCTCTAAAAGTATACATTTGTTCTTA  >  NC_000913/2878197‑2878487
                                                                                                                                                   |                                                                                                                                               
tCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACATTATCAATTACAACCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGc                                                                                                                                                 <  6:242544/148‑1 (MQ=255)
tCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACATTATCAATTACAACCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGc                                                                                                                                                 <  4:202841/148‑1 (MQ=255)
  cccGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACATTATCAATTACAACCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCgtg                                                                                                                                              >  8:297700/1‑149 (MQ=255)
  cccGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACATTATCAATTACAACCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCgtg                                                                                                                                              >  8:186571/1‑149 (MQ=255)
   ccGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTACAGACACGATACAAACGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACATTATCAATTAAAACAGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACGCAGCgtgt                                                                                                                                             >  8:167076/1‑149 (MQ=255)
     gcgtgcgaGGGTAACTCGAGACGGCGCTGAAATTAGAGAAGCGATCCAATCGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACTTTATCAATTACATCCGTCAGGGAGCGCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCGTGTTc                                                                                                                                           <  4:201788/145‑1 (MQ=255)
              gggAACTCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACATTATCAATTACAACCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCAcc                                                                                                                                  >  2:275681/1‑149 (MQ=255)
                  acgcgcgTCCGCTCAGATATTCCAGACCCGATCCACCCGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACATTCTCAATTACACCCGTCAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCTTTg                                                                                                                              <  8:364081/146‑1 (MQ=255)
                    tCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACATTATCAATTACAACCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCttt                                                                                                                               >  1:83863/1‑146 (MQ=255)
                    tCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACATTATCAATTACAACCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCttt                                                                                                                               <  2:83863/146‑1 (MQ=255)
                      gcgACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACATTATCAATTACAACCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCTTTGGCtt                                                                                                                          >  7:364630/1‑149 (MQ=255)
                       cgACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACATTATCAATTACAACCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCTTTGGCTTc                                                                                                                         <  4:258253/149‑1 (MQ=255)
                              tCAGAAATTCCAGACCCGATCCAAACGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACATTATCAATTACAACCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCg                                                                                                                  >  8:73537/1‑149 (MQ=255)
                              tCAGAAATTCCAGACCCGATCCAAACGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACATTATCAATTACAACCGACAGGGAGCACGGTTTATCCCCGCTGGCGCGGGGAACTAAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCg                                                                                                                  >  8:219108/1‑149 (MQ=255)
                              tCAGAAATTCCAGACCCGATCCAAACGGATTATCCCCGCTGGCGCGGGGAACTCTCAACATTATCAATTACAACCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCg                                                                                                                  >  7:305533/1‑149 (MQ=255)
                                    aTTCCAGACCCGATCCAAACGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACATTATCAATTACAACCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTAt                                                                                                            >  1:217040/1‑149 (MQ=255)
                                      tCCAGACCCGATCCAAACGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACATTATCAATTACAACCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATcc                                                                                                          >  7:100532/1‑149 (MQ=255)
                                        cAGACCCGATCCAAACGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACATTATCAATTACAACCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATcccc                                                                                                        >  1:347090/1‑149 (MQ=255)
                                           aCCCGATCCAAACGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACATTATCAATTACAACCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCt                                                                                                     <  8:21659/149‑1 (MQ=255)
                                             ccGATCCAAACGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACATTATCAATTACAACCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCACCGCTgg                                                                                                   >  1:183610/1‑149 (MQ=255)
                                               gATCCAAACGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACATTATCAATTACAACCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGgcg                                                                                                 >  1:353027/1‑149 (MQ=255)
                                                   cAAACGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACATTATCAATTACAACCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCggg                                                                                             <  7:7170/149‑1 (MQ=255)
                                                    aaaCGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACATTATCAATTACAACCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCggg                                                                                             >  7:322856/1‑148 (MQ=255)
                                                         gTTTATCCCCGCTGGCGCGGGGAACTCTCAACATTATCAATTACAACCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAActc                                                                                       >  4:242729/1‑149 (MQ=255)
                                                             aTCCCCGCTGGCGCGGGGAACTCTCAACATTATCAATTACAACCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGc                                                                                    <  3:242729/148‑1 (MQ=255)
                                                              tCCCCGCTGGCGCGGGGAACTCTCCACATTATCAATTACAACCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGCGt                                                                                  <  7:7115/149‑1 (MQ=255)
                                                                         gcgGGGAACTCTCAACATTATCAATTACAACCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGCGTGAGCGTATcgc                                                                       <  7:186571/149‑1 (MQ=255)
                                                                           ggggAACTCTCAACATTATCAATTACAACCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGCGTGAGCGTATcgccg                                                                     >  6:90552/1‑149 (MQ=255)
                                                                             ggAACTCTGAACATTATAAATTACAACCGACAGGGAGCCCGGTTTTTCCCCGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAAGTCTGCGTGAGc                                                                              <  3:265842/138‑1 (MQ=255)
                                                                             ggAACTCTCAACATTATAAAATACAACCGACAGCGAGCCCGGATTATCCCAGATGGCGAGGGGAACTCAGCGTGTTCGGCATAACCTTTGCCTAAGACAGCGGATTATCCCCGCTGGCCGAGGGAACACTGCGTGAgg                                                                              >  4:265842/1‑137 (MQ=255)
                                                                                    tcAACATTATCAATTACAACCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGgcgc                                                                                                <  4:277447/113‑1 (MQ=255)
                                                                                    tcAACATTATCAATTACAACCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGgcgc                                                                                                >  3:277447/1‑113 (MQ=255)
                                                                                    tcAACATTATCAATTACAACCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGCGTGAGCGTATCGCCGCGCGTCt                                                              <  5:90552/147‑1 (MQ=255)
                                                                                    tcAACATTATCAATTACAACCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGCGTGAGCGTATAGCCGCGCGTCTGc                                                            >  1:307964/1‑149 (MQ=255)
                                                                                    tcAACATTATCAATTACAACCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTAAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGCGTGAGCGTATCGCCGCGCGTCTGc                                                            >  5:25310/1‑149 (MQ=255)
                                                                                    tcAACATTATCAATTACAACCGAAAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGCGTGAGCGTATCGCCGCGCGTCTGc                                                            >  5:317840/1‑149 (MQ=255)
                                                                                      aaCATTATCAATTACAACCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCACCGTGTTCGGCATCACCtt                                                                                                                                >  3:338403/1‑79 (MQ=255)
                                                                                         aTTATCAATTACAACCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGTACTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGCGTGAGCGTATCGCCGCGCGTCTGCGAAAg                                                       <  8:305533/149‑1 (MQ=255)
                                                                                             tCAATTACAACCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGCGTGAGCGTATCGCCGCGCGTCTGCGAAAGCGGt                                                   <  2:217040/149‑1 (MQ=255)
                                                                                             tCAATTACAACCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGCGTGAGCGTATCGCCGCGCGTCTGCGAAAGCGGt                                                   <  2:307964/149‑1 (MQ=255)
                                                                                             tCAATTACAACCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCTGTGGCTGCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGCGTGAGCGTATCGCCGCGCGTCTGCGAAAGCGGt                                                   <  6:134809/149‑1 (MQ=255)
                                                                                                aTTACAACCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGCGTGAGCGTATCGCCGCGCGTCTGCGAAAGCGGTTTa                                                >  8:181197/1‑149 (MQ=255)
                                                                                                 ttACAACCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCGTGTGCGGAATCACCTTTGGCTTCGGCTGCGGCTTATCCCCGCTGGCGCGGGGAActct                                                                                      <  7:147927/110‑1 (MQ=255)
                                                                                                 ttACAACCGACAGGGAGCCCGGTTTATCCACGCTGGCGCGGGGAACTCAGCGAGTTCGGCATCACCTATGGCTTCGGGTGCGGTTTATCCCCGCTGGCGCGGGGAActct                                                                                      >  8:147927/1‑110 (MQ=38)
                                                                                                     aaCCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGTGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGCGTGAGCGTATCGCCGCGCGTCTGCGAAAGCGGTTTATcccc                                           >  4:318327/1‑149 (MQ=255)
                                                                                                     aaCCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCTTTTGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGAGTGAGCGTATCGCCGCGCGTCTGCGAAAGCGGTTTATcccc                                           >  6:254470/1‑149 (MQ=255)
                                                                                                     aaCCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGCGTGAGCGTATCGCCGCGCGTCTGCGAAAGCGGTTTATcccc                                           >  2:325582/1‑149 (MQ=255)
                                                                                                     aaCCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGCGTGAGCGTATCGCCGCGCGTCTGCGAAAGCGGTTTATcccc                                           >  6:309551/1‑149 (MQ=255)
                                                                                                     aaCCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGAGTGTGCGGAAACAGCGTAGGCTTCGGCTGAGCTTTATCAACGGTGGCGCGGGGACCTCTCAGTGAGCGTATAGGTGGGCGGCTGCGAAAGCGGTTTATca                                             >  6:1310/1‑146 (MQ=255)
                                                                                                     aaCCGACAGGGAGCACGGTTTATCCCCGCTGGCGCGGGGAACTCAGTGTGTTCGGCATCACCTTTGGCTTCGGATTCGGTTTATCCCCGCTGGCGCGGGGAACTCTCAGTGAACGTATCGCCGCGAGTCTGCGAAATCGGTTTATccca                                           >  4:303302/1‑148 (MQ=38)
                                                                                                      aCCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGCGTGAGCGTATCGCCGCGCGTCTGCGAAAGCGGTTTATCCCCg                                          >  5:82259/1‑149 (MQ=255)
                                                                                                       ccGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGCGTGAGCGTATCGCCGCGCGTCTGCGAAAGCGGTTTATCCCCGc                                         <  7:181197/149‑1 (MQ=255)
                                                                                                              ggAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGCGTGAGCGTATCGCCGCGCGTCTGCGAAAGCGGTTTATCCCCGCTGGcgcg                                  >  4:344115/1‑149 (MQ=255)
                                                                                                                aGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGTGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGCGTGAGCGTATCGCCGCGCGTCTGCGAAAGCGGTTTATCCCCGCTGGCGCggg                                <  3:303302/149‑1 (MQ=255)
                                                                                                                aGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGTGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGCGTGAGCGTATCGCCGCGCGTCTGCGAAAGCGGTTTATCCCCGCTGGCGCggg                                <  3:318327/149‑1 (MQ=255)
                                                                                                                aGCCCGGTCTCTCCCCGCTGGCGCGGGGAACTCAGCGTGTTCGGCAGCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGCGTGAGCGTATCGCCGCGCGTCTGCGAAAGCGGTTTATCCCCGCTGGCGCggg                                <  1:277772/149‑1 (MQ=255)
                                                                                                                 gCCCGGTTTATCCCCGCTGGCGCGGGGACCTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGCGTGAGCGTATCGCCGCGCGGCTGCGAAAGCGGTTTATCCCCGCTGGCGCggg                                >  5:213321/1‑148 (MQ=255)
                                                                                                                  cccGGTTTATCCCCGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGCGTGAGCGTATCGCCGCGCGTCTGCGAAAGCGGTTTATCCCCGCTGGCGCgggg                               >  6:178786/1‑148 (MQ=255)
                                                                                                                       tttATCCCCGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGCGTGAGCGTATCGCCGCGCGTCTGCGAAAGCGGTTTATCCCCGCTGGCGCGGGGAActct                         <  1:275681/149‑1 (MQ=255)
                                                                                                                            ccccGCTGGCGCGGGGAACTCAGCGTGTTCGGCATACCCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGCGTGAGCGTATCGCCGCGCGTCTGCGAAAGCGGTTTATCCCCGCTGGCGCGGGGAACTCTCTaa                     >  2:93337/1‑148 (MQ=255)
                                                                                                                             cccGCTGGCGCGGGGAACTCAGCGTGTTCGGCTTCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGCGTGCGCGTATCGCCGCGCGTCTGCGATAGCGGTTTATCCCCGCTGGCGCGGGGAACTCTCTaaaa                   <  8:212748/149‑1 (MQ=255)
                                                                                                                             cccGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGCGTGAGCGTATCGCCGCGCGTCTGCGAAAGCGGTTTATCCCCGCTGGCGCGGGGAACTCTCTaaaa                   <  2:347090/149‑1 (MQ=255)
                                                                                                                             cccGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGCGTGAGCGTATCGCCGCGCGTCTGCGAAAGCGGTTTATCCCCGCTGGCGCGGGGAACTCTCTaaaa                   <  7:219108/149‑1 (MQ=255)
                                                                                                                             cccGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGCGTGAGCGTATCGCCGCGCGTCTGCGAAAGCGGTTTATCCCCGCTGGCGCGGGGAACTCTCTaaaa                   <  7:73537/149‑1 (MQ=255)
                                                                                                                             cccGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGCGTGAGCGTATCGCCGCGCGTCTGCGAAAGCGGTTTATCCCCGCTGGCGCGGGGAACTCTCTaaaa                   <  8:100532/149‑1 (MQ=255)
                                                                                                                             cccGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGCGTGAGCGTATCGCCGCGCGTCTGCGAAAGCGGTTTATCCCCGCTGGCGCGGGGAACTCTCTaaaa                   <  6:82259/149‑1 (MQ=255)
                                                                                                                             cccGCTGGCGCGGGGAACTCAGCGTGTGCGGCAGCCCCGATCGCTTCCGCTGCGGTTTATCCCCGCGGGCGCGGGGAACTCTGCGTGAGCTTATCGCCGCGAGTCTGAGAAAGCGGTTTATCCCCGCTGGCGCGGGGAACTCTCTaaaa                   <  5:47455/149‑1 (MQ=255)
                                                                                                                              ccGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGCGTGAGCGTATCGCCGCGCGTCTGCGAAAGCGGTTTATCCCCGCTGGCGCGGGGAACTCTCTAAAAg                  <  7:167076/149‑1 (MQ=255)
                                                                                                                                    gcgcGGGGAACTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGCGTGAGCGTATCGCCGCGCGTCTGCGAAAGCGGTTTATCCCCGCTGGCGCGGGGAACTCTCTAAAAGtata              <  1:325582/147‑1 (MQ=255)
                                                                                                                                       cgGGGAACTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGCGTGAGCGTATCGCCGCGCGTCTGCGAAAGCGGTTTATCCCCGCTGGCGCGGGGAACTCTCTAAAAGTATACAttt         <  8:364630/149‑1 (MQ=255)
                                                                                                                                       cgGGGAACTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGCGTGAGCGTATCGCCGCGCGTCTGCGAAAGCGGTTTATCCCCGCTGGCGCGGGGAACTCTCTAAAAGTATACAttt         <  6:317840/149‑1 (MQ=255)
                                                                                                                                          ggAACTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGCGTGAGCGTATCGCCGCGCGTCTGCGAAAGCGGTTTATCCCCGCTGGCGCGGGGAACTCTCTAAAAGTATACATTTGtt      <  6:61689/149‑1 (MQ=255)
                                                                                                                                          ggAACTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGCGTGAGCGTATCGCCGCGCGTCTGCGAAAGCGGTTTATCCCCGCTGGCGCGGGGAACTCTCTAAAAGTATACATTTGtt      >  5:61689/1‑149 (MQ=255)
                                                                                                                                           gAACTCAGCGTGTTCGTCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGTGGAACTCTGCGTGAGCGTATCGCTGCGCGTCTGCGAAAGCGGTTTATCCCCGCTGGCGATGGGAACTCTCTAAAAGTATACATTTGTTc     <  5:1310/149‑1 (MQ=255)
                                                                                                                                              cTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGCGTGAGCGTATCGCCGCGCGt                                                                <  8:132297/87‑1 (MQ=255)
                                                                                                                                              cTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGCGTGAGCGTATCGCCGCGCGt                                                                >  7:132297/1‑87 (MQ=255)
                                                                                                                                               tCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGCGTGAGCGTATCGCCGCGCGTCTGCGAAAGCGGTTTATCCCCGCTGgc                                     >  8:240004/1‑113 (MQ=255)
                                                                                                                                               tCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGCGTGAGCGTATCGCCGCGCGTCTGCGAAAGCGGTTTATCCCCGCTGgc                                     <  7:240004/113‑1 (MQ=255)
                                                                                                                                               tCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGCGTGAGCGTATCGCCGCGCGTCTGCGAAAGCGGTTTATCCCCGCTGGCGCGGGGAACTCTCTAAAAGTATACATTTGTTCTTa  >  8:105960/1‑148 (MQ=255)
                                                                                                                                               tCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGCGTGAGCGTATCGCCGCGCGTCTGCGAAAGCGGTTTATCCCCGCTGGCGCGGGGAACTCTCTAAAAGTATACATTTGTTCTTa  >  6:115697/1‑148 (MQ=255)
                                                                                                                                               tCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGCGTGAGCGTATCGCCGCGCGTCTGCGAAAGCGGTTTATCCCCGCTGGCGCGGGGAACTCTCTAAAAGTATACATTTGTTCTTa  >  7:283106/1‑148 (MQ=255)
                                                                                                                                               tCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGCGTGAGCGTATCGCCGCGCGTCTGCGAAAGCGGTTTATCCCCGCTGGCGCGGGGAACTCTCTAAAAGTATACATTTGTTCTTa  >  3:191593/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
TCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACATTATCAATTACAACCGACAGGGAGCCCGGTTTATCCCCGCTGGCGCGGGGAACTCAGCGTGTTCGGCATCACCTTTGGCTTCGGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCTGCGTGAGCGTATCGCCGCGCGTCTGCGAAAGCGGTTTATCCCCGCTGGCGCGGGGAACTCTCTAAAAGTATACATTTGTTCTTA  >  NC_000913/2878197‑2878487

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: