Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 342,393 C→T 8.2% D423D (GAC→GAT yahJ → putative metallo‑dependent hydrolase domain deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913342,3930CT 8.2% 222.7 / 7.6 74D423D (GAC→GATyahJputative metallo‑dependent hydrolase domain deaminase
Reads supporting (aligned to +/- strand):  ref base C (34/33);  new base T (3/3);  total (38/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.75e-01

ACATGCTGGAAAAAGCCAATCTGTACGCGCAGCTCTATATTCGTCCTAACGAACAGAACCTCTCCCGCTCGCTGTTTTTAGCCACTGGCGATGTATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGCGCGTATCTCGCCGAGAACCGCAACGTTCCATAAAGGGCAACTGGTGTGGGGGAGTGTGGCAGGTTGATGCGCATTTGAAACGCCTGCAACGGTGAA  >  NC_000913/342248‑342536
                                                                                                                                                 |                                                                                                                                               
aCATGCTGGAAAAAGCCAATCTGTACGCGCAGCTCTATATTCGTCCTAACGAACAGAACCTCTCCCGCTCGCTGTTTTTAGCCACTGGCGATGTATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGc                                                                                                                                               <  6:19430/148‑1 (MQ=255)
    gCTGGAAAAAGCCAATCTGTACGCGCAGCTCTATATTCGTCCTAACGAACAGAACCTCTCCCGCTCGCTGTTTTTAGCCACTGGCGATGTATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGc                                                                                                                                           >  8:391577/1‑148 (MQ=255)
    gCTGGAAAAAGCCAATCTGTACGCGCAGCTCTATATTCGTCCTAACCAACAGACCCTCACCCGCTCGCTGTTTTTAGCCACTGGCGATGTATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGAAGCCAGc                                                                                                                                           >  7:122387/1‑148 (MQ=255)
      tGGAAAAAGCCAATCTGTACGCGCAGCTCTATATTCGTCCTAACGAACAGAACCTCTCCCGCTCGCTGTTTTTAGCCACTGGCGATGTATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCtt                                                                                                                                         >  8:27995/1‑148 (MQ=255)
         aaaaaGCCAATCTGTACGCGCAGCTCTATATTCGTCCTAACGAACAGAACCTCTCCCGCTCGCTGTTTTTAGCCACTGGCGATGTATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTtgt                                                                                                                                      <  8:340421/148‑1 (MQ=255)
            aaGCCAATCTGTACGCGCAGCTCTATATTCGTCCTAACGAACAGAACCTCTCCCGCTCGCTGTTTTTAGCCACTGGCGATGTATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCt                                                                                                                                   <  7:391577/148‑1 (MQ=255)
                cAATCTGTACGCGCAGCTCTATATTCGTCCTAACGAACAGAACCTCTCCCGCTCGCTGTTTTTAGCCACTGGCGATGTATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCtggtg                                                                                                                               >  4:309895/1‑148 (MQ=255)
                    cTGTACGCGCAGCTCTCTATTCGTCCTAACGAACAGAACCTCTCCCGCTCGCTGGTTTTCGCCACTGGCGATGTATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAACGCGCAGGATGACGCCAGCTTTGTGCTGGGGGACg                                                                                                                           <  3:5492/148‑1 (MQ=255)
                    cTGTACGCGCAGCTCTATATTCGTCCTAACGAACAGAACCTCTCCCGCTCGCTGTTTTTAGCCACTGGCGATGTATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGATGCCAGCTTTGTGCTGGTGGACg                                                                                                                           >  5:385185/1‑148 (MQ=255)
                    cTGTACGCGCAGCTCTATATTCGTCCTAACGAACAGAACCTCTCCCGCTCGCTGTTTTTAGCCACTGGCGATGTATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGATGCCAGCTTTGTGCTGGTGGACg                                                                                                                           >  5:19858/1‑148 (MQ=255)
                       tACGCGCAGCTCTATATTCGTCCTAACGAACAGAACCTCTCCCGCTCGCTGTTTTTAGCCACTGGCGATGTATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACg                                                                                                                           <  7:27995/145‑1 (MQ=255)
                                 tctATATTCGTCCTAACGAACAGAACCTCTCCCGCTCGCTGTTTTTAGCCACTGGCGATGTATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTccgc                                                                                                              >  8:16771/1‑148 (MQ=255)
                                    atatTCGTCCTAACGAACAGAACCTCTCCCGCTCGCTGTTTTTAGCCACTGGCGATGTATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTccgcc                                                                                                             >  5:153548/1‑146 (MQ=255)
                                        tCGTCCTAACGAACAGAACCTCTCCCGCTCGCTGTTTTTAGCCACTGGCGATGTATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTccgcc                                                                                                             <  8:154844/142‑1 (MQ=255)
                                        tCGTCCTAACGAACAGAACCTCTCCCGCTCGCTGTTTTTAGCCACTGGCGATGTATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTccgcc                                                                                                             >  7:154844/1‑142 (MQ=255)
                                               aaCGAACAGAACCTCTCCCGCTCGCTGTTTTTAGCCACTGGCGATGTATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGgcgc                                                                                                >  2:340254/1‑148 (MQ=255)
                                                aCGAACAGAACCTCTCCCGCTCGCTGTTTTTAGCCACTGGCGATGTATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGgcgc                                                                                                <  6:86650/147‑1 (MQ=255)
                                                 cGAACAGAACCTCTCCCGCTCGCTGTTTTTAGCCACTGGCGATGTATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGCGCGt                                                                                              <  7:16771/148‑1 (MQ=255)
                                                    aCAGAACCTCTCCCGCTCGCTGTTTTTAGCCACTGGCGATGTATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGTGCGTAtc                                                                                           >  6:89290/1‑148 (MQ=255)
                                                       gAACCTCTCCCGCTCGCTGTTTTTAGCCACTGGCGATGTATTGCCGCTTAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCtg                                                                                                                                  <  6:53608/106‑1 (MQ=255)
                                                       gAACCTCTCCCGCTCGCTGTTTTTAGCCACTGGCGATGTATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCtg                                                                                                                                  >  5:53608/1‑106 (MQ=255)
                                                       gAACCTCTCCCGCTCGCTGTTTTTAGCCACTGGCGATGTATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGCGCGTATCTCg                                                                                        <  6:153548/148‑1 (MQ=255)
                                                         aCCTCTCCCGCTCGCTGTTTTTAGCCACTGGCGATGTATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGCGCGTATCTCGcc                                                                                      >  7:393753/1‑148 (MQ=255)
                                                                ccGCTCGCTGTTTTTAGCCACTGGCGATGTATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTtgt                                                                                                                                      >  4:27718/1‑93 (MQ=255)
                                                                ccGCTCGCTGTTTTTAGCCACTGGCGATGTATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTtgt                                                                                                                                      <  3:27718/93‑1 (MQ=255)
                                                                ccGCTCGCTGTTTTTAGCCACTGGCGATGTATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGCGCGTATCTCGCCGAGaa                                                                                 <  1:340254/146‑1 (MQ=255)
                                                                ccGCTCGCTGTTTTTAGCCACTGGCGATGTATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGCGCGTATCTCGCCGAGAAc                                                                                <  5:384661/147‑1 (MQ=255)
                                                                       ctGTTTTTAGCCACTGGCGATGTATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGCGCGTATCTCGCCGAGAACCGCAACGt                                                                        >  4:268354/1‑148 (MQ=255)
                                                                         gTTTTTAGCCACTGGCGATGTATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCg                                                                                                       <  1:203624/115‑1 (MQ=255)
                                                                         gTTTTTAGCCACTGGCGATGTATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCg                                                                                                       >  2:203624/1‑115 (MQ=255)
                                                                          tttttAGCCACTGGCGATGTATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGCGCGTATCTCGCCGAGAACCGCAACGTTc                                                                      >  8:2470/1‑147 (MQ=255)
                                                                               aGCCACTGGCGATGTATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGCGCGTATCTCGCCGAGAACCGCAACGTTCCATaa                                                                 <  7:2470/147‑1 (MQ=255)
                                                                               aGCCACTGGCGATGTATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGCGCGTATCTCGCCGAGAACCGCAACGTTCCATaa                                                                 <  3:268354/147‑1 (MQ=255)
                                                                                    cTGGCGATGTATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGCGCGTATCTCGCCGAGAACCGCAACGTTCCATAAAGGGCa                                                           <  3:309895/148‑1 (MQ=255)
                                                                                     tGGCGATGTATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGATGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGCGCGTATCTCGCCGAGAACCGCAACGTTCCATAAAGGGCaa                                                          >  5:241475/1‑148 (MQ=255)
                                                                                     tGGCGATGTATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGCGCGTATCTCGCCGAGAACCGCAACGTTCCATAAAGGGCaa                                                          >  7:80932/1‑148 (MQ=255)
                                                                                     tGGCGATGTATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGCGCGTATCTCGCCGAGAACCGCAACGTTCCATAAAGGGCaa                                                          >  6:300548/1‑148 (MQ=255)
                                                                                     cggCGATGTATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGCGCGTATCTCGCCGAGAACCGCAACGTTCCa                                                                    <  3:119439/137‑1 (MQ=255)
                                                                                     cggCGATGTATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGCGCGTATCTCGCCGAGAACCGCAACGTTCCa                                                                    >  4:119439/2‑138 (MQ=255)
                                                                                      ggCTATGTATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGCGCGTATCTCGCCGAGAACCGCAACGTTCCATAAAGGGCAAc                                                         >  5:360085/1‑148 (MQ=255)
                                                                                      ggCGATGTATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGCGCGTATCTCGCCGAGAACCGCAACGTTCCATAAAGGGCAAc                                                         >  1:62709/1‑148 (MQ=255)
                                                                                      ggCGATGTATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGCGCGTATCTCGCCGAGAACCGCAACGTTCCATAAAGGGCAAc                                                         >  3:223881/1‑148 (MQ=255)
                                                                                      ggCGATGTATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGCGCGTATCTCGCCGAGAACCGCAACGTTCCATAAAGGGCAAc                                                         >  4:81349/1‑148 (MQ=255)
                                                                                      ggCGATGTATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGCGCGTATCTCGCCGAGAACCGCAACGTTCCATAAAGGGCAAc                                                         >  3:346657/1‑148 (MQ=255)
                                                                                          aTGTATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGCGCGTATCTCGCCGAGAACCGCAACGTTCCATAAAGGGCAACTGgt                                                     <  8:393753/148‑1 (MQ=255)
                                                                                          aTGTATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGCGCGTATCTCGCCGAGAACCGCAACGTTCCATAAAGGGCAACTGgt                                                     <  8:80932/148‑1 (MQ=255)
                                                                                             tATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGATGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGCGCGTATCTCGCCGAGAACCGCAACGTTCCATAAAGGGCAACTGgtgt                                                   <  6:19858/147‑1 (MQ=255)
                                                                                             tATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGATGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGCGCGTATCTCGCCGAGAACCGCAACGTTCCATAAAGGGCAACTGgtgt                                                   <  6:385185/147‑1 (MQ=255)
                                                                                             tATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGTGCGTATCTCGCCGAGAACCGCAACGTTCCATAAAGGGCAACTGGtgtg                                                  <  5:89290/148‑1 (MQ=255)
                                                                                             tATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGCGCGTATCTCGCCGAGAACCGCAACGTTCCATAAAGGGCAACTGgtgt                                                   <  3:81349/147‑1 (MQ=255)
                                                                                              aTTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGCGCGTAtc                                                                                           <  1:335722/106‑1 (MQ=255)
                                                                                              aTTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGCGCGTAtc                                                                                           >  2:335722/1‑106 (MQ=255)
                                                                                                  ccGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGCGCGTATCTCGCCGAGAACCGCAACGTTCCATAAAGGGCAACTGGTGTGGGGGa                                             <  4:223881/148‑1 (MQ=255)
                                                                                                       gaatgaaAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGCGCGTATCTCGCCGAGAACCGCAACGTTCCATAAAGGGCAACTGGTGTGGGGGAGtgtg                                        >  5:346317/1‑148 (MQ=255)
                                                                                                            aaaaaGGCGAGCGTGTATGGACAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGTACGCCTCCTGTTCCGCCGAGGCGGTGGCGCGTATCTCGCCGAGCACCGCAACGTTCCATAAAGGGCAACTGGTGTGGGGGAGTGTTGCAgg                                   >  2:322801/1‑148 (MQ=255)
                                                                                                               aaGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGCGCGTATCTCGCCGAGATCCGCAACGTTCCATAAAGGGCAACTGGTGTGGGGGAGTGTGGCAGGTTg                                <  4:346657/148‑1 (MQ=255)
                                                                                                                 ggCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGCGCGTATCTCGCCGAGAACCGCAACGTTCCATAAAGGGCAACTGGTGTGGGGGAGTGTGGCAGGTTGAt                              <  5:300548/148‑1 (MQ=255)
                                                                                                                  gCGAGCGTGTATGGTCAAAAGCGCAGGATGTCGCCAGCTTTGTGCTGGTGGACGCCTCTTGTTCCGCCGAGGCGGTGGCGCGTATCTCGCCGAGATCCGCAACGTTCCATAAAGGGAATCTGGTGTGGGGGAGTGTGGCAGGTTgtt                              <  8:122387/147‑3 (MQ=255)
                                                                                                                     aGCGTGTATGGCCAACAGCGCAGGATGACCCCAGCTTTGTGCTGGCGGACGCCTCCGGTACCGCCGCGGCGGTGGCGCGTATCTCGCCGAGAACCGAAACTTTCCAAAAAGGGCAACTGGTGTGGGGTAGTTTCCCAGGTTGATccgc                          <  6:360085/148‑1 (MQ=255)
                                                                                                                      gCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGCGCGTATCTCGCCGAGAACCGCAACGTTCCATAAAGGGCAACTGgtgt                                                   >  4:312817/1‑122 (MQ=255)
                                                                                                                      gCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGCGCGTATCTCGCCGAGAACCGCAACGTTCCATAAAGGGCAACTGgtgt                                                   <  3:312817/122‑1 (MQ=255)
                                                                                                                      gCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGCGCGTATCTCGCCGAGAACCGCAACGTTCCATAAAGGGCAACTGGTGTGGGGGAGTGTGGCAGGTTGATGCGCa                         >  5:13736/1‑148 (MQ=255)
                                                                                                                       cGTGTATGGCCAAAAGCGCAGGATGATGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGCGCGTATCTCGCCGAGAACCGCAACGTTCCATAAAGGGCAACTGGTGTGGGGGAGTGTGGCAGGTTGATGCGCAt                        <  6:241475/148‑1 (MQ=255)
                                                                                                                          gtATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGCGCGTAt                                                                                            >  2:152033/1‑77 (MQ=255)
                                                                                                                          gtATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGCGCGTAt                                                                                            <  1:152033/77‑1 (MQ=255)
                                                                                                                            aTGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGCGCGTATCTCGCCGAGAACCGCAACGTTCCATAAAGGGCAACTGGTGTGGGGGAGTGTGGCAGGTTGATGCGCATTTGa                    >  1:306079/1‑147 (MQ=255)
                                                                                                                             tGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGCGCGTATCTCGCCGAGAACCGCAACGTTCCATAAAGGGCAACTGGTGTGGGGGAGTGTGGCAGGTTGATGCGCATTTGaa                   <  2:62709/147‑1 (MQ=255)
                                                                                                                                ccAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGCGCGTATCTCGCCGAGAACCGCAACGTTCCATAAAGGGCa                                                           <  4:266052/104‑1 (MQ=255)
                                                                                                                                ccAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGCGCGTATCTCGCCGAGAACCGCAACGTTCCATAAAGGGCa                                                           >  3:266052/1‑104 (MQ=255)
                                                                                                                                             aTGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGCGCGTATCTCGCCGAGAACCGCAACGTTCCATAAAGGGCAACTGGTGTGGGGGAGTGTGGCAGGTTGATGCGCATTTGAAACGCCTGCAACGGTGaa  >  5:125255/1‑148 (MQ=255)
                                                                                                                                             aTGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGCGCGTATCTCGCCGAGAACCGCAACGTTCCATAAAGGGCAACTGGTGTGGGGGAGTGTGGCAGGTTGATGCGCATTTGAAACGCCTGCAACGGTGaa  >  6:83855/1‑148 (MQ=255)
                                                                                                                                              tGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGCGCGTATCTCGCCGAGAACCGCAACGtt                                                                       <  5:94411/78‑1 (MQ=255)
                                                                                                                                              tGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGCGCGTATCTCGCCGAGAACCGCAACGtt                                                                       >  6:94411/1‑78 (MQ=255)
                                                                                                                                              tGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGCGCGTATCTCGCCGAGAACCGCAACGTTCCATAAAGGGCAACTGGTGTGGGGGAGTGTGGCAGGTTGATGCGCATTTGAAACGCCTGCAACGGTGaa  >  5:128222/1‑147 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
ACATGCTGGAAAAAGCCAATCTGTACGCGCAGCTCTATATTCGTCCTAACGAACAGAACCTCTCCCGCTCGCTGTTTTTAGCCACTGGCGATGTATTGCCGCTGAATGAAAAAGGCGAGCGTGTATGGCCAAAAGCGCAGGATGACGCCAGCTTTGTGCTGGTGGACGCCTCCTGTTCCGCCGAGGCGGTGGCGCGTATCTCGCCGAGAACCGCAACGTTCCATAAAGGGCAACTGGTGTGGGGGAGTGTGGCAGGTTGATGCGCATTTGAAACGCCTGCAACGGTGAA  >  NC_000913/342248‑342536

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: