Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,443,811 T→C 6.0% intergenic (+41/+80) mnmC → / ← yfcL fused 5‑methylaminomethyl‑2‑thiouridine‑forming enzyme methyltransferase and FAD‑dependent demodification enzyme/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,443,8110TC 6.0% 205.9 / 3.2 66intergenic (+41/+80)mnmC/yfcLfused 5‑methylaminomethyl‑2‑thiouridine‑forming enzyme methyltransferase and FAD‑dependent demodification enzyme/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (32/30);  new base C (2/2);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

AACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGTGTGACGTGTAGGCCTGATAAGACGCGATAAGCGTCGCATCAGGCATGCTGCTCCAGACGCCGCAACTTACTGCTGTGACGCCTGCTGAAACAGGTTTTCCCACATATCGGTCACCAGCGCCTGGTCACGCGGC  >  NC_000913/2443673‑2443957
                                                                                                                                          |                                                                                                                                                  
aaCCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGAcc                                                                                                                                            >  5:205874/1‑147 (MQ=255)
   cGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACg                                                                                                                                        <  7:303171/148‑1 (MQ=255)
     aTTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACgtg                                                                                                                                      >  4:173531/1‑148 (MQ=255)
      ttCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACgtgt                                                                                                                                     <  1:205241/148‑1 (MQ=255)
        ccGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGtgtg                                                                                                                                    <  2:27236/147‑1 (MQ=255)
         cGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGGTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGTGTGa                                                                                                                                   <  8:66855/147‑1 (MQ=255)
           atggatgCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGTGTGACgt                                                                                                                                >  8:97024/1‑148 (MQ=255)
            tggatgCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGTGTGACgtg                                                                                                                               <  7:97024/148‑1 (MQ=255)
              gatgCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGTAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGTGTGACgtgt                                                                                                                              <  3:173531/147‑1 (MQ=255)
              gatgCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGTGTGACgtgt                                                                                                                              <  6:368612/147‑1 (MQ=255)
                  ccAGTAAGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAACTGTTGAAGGGTACAGCGGTTAATGCGGTGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGTGTGACGTGGAGGc                                                                                                                          <  3:221436/147‑1 (MQ=39)
                      tACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGa                                                                                                                                              >  5:387172/1‑123 (MQ=255)
                      tACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGa                                                                                                                                              <  6:387172/123‑1 (MQ=255)
                       aCGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGTGTGACgtgt                                                                                                                              <  6:316914/138‑1 (MQ=255)
                       aCGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGTGTGACgtgt                                                                                                                              >  5:316914/1‑138 (MQ=255)
                            ggcggcGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGTGTGACGTGTAGGCCTGATAAGAcg                                                                                                               >  1:241222/1‑148 (MQ=255)
                             gcggcgTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGTGTGACGTGTAGGCCTGATAAGAcg                                                                                                               <  2:198018/147‑1 (MQ=255)
                                gcgTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGTGTGACGTGTAGGCCTGATAAGACGCGa                                                                                                            >  4:176107/1‑147 (MQ=255)
                                 cgTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGTGTGACGTGTAGGCCTGATGCGAcgcttat                                                                                                          <  2:147005/148‑5 (MQ=255)
                                   ttAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGTGTGACGTGTAGGCCTGATAAGACGCGATAAg                                                                                                        >  4:286528/1‑148 (MQ=255)
                                      aaCCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGTGTGACGTGTAGGCCTGATAAGACGCGATAAGCGt                                                                                                     <  3:176107/148‑1 (MQ=255)
                                          cGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGTGTGACGTGTAGGCCTGATAAGACGCGATAAGCGTCGc                                                                                                  >  1:280177/1‑147 (MQ=255)
                                          cGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGTGTGACGTGTAGGCCTGATAAGACGCGATAAGCGTCGCa                                                                                                 >  1:218112/1‑148 (MQ=255)
                                          cGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTCTTCCGACCTACGTGTGACGTGTAGGCCTGATAAGACGCGATAAGCGTCGCa                                                                                                 >  8:251085/1‑148 (MQ=255)
                                             aTCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGTGTGACGTGTAGGCCTGATAAGACGCGATAAGCGTCGCATCa                                                                                              <  1:161327/148‑1 (MQ=255)
                                                ggTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGTGTGACGTGTAGGCCTGATAAGACGCGATAAGCGTCGCATCAgg                                                                                            <  3:390203/147‑1 (MQ=255)
                                                  ttATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGTGTGACGTGTAGGCCTGATAAGACGCGATAAGCGTCGCATCAGGCAt                                                                                         >  5:257036/1‑148 (MQ=255)
                                                  ttATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGTGTGACGTGTAGGCCTGATAAGACGCGATAAGCGTCGCATCAGGCAt                                                                                         >  1:382778/1‑148 (MQ=255)
                                                         tGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGTGTGACGTGTAGGCCTGATAAGACGCGATAAGCGTCGCATCAGGCATGTTGCTc                                                                                  <  2:241222/148‑1 (MQ=255)
                                                         tGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGGGAACTCCTTATCCGACCTACGTGTGACGTGTAGGCCTGATAGGCCGCGATCAGCGTCGCAGCAGGCAtgctgagc                                                                                  >  2:440724/1‑145 (MQ=255)
                                                          gCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGTGTGACgtgt                                                                                                                              >  8:121902/1‑103 (MQ=255)
                                                          gCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGTGTGACgtgt                                                                                                                              <  7:121902/103‑1 (MQ=255)
                                                           cGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGTGTGACGTGTAGGCCTGATAAGACGCGATAAGCGTCGCATCAGGCATGCTGCTc                                                                                  >  8:284084/1‑146 (MQ=255)
                                                              aaaTTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGTGTGACGTGTAGGCCTGATAAGACGCGATAAGCGTCGCATCAGGCATGCTGCTCCAGAc                                                                             >  3:158577/1‑148 (MQ=255)
                                                                aTTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGTGTGACGTGTAGGCCTGATAAGACGCGATAAGCGTCGCATCAGGCATGCTGCTCCAGAcgc                                                                           <  7:284084/148‑1 (MQ=255)
                                                                 ttgttgAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGTGTGACGTGTAGGCCTGATAAGACGCGATAAGCGTCGCATCAGGCATGCTGCTCCAGAcgcc                                                                          <  6:205874/148‑1 (MQ=255)
                                                                  tgttgAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTACCAGTGTATGCCGGATCCGGCGTGAACGCCTTAACCGACCTACGTGTGACTTGTAGGCCTCATAAGACGCGATAAGCGTCTCATCAGGCATGCTGCTCCAGAcgccg                                                                         >  7:196217/1‑148 (MQ=255)
                                                                     tgAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGTGTGACGTGTAGGCCTGATAAGACGCGATAAGCGTCGCATCAGGCATGCTGCTCCAGACGCCGCaa                                                                      >  7:372987/1‑148 (MQ=255)
                                                                      gAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGTGTGACGTGTAGGCCTGATAAGACGCGATAAGCGTCGCATCAGGCATGCTGCTCCAGACGCCGCAAc                                                                     >  4:193537/1‑148 (MQ=255)
                                                                      gAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGTGTGACGTGTAGGCCTGATAAGACGCGATAAGCGTCGCATCAGGCATGCTGCTCCAGACGCCGCAAc                                                                     >  8:249224/1‑148 (MQ=255)
                                                                        aGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGTGTGACGTGTAGGCCTGATAAGACGCGATAAGCGTCGCATCAGGCATGCTGCTCCAGACGCCGCAACtt                                                                   <  3:286528/148‑1 (MQ=255)
                                                                               aGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGGTTGCCGGATGCGGCGTGCACTCATTATCCGACGTAGGTGTGACGTGTAGGCCTGATAAGACGCGAAAAGCGTCGCATAAGGCATGCTGCTCCAGACGCCGCAACTTActgctg                                                            <  5:386033/148‑1 (MQ=255)
                                                                                 cGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGTGTGACGTGTAGGCCTGATAAGACGCGATAAGCGTCGCATCAGGCATGCTGCTCCAGACGCCGCAACTTACTGCtgtg                                                          <  2:382778/148‑1 (MQ=255)
                                                                                 cGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGTGTGACGTGTAGGCCTGATAAGACGCGATAAGCGTCGCATCAGGCATGCTGCTCCAGACGCCGCAACTTACTGCtgtg                                                          <  4:158577/148‑1 (MQ=255)
                                                                                    ttAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGTGTGACGTGTAGGCCTGATAAGACGCGATAAGCGTCGCATCAGGCATGCTGCTCCAGACGCCGCAACTTACTGCTGTGACg                                                       >  1:334563/1‑148 (MQ=255)
                                                                                               tAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGTGTGACGTGTAGGCCTGATAAGACGCGATAAGCGTCGCATCAGGCATGCTGCTCCAGACGCCGCAACTTACTGCTGTGACGCCTGCTGAAAc                                            <  6:257036/148‑1 (MQ=255)
                                                                                                 aTCGGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCGACGTGTGACGTGTAGGCCTGCTAAGACGCGATAAGCGTCGCATCAGGCATGCGGCTCCAGACGCCGCAACTTACTGCTGTGACGCCTGCTGAAACa                                           <  6:106514/147‑1 (MQ=255)
                                                                                                  tCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGTGTGACGTGTAGGCCTGATAAGACGCTATAAGCGTCGCATCAGGCATGCTGCTCCAGACGCCGCAACTTACTGCTGTGCCGCCTGCTGTAACAgg                                         >  4:340077/1‑148 (MQ=255)
                                                                                                   cTGCACTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGCCCTACGTGTGCCGTGTAGGCCTGATAAGACGCGATAAGCGTCGCATCAGGCATGCTGCTCCAGACGCCGCAACTTACTGCTGTGACGCCTGCTGAAACAGGt                                        >  4:312921/1‑148 (MQ=255)
                                                                                                     gCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGTGTGACGTGTAGGCCTGATAAg                                                                                                                  >  7:9990/1‑72 (MQ=255)
                                                                                                     gCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGTGTGACGTGTAGGCCTGATAAg                                                                                                                  <  8:9990/72‑1 (MQ=255)
                                                                                                      ctctCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGTGTGACGTGTAGGCCTGATAAGACGCGATAAGCGTCGCATCAGGCATGCTGCTCCAGACGCCGCAACTTACTGCTGTGACGCCTGCTGAAACAgg                                         <  2:218112/144‑1 (MQ=255)
                                                                                                       tctcGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGTGTGACGTGTAGGCCTGATAAGACGCGATAAGCGTCGCATCAGGCATGCTGCTCCAGACGCCGCAACTTACTGCTGTGACGCCTGCTGAAACAGGTGTTc                                    >  2:414507/1‑148 (MQ=255)
                                                                                                        ctcGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGTGTGACGTGTAGGCCTGATAAGACGCGATAAGCGTCGCATCAGGCATGCTGCTCCAGACGCCGCAACTTACTGCTGTGACGCCTGCTGAAACAGGTTTTcc                                   >  3:84338/1‑148 (MQ=255)
                                                                                                         tcGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGTGTGACGTGTAGGCCTGATAAGACGCGATAAGCGTCGCATCAGGCATGCTGCTCCAGACGCCGCAACTTACTGCTGTGACGCCTGCTGAAACAGGTTTTc                                    <  3:314518/146‑1 (MQ=255)
                                                                                                           gCAGTGTTTGCCGGATGCGGCGTGAACGCCTCTTCCGACCTACGTGTGACGTGTAGGCCTGATAAGACGCGATAAGCGTCGCATCAGGCATGCTGCTCCAGACGCCGCAACTTACTGCTGTGACGCCTGCTGAAACAGGTTTTcc                                   <  7:251085/145‑1 (MQ=255)
                                                                                                                gtTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGTGTGACGTGTAGGCCTGATAAGACGCGATAAGCGTCTCATCAGGCATGCTGCTCCAGACGCCGCAACTTACTGCtgtg                                                          >  4:258104/1‑117 (MQ=255)
                                                                                                                gtTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGTGTGACGTGTAGGCCTGATAAGACGCGATAAGCGTCTCATCAGGCATGCTGCTCCAGACGCCGCAACTTACTGCtgtg                                                          <  3:258104/117‑1 (MQ=255)
                                                                                                                   tGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGTGTGACGTGTAGGCCTGATAAGACGCGATAAGCGTCGCATCAGGCATGCTGCTCCAGACGCCGCAACTTACTGCTGTGACGCCTGCTGAAACAGGTTTTCCCACATATCgg                         <  1:418091/147‑1 (MQ=255)
                                                                                                                      cGGATGCGGCGTGAACGCCTTATCCGACCTACGTGTGACGTGTAGGCCTGATAAGACGCGATAAGCGTCGCATCAGGCATGCTGCTCCAGACGCCGCAACTTACTGCTGTGACGCCTGCTGAAACAGGTTTTCCCACATATCGGTCAc                     >  1:23643/1‑148 (MQ=255)
                                                                                                                       ggATGCGGCGTGAACGCCTTATCCGACCTACGTGTGACGTGTAGGCCTGATAAGACGCGATAAGCGTCGCATCAGGCATGCTGCTCCAGACGCCGCATCTTACTGCTGTGACGCCTGCTGAAACAGGTTTTCCCACATATCGGTCAc                     >  2:103563/1‑147 (MQ=255)
                                                                                                                         ataaggcgTTCACGCCGCATCCGACCTACGTGTGACGTGTAGGCCTGATAAGACGCGATAAGCGTCGCATCAGGCATGCTGCTCCAGACGCCGCAACTTACTGCTGTGACGCCTGCTGAAACAGGTTTTCCCACATATCGGTCACCAg                  >  2:230260/5‑148 (MQ=39)
                                                                                                                              gcgttcaCGCCGCATCCGACCTACGTGTGACGTGTAGGCCTGATAAGACGCGATAAGCGTCGCATCAGGCATGCTGCTCCAGACGCCGCAACTTACTGCTGTGACGCCTGCTGAAACAGGTTTTCCCACATATCGGTCACCAGCGCCt             <  1:230260/142‑1 (MQ=255)
                                                                                                                              gcgTGAACGCCTTATCCGACCTACGTGTGACGTGTAGGCCTGATAAGACGCGATAAGCGTCGCATCAGGCATGCTGCTCCAGACGCCGCAACTTACTGCTGTGACGCCTGCTGAAACAGGTTTTCCCACATATCGGTCACCAGCGCCt             >  1:187994/1‑148 (MQ=255)
                                                                                                                                     cGCCTTATCCGACCTACGTGTGACGTGTAGGCCTGATAAGACGCGATAAGCGTCGCATCAGGCATGCTGCTCCAGACGCCGCAACTTACTGCTGTGCCGCCTGCTGAAACAGGTTTTCCCACATATCGGTCACCAGCGCCTGGTCAcg      >  5:130469/1‑148 (MQ=255)
                                                                                                                                      gCCTTATCCGACCTACGTGTGACGTGTAGGCCTGATAAGACGCGATAAGCGTCGCATCAGGCATGCTGCTCCAGACGCCGCAACTTACTGCTGTGACGCCTGCTGAAACAGGTTTTCCCACATATCGGTCACCAGCGCCTGGTCAcgc     >  2:129489/1‑148 (MQ=255)
                                                                                                                                          tATCCGACCTACGTGTGACGTGTAGGCCTGATAAGACGCGATAAGCGTCGCATCAGGCATGCTGCTCCAGACGCCGCAACTTACTGCTGTGACTCCTGCTGAAACAGGTTTTCCCACATATCGGTCACCAGCGCCTGGTCACgcggc  >  4:207736/1‑147 (MQ=255)
                                                                                                                                          tATCCGACCTACGTGTGACGTGTAGGCCTGATAAGACGCGATAAGCGTCGCATCAGGCATGCTGCTCCAGACGCCGCAACTTACTGCTGTGACGCCTGCTGAAACAGGTTTTCCCACATATCGGTCACCAGCGCCTGGTCACgcggc  >  2:367311/1‑147 (MQ=255)
                                                                                                                                          |                                                                                                                                                  
AACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAATCGGTTATGGGTGCGGAAATTGTTGAAGGGTAAAGCGGTTAAGGCGGGGTAATCTGCTCTCGCAGTGTTTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGTGTGACGTGTAGGCCTGATAAGACGCGATAAGCGTCGCATCAGGCATGCTGCTCCAGACGCCGCAACTTACTGCTGTGACGCCTGCTGAAACAGGTTTTCCCACATATCGGTCACCAGCGCCTGGTCACGCGGC  >  NC_000913/2443673‑2443957

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