Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,402,942 T→A 5.7% N105I (AAT→ATT)  abgB ← p‑aminobenzoyl‑glutamate hydrolase, B subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,402,9420TA 5.7% 215.0 / 2.8 70N105I (AAT→ATT) abgBp‑aminobenzoyl‑glutamate hydrolase, B subunit
Reads supporting (aligned to +/- strand):  ref base T (27/38);  new base A (2/2);  total (29/41)
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.92e-01

AACGTTTTACCCGAGCCGCCTTCTTCGCCAGGACAACCATAAAAGCGCACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTG  >  NC_000913/1402795‑1403088
                                                                                                                                                   |                                                                                                                                                  
aaCGTTTTACCCGAGCCGCCTTCTTCGCCAGGACAACCATAAAAGCGCACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCAt                                                                                                                                                    >  5:169084/1‑148 (MQ=255)
aaCGTTTTACCCGAGCCGCCTTCTTCGCCAGGACAACCATAAAAGCGCACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCAt                                                                                                                                                    >  7:28539/1‑148 (MQ=255)
aaCGTTTTACCCGAGCCGCCTTCTTCGCCAGGACAACCATAAAAGCGCACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCAt                                                                                                                                                    <  8:300383/148‑1 (MQ=255)
    ttttACACGAGGCGCCTTCTTCGCCAGGACAACCATAAAAGCGCACCGTGACGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGGTATTGCTGCGGCAAATGCGGCGGTTCGCAGCAAATTGTGTCCGCAACCGTGACCATTTTc                                                                                                                                                <  8:115011/148‑1 (MQ=255)
       tACCCGAGCCGCCTTCTTCGCCAGGACAACCATAAAAGCGCACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCAtt                                                                                                                                                   >  3:12307/1‑142 (MQ=255)
       tACCCGAGCCGCCTTCTTCGCCAGGACAACCATAAAAGCGCACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCAtt                                                                                                                                                   <  4:12307/142‑1 (MQ=255)
        aCCCGAGCCGCCTTCTTCGCCAGGACAACCATAAAAGCGCACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGAAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCAcc                                                                                                                                             <  5:363332/147‑1 (MQ=255)
        aCCCGAGCCGCCTTCTTCGCCAGGACAACCATAAAAGCGCACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTAGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCAATTTCACCg                                                                                                                                            >  5:262607/1‑148 (MQ=255)
        aCCCGAGCCGCCTTCTTCGCCAGGACAACCATAAAAGCGCACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTAGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCAATTTCACCg                                                                                                                                            >  5:257070/1‑148 (MQ=255)
        aCCCGAGCCGCCTTATTCGCAAGGACAACCATAAAAGCGCACCTTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCTGCAAAGTCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCAcc                                                                                                                                             <  5:103152/147‑1 (MQ=255)
              gccgccTTCTTCGCCAGGACAACCATAAAAGCGCACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGt                                                                                                                                       >  5:204858/1‑147 (MQ=255)
                  ccTTCTTCGCCAGGACAACCATAAAAGCGCACCGTGCCGCCTTGCCGATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACgg                                                                                                                                  <  5:352653/148‑1 (MQ=255)
                  ccTTCTTCGCCAGGACAACCATAAAAGCGCACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACg                                                                                                                                   >  5:322517/1‑147 (MQ=255)
                   cttcttCGCCAGGACAACCATAAAAGCGCACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGa                                                                                                                                 <  5:90742/148‑1 (MQ=255)
                     tcttcGCCAGGACAACCATAAAAGCGCACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATg                                                                                                                               >  8:133571/1‑148 (MQ=255)
                     tcttaGCGAGGAGAACCATATAAGCGCAACGTGCCGCCTTGCCCATTTTTTTCCAGCCATTTCTTGACGGCTATTGCAGCGTCATAGGCGTCGGTTCCAAGCAAATTGTGTCCGCAACCGTGACCAtt                                                                                                                                                   <  3:117362/128‑1 (MQ=255)
                      cttcGCCAGGACAACCATAAAAGCGCACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGAt                                                                                                                                <  6:169084/146‑1 (MQ=255)
                      cttcGCCAGGACAACCATAAAAGCGCACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGAt                                                                                                                                <  3:138550/146‑1 (MQ=255)
                      cttcGCCAGGACAACCATAAAAGCGCACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATg                                                                                                                               <  8:28539/147‑1 (MQ=255)
                      cttcGCCAGGAAAACCATAAAAGCGCACCTTGCCGCATTGCCCATATTGTTCCAGCCATTTCATGACGGCTATTGCAGCGGCAAAGTAGTCGGTTCCCAGAAAATTGAGTCCGCAACCGTTACCAtat                                                                                                                                                  >  4:117362/1‑126 (MQ=255)
                         cGCCAGGACAACCATAAAAGCGCACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGt                                                                                                                              <  7:60071/145‑1 (MQ=255)
                         cGCCAGGACAACCATAAAAGCGCACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGt                                                                                                                              <  7:244924/145‑1 (MQ=255)
                         cGCCAGGACAACCATAAAAGCGCACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGt                                                                                                                              <  6:204858/145‑1 (MQ=255)
                             aGGACAACCATAAAAGCGCACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTGCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTgc                                                                                                                       >  4:353325/1‑148 (MQ=255)
                              ggACAACCATAAAAGCGCACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTgcg                                                                                                                      <  8:179571/148‑1 (MQ=255)
                              ggACAACCATAAAAGCGCACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTgcg                                                                                                                      <  5:345416/148‑1 (MQ=255)
                                aCAACCATAAAAGCGCACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCAttt                                                                                                                                                  >  3:370786/1‑118 (MQ=255)
                                aCAACCATAAAAGCGCACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCAttt                                                                                                                                                  <  4:370786/118‑1 (MQ=255)
                                aCAACCATAAAAGCGCACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTgcgcg                                                                                                                    >  5:150812/1‑148 (MQ=255)
                                   aCCATAAAAGCGCACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCaa                                                                                                                 <  5:319196/148‑1 (MQ=255)
                                   aCCATAAAAGCGCACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGGCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCaa                                                                                                                 <  6:105488/148‑1 (MQ=255)
                                     cATAAAAGCGCACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGTGCGCAAcc                                                                                                               <  5:191686/148‑1 (MQ=255)
                                      aTAAAAGCGCACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCt                                                                                                              >  1:175565/1‑148 (MQ=255)
                                       tAAAAGCGCACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCtg                                                                                                             <  1:247827/148‑1 (MQ=255)
                                         aaaGCGCACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCtgc                                                                                                            <  3:321436/147‑1 (MQ=255)
                                         aaaGCGCACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCtgc                                                                                                            <  6:322517/147‑1 (MQ=255)
                                         aaaGCGCACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCtgc                                                                                                            <  1:261599/147‑1 (MQ=255)
                                            gcgcACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCtgcttgc                                                                                                        <  3:339111/148‑1 (MQ=255)
                                            gcgcACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCtgcttgc                                                                                                        <  8:390230/148‑1 (MQ=255)
                                              gcACCGTGCCGCCTTGCGCATATTGTTCCAGCCATTTCTAGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCAATTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTg                                                                                                      <  6:262607/148‑1 (MQ=255)
                                              gcACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTAGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCAATTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTg                                                                                                      <  6:257070/148‑1 (MQ=255)
                                                   gTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCtgcttgct                                                                                                       >  2:154357/1‑142 (MQ=255)
                                                   gTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCtgcttgct                                                                                                       <  1:154357/142‑1 (MQ=255)
                                                            tGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGc                                                                                          >  1:229169/1‑146 (MQ=255)
                                                             gCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCa                                                                                        <  4:122148/147‑1 (MQ=255)
                                                             gCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCGGCTTGCTGACTTAAACCTGCCa                                                                                        <  6:350559/147‑1 (MQ=255)
                                                                 atatTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCg                                                                                   >  3:99544/1‑148 (MQ=255)
                                                                  tatTGTTCCAGCCATTTCTTGACTGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGAAACCGTGACCATTTTCACCGTGCGTCACGGa                                                                                                                                 <  2:221583/101‑1 (MQ=255)
                                                                   atTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGt                                                                                  >  6:170931/1‑147 (MQ=255)
                                                                   atTGTTCCAGCCATTTCTTGAAGGCTATTGCAGCGGCAAAGGCTTCGTTTACCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGAt                                                                                                                                >  1:221583/1‑101 (MQ=255)
                                                                   atTGTTCAAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGCCCCTTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTTCTGACTTAAACATGCCAGGGCGt                                                                                  >  2:318299/1‑147 (MQ=255)
                                                                      gTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCat                                                                               <  7:133571/147‑1 (MQ=255)
                                                                                tttCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCcagc                                                                    <  2:175565/148‑1 (MQ=255)
                                                                                 ttCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCTCGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCcag                                                                     >  8:396457/1‑146 (MQ=255)
                                                                                 ttCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTaaa                                                                                               <  6:292820/120‑1 (MQ=255)
                                                                                 ttCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTaaa                                                                                               >  5:292820/1‑120 (MQ=255)
                                                                                 ttCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCcag                                                                     >  1:405759/1‑146 (MQ=255)
                                                                                   cTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCat                                                                               <  6:21739/134‑1 (MQ=255)
                                                                                   cTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCat                                                                               >  5:21739/1‑134 (MQ=255)
                                                                                      gACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAgg                                                                 <  4:329995/145‑1 (MQ=255)
                                                                                      gACGGCTATTGCAGCGGCAAAGGCGGCGGTGCCCAGCAACTTGTGTCCGCAACCGTGACCAGTTTCACCGGGCGGCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAggg                                                                <  6:275776/146‑1 (MQ=255)
                                                                                               tGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGt                                                     >  6:208310/1‑148 (MQ=255)
                                                                                                         aaGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGGCAGGGCGTCat                                                                               <  6:273663/112‑1 (MQ=255)
                                                                                                         aaGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCat                                                                               >  5:273663/1‑112 (MQ=255)
                                                                                                             cggcggTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCaaa                                       >  8:212334/1‑148 (MQ=255)
                                                                                                                cggGTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGa                                    <  6:333052/148‑1 (MQ=255)
                                                                                                                          cAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAg                          >  6:81669/1‑148 (MQ=255)
                                                                                                                           aaaTTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAgg                         >  8:182120/1‑148 (MQ=255)
                                                                                                                             aTTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCa                       <  3:301284/148‑1 (MQ=255)
                                                                                                                              ttGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCAt                      <  3:353325/148‑1 (MQ=255)
                                                                                                                                 tgtCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCAttt                    <  5:170931/147‑1 (MQ=255)
                                                                                                                                  gtCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTgg                  >  8:238552/1‑148 (MQ=255)
                                                                                                                                   tCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTgg                  <  2:405759/147‑1 (MQ=255)
                                                                                                                                    ccGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGa                >  6:157259/1‑148 (MQ=255)
                                                                                                                                      gCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGata              <  5:81669/148‑1 (MQ=255)
                                                                                                                                           cGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACATGCCAGGGCGTCATATTCTCCCAGCAGGGCGACAACCGGTTTGCCTTGACCATACGAAGCAATAAAGGCATTTGGGATATGGcc         >  1:302587/1‑148 (MQ=255)
                                                                                                                                              gACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTACCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTAc      >  1:208494/1‑148 (MQ=255)
                                                                                                                                                  aTTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTg  >  2:67359/1‑148 (MQ=255)
                                                                                                                                                  aTTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTg  >  7:351007/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AACGTTTTACCCGAGCCGCCTTCTTCGCCAGGACAACCATAAAAGCGCACCGTGCCGCCTTGCCCATATTGTTCCAGCCATTTCTTGACGGCTATTGCAGCGGCAAAGGCGGCGGTTCCCAGCAAATTGTGTCCGCAACCGTGACCATTTTCACCGGGCGTCACGGATGTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTG  >  NC_000913/1402795‑1403088

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: