Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,736,323 A→G 5.1% intergenic (+166/+30) xylR → / ← bax xylose divergent operon transcriptional activator/putative glucosaminidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,736,3230AG 5.1% 1194.4 / 33.5 352intergenic (+166/+30)xylR/baxxylose divergent operon transcriptional activator/putative glucosaminidase
Reads supporting (aligned to +/- strand):  ref base A (175/157);  new base G (10/8);  total (186/166)
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.30e-01

Only 100 of 362 total aligned reads displayed.

GGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTTATAACTCTTCCCTTTGGTCGAGTAGCCCTTCAGCTTGTGAATCATCGCTGTGGCAGTCA  >  NC_000913/3736175‑3736471
                                                                                                                                                    |                                                                                                                                                    
 gCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGAc                                                                                                                                                     <  3:554895/149‑1 (MQ=255)
   gTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACaa                                                                                                                                                   <  1:1019637/149‑1 (MQ=255)
   gTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACaa                                                                                                                                                   <  1:827773/149‑1 (MQ=255)
      acaTCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTc                                                                                                                                                <  2:684018/149‑1 (MQ=255)
        aTCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCgg                                                                                                                                              >  2:539660/1‑149 (MQ=255)
        aTCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCgg                                                                                                                                              >  2:932308/1‑149 (MQ=255)
         tCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGa                                                                                                                                             <  1:247310/149‑1 (MQ=255)
          cAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGAt                                                                                                                                            <  1:539660/149‑1 (MQ=255)
                 cGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTc                                                                                                                                     <  1:594332/149‑1 (MQ=255)
                       gTGGTTGTCGGATGCGGCGTAAACCCCTTATCCGACCTCCCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCAACTGCCTGAGAAAAAGTAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGaaa                                                                                                                              <  2:20602/149‑1 (MQ=255)
                         ggTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAAt                                                                                                                             <  4:1131006/149‑1 (MQ=255)
                            tGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCa                                                                                                                          >  1:613104/1‑149 (MQ=255)
                             gTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAg                                                                                                                                         <  1:851030/133‑1 (MQ=255)
                                      ggcgTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGa                                                                                                                                                      <  2:964085/111‑1 (MQ=255)
                                         gTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATgcgc                                                                                                             <  2:764235/149‑1 (MQ=255)
                                           gAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCaa                                                                                                                                   >  2:864257/1‑125 (MQ=255)
                                            aaCGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCg                                                                                                                                               <  3:531356/112‑1 (MQ=255)
                                             aCGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTTGTCAGTCGGATaa                                                                                                                                          <  4:1005178/116‑1 (MQ=255)
                                             aCGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGAc                                                                                                                                                     <  2:586807/105‑1 (MQ=255)
                                               gCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGAt                                                                                                       >  1:1111087/1‑149 (MQ=255)
                                               gCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGAt                                                                                                       >  3:311762/1‑149 (MQ=255)
                                                ccTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTTGTCAGTCGGATaa                                                                                                                                          >  3:870245/1‑113 (MQ=39)
                                                ccTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTTGTCAGTCGGATaa                                                                                                                                          <  4:870245/113‑1 (MQ=39)
                                                     tCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAAtt                                                                                                                       >  3:227040/1‑127 (MQ=255)
                                                     tCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAAtt                                                                                                                       <  4:227040/127‑1 (MQ=255)
                                                     tCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCg                                                                                                 >  2:455412/1‑149 (MQ=255)
                                                       cGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGtt                                                                                               >  3:725748/1‑149 (MQ=255)
                                                            tACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGaaa                                                                                                                              <  4:204792/113‑1 (MQ=255)
                                                            tACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTa                                                                                          >  3:153643/1‑149 (MQ=255)
                                                              cccGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATc                                                                                        <  2:613104/149‑1 (MQ=255)
                                                               ccGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATcc                                                                                       >  2:813101/1‑149 (MQ=255)
                                                                 gCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTg                                                                                     <  3:851807/149‑1 (MQ=255)
                                                                   cAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGt                                                                                   <  4:1067721/149‑1 (MQ=255)
                                                                    aGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGt                                                                                   >  1:1062763/1‑148 (MQ=255)
                                                                     ggCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTg                                                                                              >  4:650913/1‑136 (MQ=255)
                                                                       cATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACa                                                                                                                    <  4:165724/112‑1 (MQ=255)
                                                                        aTGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACAtt                                                                              >  4:473239/1‑149 (MQ=255)
                                                                           aTAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCa                                                                                                                          <  3:659461/102‑1 (MQ=255)
                                                                              aaaCGCACCAACAACGCTTCAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAAc                                                                        >  2:895033/1‑149 (MQ=255)
                                                                                  gCACCAACAACGCTTCAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGtt                                                                                           <  1:435154/126‑1 (MQ=255)
                                                                                  gCACCAACAACGCTTCAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGt                                                                    >  1:808534/1‑149 (MQ=255)
                                                                                    aCCAACAACGCTTCAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAg                                                                  >  3:489086/1‑149 (MQ=255)
                                                                                       aacaacGCTTCAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAgttg                                                               >  1:1124485/1‑149 (MQ=255)
                                                                                       aacaacGCTTCAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATACGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAgttg                                                               >  4:898244/1‑149 (MQ=255)
                                                                                          aacGCTTCAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTta                                                            >  4:407576/1‑149 (MQ=255)
                                                                                          aacGCTTCAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTta                                                            >  3:120352/1‑149 (MQ=255)
                                                                                             gCTTCAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACACGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTta                                                            >  1:915325/1‑146 (MQ=255)
                                                                                               ttCAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTgg                                                                                             >  3:1042248/1‑111 (MQ=255)
                                                                                                 cAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTGACCGTTGGTGATCCTGGTACATTGCGAACAGGTAGTTGTTATAActct                                                     >  4:89133/1‑149 (MQ=255)
                                                                                                 cAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTTATAActct                                                     <  1:267617/149‑1 (MQ=255)
                                                                                                  aGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTTATAACTCtt                                                    <  3:466597/149‑1 (MQ=255)
                                                                                                    gCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTTATAACTCTTcc                                                  <  4:617388/149‑1 (MQ=255)
                                                                                                      tCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTTATAACTCTTCCCt                                                >  1:465922/1‑149 (MQ=255)
                                                                                                      tCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTTATAACTCTTCCCt                                                >  1:407459/1‑149 (MQ=255)
                                                                                                        gTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTTATAACTCTTCCCttt                                              <  1:728029/149‑1 (MQ=255)
                                                                                                         tAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTTACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTTATAACTCTTCCCTTTg                                             <  2:176294/149‑1 (MQ=255)
                                                                                                          aGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGATACATTGCGAACAGGTAGTTGTTATAACTCTTCCCTTTgg                                            <  4:666720/149‑1 (MQ=255)
                                                                                                           gCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTTATAACTCTTCCCTTTGGt                                           >  4:514719/1‑149 (MQ=255)
                                                                                                              cAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTg                                                                                              <  1:755022/95‑1 (MQ=255)
                                                                                                               aaCTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTTATAACTCTTCCCTTTGGTCGAg                                       <  2:486076/149‑1 (MQ=255)
                                                                                                                 cTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTTATAACTCTTCCCTTTGGTCGAGTa                                     <  2:255472/149‑1 (MQ=255)
                                                                                                                  tGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTTATAACTCTTCCCTTTGGTCGAGTAg                                    <  2:1124485/149‑1 (MQ=255)
                                                                                                                   gCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTTATAACTCTTCCCTTTGGTCGAGTAGc                                   <  4:153643/149‑1 (MQ=255)
                                                                                                                   gCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTTATAACTCTTCCCTTTGGTCGAGTAGc                                   <  4:96066/149‑1 (MQ=255)
                                                                                                                      tGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTTATAACTCTTCCCTTTGGTCGAGTAGCCCt                                <  2:1111087/149‑1 (MQ=255)
                                                                                                                            aaaGTAAAAAGCCTTATCCGACTGACAAGTCGGATaa                                                                                                                                          <  2:231242/37‑1 (MQ=255)
                                                                                                                            aaaGTAAAAAGCCTTATCCGACTGACAAGTCGGATaa                                                                                                                                          >  1:231242/1‑37 (MQ=255)
                                                                                                                            aaaGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAg                                                                      >  4:784147/1‑105 (MQ=255)
                                                                                                                            aaaGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAg                                                                      <  3:784147/105‑1 (MQ=255)
                                                                                                                            aaaGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTTATAACTCTTCCCtt                                               <  1:887817/128‑1 (MQ=255)
                                                                                                                             aaGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTTATAACTCTTCCCTTTGGTCGAGTAGCCCTTCAGc                           >  3:894531/1‑147 (MQ=255)
                                                                                                                              aGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTTATAACTCTTCCCTTTGGTCGAGTAGCCCTTCAGCTtg                        >  2:730258/1‑149 (MQ=255)
                                                                                                                              aGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTTATAACTCTTCCCTTTGGTCGAGTAGCCCTTCAGCTtg                        <  3:18862/149‑1 (MQ=255)
                                                                                                                               gTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTTATAACTCTTcc                                                  >  1:962418/1‑122 (MQ=255)
                                                                                                                               gTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTTATAACTCTTCCCTTTGGTCGAGTAGCCCTTCAGCTtgt                       <  3:701731/149‑1 (MQ=255)
                                                                                                                                   aaaGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTTATAACTCTTCCCTTTg                                             <  4:311823/123‑1 (MQ=255)
                                                                                                                                     aGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTTATAACTCTTCCCTTTGGTCGAGTAGCCCTTCAGCtt                         <  1:12772/141‑1 (MQ=255)
                                                                                                                                      gCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACata                                                                                                                  <  3:930794/51‑1 (MQ=255)
                                                                                                                                       ccTTATCCGACTTGTCAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTTATAACTCTTCCCTTTGGTCGAGTAGCCCTTCAGCTTGTGAatcat                >  2:595719/1‑148 (MQ=255)
                                                                                                                                       ccTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTTATAACTCTTCCCTTTGGTCGAGTAGCCCTTCAGCTTGTGAatcatc               >  4:454617/1‑149 (MQ=255)
                                                                                                                                        cTTATCCGACTTGTCAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTTATAACTCTTCCCTTTGGTCGAGTAGCCCTTCAGCTTGTGAATCATCg              >  2:843672/1‑149 (MQ=255)
                                                                                                                                         ttATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTTATAACTCTTCCCtt                                               >  3:409983/1‑115 (MQ=255)
                                                                                                                                         ttATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTTATAACTCTTCCCTTTGGTCGa                                        <  1:202020/122‑1 (MQ=255)
                                                                                                                                         ttATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTTATAACTCTTCCCTTTGGTCGAGTAGCCCTTCAGCTTGTGAatcatc               >  4:843551/1‑147 (MQ=255)
                                                                                                                                         ttATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTTATAACTCTTCCCTTTGGTCGAGTAGCCCTTCAGCTTGTGAATCATCGc             >  4:121225/1‑149 (MQ=255)
                                                                                                                                          tATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAgttg                                                               <  4:601821/98‑1 (MQ=255)
                                                                                                                                          tATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTTATAACTCTTCCCTTTGGTCGAGTAGCCCTTCAGCTTGTGAATCATCGCt            >  2:655357/1‑149 (MQ=255)
                                                                                                                                           aTCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTa                                                                                                     <  3:388077/59‑1 (MQ=255)
                                                                                                                                           aTCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTa                                                                                                     >  4:388077/1‑59 (MQ=255)
                                                                                                                                            tCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTTATAACTCTTCCCTTTGGTCGAGTAGCCCTTCAGCTTGTGAATCATCGCtgt          >  1:583434/1‑149 (MQ=255)
                                                                                                                                               gACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTTATaa                                                         >  2:347021/1‑99 (MQ=255)
                                                                                                                                               gACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTTATaa                                                         <  1:347021/99‑1 (MQ=255)
                                                                                                                                               gACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTTATAACTCTTCCCTTTGGTCGAGTAGCCCTTCAGCTTGTGAatca                 <  3:606690/139‑1 (MQ=255)
                                                                                                                                               gACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTTATAACTCTTCCCTTTGGTCGAGTAGCCCTTCAGCTTGTGAATCATCGCTGTgg        >  4:72727/1‑148 (MQ=255)
                                                                                                                                               gACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTTATAACTCTTCCCTTTGGTCGAGTAGCCCTTCAGCTTGTGAATCATCGCTGTGGc       >  1:508373/1‑149 (MQ=255)
                                                                                                                                               gACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTTATAACTCTTCCCTTTGGTCGAGAAGCCCTTCAGCTTGTGAATCATTGCTGTGGc       >  3:108792/1‑149 (MQ=255)
                                                                                                                                                 cTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTTATAACTCTTCCCTTTGGTCGAGTAGCCCTTCAGCTTGTGAATCATCGCTGTGGCAg     <  2:407459/149‑1 (MQ=255)
                                                                                                                                                   gACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTTATAACTCTTCCCTTTGGTCGAGTAGCCCTTCAGCTTGTGAATCATCGCTGTGGCAGTc   >  2:1037960/1‑149 (MQ=255)
                                                                                                                                                    aCAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTTATAACTCTTCCCTTTGGTCGAGTAGCCCTTCAGCTTGTGAATCATCGCTGTGGCAGTCa  >  2:1003052/1‑149 (MQ=255)
                                                                                                                                                    aCAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTTATAACTCTTCCCCTTGGTc                                          >  3:629037/1‑109 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTTATAACTCTTCCCTTTGGTCGAGTAGCCCTTCAGCTTGTGAATCATCGCTGTGGCAGTCA  >  NC_000913/3736175‑3736471

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