Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,962,966 T→A 6.0% Q64L (CAG→CTG)  ygdB ← DUF2509 domain‑containing protein YgdB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,962,9660TA 6.0% 268.5 / 5.3 83Q64L (CAG→CTG) ygdBDUF2509 domain‑containing protein YgdB
Reads supporting (aligned to +/- strand):  ref base T (46/32);  new base A (3/2);  total (49/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

GTTTCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTACCCAAAACCAGCAGCATC  >  NC_000913/2962821‑2963115
                                                                                                                                                 |                                                                                                                                                     
gTTTCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCca                                                                                                                                                   >  1:1137958/1‑150 (MQ=255)
  ttCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCcagc                                                                                                                                                 <  1:1789123/150‑1 (MQ=255)
  ttCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGCTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCcagc                                                                                                                                                 <  2:1739125/150‑1 (MQ=255)
   tCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCcagca                                                                                                                                                <  2:1729757/150‑1 (MQ=255)
    ccATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGTCTTCGGTTTCAGCGTACTGGGAGCACTGAACTGCTGGCTGCCTCTGGcagcag                                                                                                                                               >  2:1563466/1‑150 (MQ=255)
     cATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGt                                                                                                                                              >  2:378465/1‑150 (MQ=255)
         aaTGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCcagc                                                                                                                                                 <  1:662148/143‑1 (MQ=255)
         aaTGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCcagc                                                                                                                                                 >  2:662148/1‑143 (MQ=255)
           tGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCAGACGCAGCCAGCAGTTCagtg                                                                                                                                        <  2:1373098/150‑4 (MQ=255)
           tGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCAGACGCAGCCAGCAGTTCagtg                                                                                                                                        <  2:365784/150‑4 (MQ=255)
                tCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTccc                                                                                                                                   <  2:1423747/150‑1 (MQ=255)
                   ccTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCAc                                                                                                                                >  2:1447643/1‑150 (MQ=255)
                      gTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGcc                                                                                                                             >  1:492114/1‑150 (MQ=255)
                         cGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAgc                                                                                                                          >  1:657710/1‑150 (MQ=255)
                           ccACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTAACTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTCCATTTTACCCCACGCCAgcgc                                                                                                                        >  2:725695/1‑150 (MQ=255)
                            cacaACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGcc                                                                                                                       >  2:234084/1‑150 (MQ=255)
                             acaACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCg                                                                                                                      <  1:1540165/150‑1 (MQ=255)
                              caACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGa                                                                                                                     >  1:229739/1‑150 (MQ=255)
                                    aaCGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAAc                                                                                                               >  2:582085/1‑150 (MQ=255)
                                       gCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGAt                                                                                                            >  2:1877526/1‑150 (MQ=255)
                                         cTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATgg                                                                                                          <  2:229739/150‑1 (MQ=255)
                                         cTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATgg                                                                                                          >  1:448234/1‑150 (MQ=255)
                                           tCATAACCGGCAAGCAATACGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCAGCGGTTTCAGCGTACTGCGAGCACGGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGcc                                                                                                        <  2:1842899/150‑1 (MQ=255)
                                             aTAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTg                                                                                                      >  2:1112517/1‑150 (MQ=255)
                                                aCCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCAACGGTTTCAGCGTACTGCGAGAACTGAACTGCTGGCTGCGTCTGCAAGCAGTCCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTggcg                                                                                                   >  2:1427344/1‑150 (MQ=255)
                                                   ggCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGgcggcg                                                                                                >  1:951951/1‑150 (MQ=255)
                                                   ggCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGgcggcg                                                                                                <  1:234084/150‑1 (MQ=255)
                                                   ggCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGgcggcg                                                                                                <  1:1112517/150‑1 (MQ=255)
                                                    gCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGgcgc                                                                                               <  1:378465/150‑1 (MQ=255)
                                                      aatcaatAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCaa                                                                                             >  2:1985993/1‑150 (MQ=255)
                                                          aatAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAAtggc                                                                                         >  1:344326/1‑148 (MQ=255)
                                                          aatAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAAtgac                                                                                         <  1:475437/150‑1 (MQ=255)
                                                          aatAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAAtgac                                                                                         <  1:306012/150‑1 (MQ=255)
                                                             aaGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAAtgactga                                                                                      <  1:1877526/150‑1 (MQ=255)
                                                               ggCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATgactgact                                                                                    <  2:1958468/150‑1 (MQ=255)
                                                               ggCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATgactgact                                                                                    >  1:1654987/1‑150 (MQ=255)
                                                                  ttCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAg                                                                                                                           <  2:1515087/108‑1 (MQ=255)
                                                                  ttCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAg                                                                                                                           >  1:1515087/1‑108 (MQ=255)
                                                                    cATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCa                                                                               <  1:1328494/150‑1 (MQ=255)
                                                                        aTCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCt                                                                           >  2:212708/1‑150 (MQ=255)
                                                                         tCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTc                                                                          >  2:14302/1‑150 (MQ=255)
                                                                          cTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCa                                                                         >  2:976888/1‑150 (MQ=255)
                                                                          cTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCa                                                                         >  2:385684/1‑150 (MQ=255)
                                                                          cTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCa                                                                         >  2:374515/1‑150 (MQ=255)
                                                                           tGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTcccc                                                                                                                                  <  2:306644/92‑1 (MQ=255)
                                                                           tGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTcccc                                                                                                                                  >  1:306644/1‑92 (MQ=255)
                                                                           tGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCAc                                                                        >  1:1777113/1‑150 (MQ=255)
                                                                           tGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCAGACGCAGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCAc                                                                        >  1:1552142/1‑150 (MQ=255)
                                                                             gcAGTAAACGCAATCAAACGTGGGGATCGGTTTCAGCTTAGTGGGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATTGCGTGGCGGCGCAATGAGTGACTTTCCATGCTCAcgc                                                                      <  2:1203381/149‑1 (MQ=255)
                                                                              cAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCAcgcg                                                                     <  1:1447643/150‑1 (MQ=255)
                                                                                       caaacaaaCCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCaaa                                                            >  1:1698368/1‑150 (MQ=255)
                                                                                        aaacaaacCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCaaaa                                                           >  1:862189/1‑150 (MQ=255)
                                                                                           caaacCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTg                                                        >  1:255947/1‑150 (MQ=255)
                                                                                            aaacCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGc                                                       <  1:582085/150‑1 (MQ=255)
                                                                                             aacCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCg                                                      >  2:428733/1‑150 (MQ=255)
                                                                                               ccTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGAt                                                    >  1:1724422/1‑150 (MQ=255)
                                                                                                cTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTACCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATc                                                   <  1:725695/150‑1 (MQ=255)
                                                                                                 tGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATcc                                                  <  1:29588/150‑1 (MQ=255)
                                                                                                     cATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGCCTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGtt                                              >  2:1758121/1‑150 (MQ=255)
                                                                                                     cATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGtt                                              <  1:1985993/150‑1 (MQ=255)
                                                                                                      aTCTGTTTCAGCGTACTGCTAGCACTGAACTGCTGGCTTCGTCTGCCAGCAGTGCATTTTTCCCCACGCAAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTg                                             <  2:1854645/150‑1 (MQ=255)
                                                                                                      aTCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTg                                             >  2:489359/1‑150 (MQ=255)
                                                                                                       tCGGTTGCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTgg                                            <  1:878177/150‑1 (MQ=255)
                                                                                                        cGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCTCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGc                                           <  1:457288/150‑1 (MQ=255)
                                                                                                         ggTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAAtg                                                                                           >  2:1870590/1‑101 (MQ=255)
                                                                                                         ggTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAAtg                                                                                           <  1:1870590/101‑1 (MQ=255)
                                                                                                         ggTTTCAGCGTACAGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCt                                          <  2:650007/150‑1 (MQ=255)
                                                                                                            ttCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCAt                                       >  2:598937/1‑150 (MQ=255)
                                                                                                                gCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTcc                                    >  2:219293/1‑149 (MQ=255)
                                                                                                                gCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTcc                                    <  1:219293/149‑1 (MQ=255)
                                                                                                                   tACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGt                                >  1:1474082/1‑150 (MQ=255)
                                                                                                                    aCTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTa                               <  2:1529549/150‑1 (MQ=255)
                                                                                                                             aCTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAAc                      >  1:1723911/1‑150 (MQ=255)
                                                                                                                             aCTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAAc                      >  2:94207/1‑150 (MQ=255)
                                                                                                                             aCTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAAc                      <  2:190820/150‑1 (MQ=255)
                                                                                                                              cTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACt                     >  1:729848/1‑150 (MQ=255)
                                                                                                                              cTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACt                     >  1:736736/1‑150 (MQ=255)
                                                                                                                              atgcaCTGCTGGCAGACGCAGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACt                     >  1:724817/5‑150 (MQ=255)
                                                                                                                                gAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTAc                   >  1:717517/1‑150 (MQ=255)
                                                                                                                                 caCTGCTGGCAGACGCAGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTAcc                  >  2:1076557/2‑150 (MQ=255)
                                                                                                                                 aaCTGCTGGTTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTAcc                  >  1:595906/1‑150 (MQ=255)
                                                                                                                                     gctggctgCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTACCCaaa              >  1:759594/1‑150 (MQ=255)
                                                                                                                                               tCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTACCCAAAACCAgcagca    >  1:644492/1‑150 (MQ=255)
                                                                                                                                                 tGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTACCCAAAACCAGCAGCATc  >  2:919121/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GTTTCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTACCCAAAACCAGCAGCATC  >  NC_000913/2962821‑2963115

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