Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 5,171 C→T intergenic (+151/‑63) thrC → / → yaaX threonine synthase/DUF2502 domain‑containing protein YaaX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009135,1710CT98.1% 173.4 / NA 53intergenic (+151/‑63)thrC/yaaXthreonine synthase/DUF2502 domain‑containing protein YaaX
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (27/25);  total (28/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAA  >  NC_000913/4928‑5406
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                            
gagCTGGCAGAACGTGTTGATTTTCCCTTTTTTTCTCAATATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTATATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAAGGCGGTAACTTAGAGAGTAGGAGTGCGGTGAATAACAACCGTCGTTCTCa                                                                                                                                                                                                                                      <  2:315340/251‑1 (MQ=255)
             gTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTcc                                                                                                                                                                                                                         >  2:525338/1‑251 (MQ=255)
                   aTTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGatat                                                                                                                                                                                                                   >  2:473191/1‑251 (MQ=255)
                       aCCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGAc                                                                                                                                                                                                               <  1:409543/251‑1 (MQ=255)
                         ccTTGTTTCCACAAAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCAGCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGTAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGAccc                                                                                                                                                                                                             <  2:623051/251‑1 (MQ=255)
                            tgtttttaaaaagtgtggccgccgATTTTGCTGCGTTGCGTAAATTGATGATGAAGCATCAGTAAAATCTATTCATTATCACAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAGCGACAGGGAAAAAGGAGAAAGTCTCAATAAATGCGGTAACTTAGGGAGTTGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCAt                                                                                                                                                                                                           <  2:236404/232‑1 (MQ=255)
                                      aaaTCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATGAGGAGTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATGACGGGCAATg                                                                                                                                                                                                <  2:595073/250‑1 (MQ=255)
                                        aTCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGAt                                                                                                                                                                                              >  1:537329/1‑251 (MQ=255)
                                            gCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaaa                                                                                                                                                                                          <  1:558458/251‑1 (MQ=255)
                                              ccgccgATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTAATCATTATCTCAATCAGGCCGGGTGTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAgg                                                                                                                                                                                        <  1:229716/251‑1 (MQ=255)
                                                 ccgATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAACCGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGt                                                                                                                                                                                     <  1:473191/251‑1 (MQ=255)
                                                   gATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTa                                                                                                                                                                                    >  1:581382/1‑250 (MQ=255)
                                                      tttGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑t                                                                                                                                                                               >  1:487551/1‑251 (MQ=255)
                                                        tgctgcGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑tat                                                                                                                                                                             >  1:329579/1‑249 (MQ=255)
                                                                 gcgtAAATTGAGGAGGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAATAACGACAGGGAAAAAGGAGAAAGTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTGCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGAt                                                                                                                                                                    <  1:36974/251‑1 (MQ=255)
                                                                   gtAAATTGATGATGAATAATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGc                                                                                                                                                                  >  2:87016/1‑251 (MQ=255)
                                                                     aaaTTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCa                                                                                                                                                                 >  1:346485/1‑250 (MQ=255)
                                                                          gatgatgaATCATCAGTAAAAACTATTCATTATCTCAATCAGGCCGGGTTTGCTTGTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCAACGCGTAATCTCCGGATAACGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATcga                                                                                                                                                            <  2:525317/250‑3 (MQ=255)
                                                                             gatgaATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCg                                                                                                                                                         <  1:400548/250‑1 (MQ=255)
                                                                                gaATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTAc                                                                                                                                                      <  1:164638/250‑1 (MQ=255)
                                                                                    catcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCa                                                                                                                                                 <  2:166611/251‑1 (MQ=255)
                                                                                       caGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACtt                                                                                                                                              <  2:2479/251‑1 (MQ=255)
                                                                                             actCTAGTCATTATCTCACTCAGGCCGGGTTTACTTATATGCAGCC‑GGCTTTTTTTTCAAGAAATTATGGAGAAAAACGACAGGGAAATAGGTGAAATTCTCAATAAATGCGGTAACTAAGAGATTAGGATTGCGGAGAATAACAACCGCCGTACTCATCGCGTAATCTCCGGATATAAACCCATAACGGGCAATGCTAACAGGA‑TCACCATATGAAAAAGCTGAAATCTATCGTACTCGAACTTTCCCTg                                                                                                                                        >  2:109328/3‑251 (MQ=255)
                                                                                              aTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTgg                                                                                                                                       >  2:178695/1‑251 (MQ=255)
                                                                                                    tCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTgg                                                                                                                                 >  1:411235/1‑251 (MQ=255)
                                                                                                      aTTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTGTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTc                                                                                                                               <  2:212253/251‑1 (MQ=255)
                                                                                                                 cAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgc                                                                                                                    >  2:358045/1‑251 (MQ=255)
                                                                                                                     ccGGGTGTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGATAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGAGTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAATATGCAATCTATCGTACTCGCACTTTCCCTCGTTCTGGCCGCTCCCATGgcagca                                                                                                                <  2:207439/251‑1 (MQ=255)
                                                                                                                               ttttATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCgg                                                                                                      >  1:448274/1‑251 (MQ=255)
                                                                                                                                tttATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCgg                                                                                                      >  1:329905/1‑250 (MQ=255)
                                                                                                                                             gCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGt                                                                                         <  2:305462/251‑1 (MQ=255)
                                                                                                                                             gCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGACAGCACAGGCTGCGGAAATTACGTTAGt                                                                                         >  1:329642/1‑251 (MQ=255)
                                                                                                                                                     aTGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTc                                                                                   >  1:143343/1‑249 (MQ=255)
                                                                                                                                                            aaTTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAAtt                                                                          <  1:178695/251‑1 (MQ=255)
                                                                                                                                                                 tGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTAc                                                                        >  1:550629/1‑248 (MQ=255)
                                                                                                                                                                      aaaaaCGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGc                                                                <  2:512233/251‑1 (MQ=255)
                                                                                                                                                                      aaaaaCGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGc                                                                <  1:59013/251‑1 (MQ=255)
                                                                                                                                                                          acgaccgggAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGCTATCGACCCATAACGGGCCATGATAAAAGGAGTAACC‑TATGCAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCAATGGCAGCACGGGCTGCGGAAATTAAGTTAGTCCCGTCAGTAAAATTAAAGATAGGCGATc                                                            >  1:312192/7‑251 (MQ=255)
                                                                                                                                                                             aCAGGGAAAAAGGAGACATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTg                                                         >  1:130876/1‑251 (MQ=255)
                                                                                                                                                                                   aaaaaGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATc                                                   >  2:179172/1‑251 (MQ=255)
                                                                                                                                                                                     aaaGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGAGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATACCGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGt                                                 >  1:89584/1‑251 (MQ=255)
                                                                                                                                                                                      aaGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTg                                                <  1:132884/251‑1 (MQ=255)
                                                                                                                                                                                             aaTTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTAtt                                          >  2:181712/1‑250 (MQ=255)
                                                                                                                                                                                              aTTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTa                                         >  1:219553/1‑250 (MQ=255)
                                                                                                                                                                                                   cAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TAGGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTg                                      <  1:363592/248‑1 (MQ=255)
                                                                                                                                                                                                          tGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGt                            >  1:583388/1‑251 (MQ=255)
                                                                                                                                                                                                                  cTTAGAGCTTAGGACTGCGGAGAATAACAACAGTCGTTCTCATAGCGTAATCTCCGGATATCGACCAATAACGGGCACTGATAAAAAGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCAATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGACAAATTACAGCACGGCGAACGTGATAATCGTGGCTCTTACTGGGATGGAGGCCACTGgcg                    >  2:275291/1‑251 (MQ=255)
                                                                                                                                                                                                                         aTTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATACTCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCAc             >  2:499542/1‑251 (MQ=255)
                                                                                                                                                                                                                             ggATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCTGATATCGACCCATAACGGGCAATGATAAAATGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCTCACTTTCCCTTGTTCTTGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCTTCAGTAAAATTACAGATAGGCGATCGTg                                                         <  2:60704/203‑1 (MQ=255)
                                                                                                                                                                                                                             ggATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGAACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTg                                                         >  1:60704/1‑203 (MQ=255)
                                                                                                                                                                                                                                  gCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACGGGGATGGAGGTCACTGGCGCGACCACGGCtggtg     <  1:165200/250‑1 (MQ=255)
                                                                                                                                                                                                                                   cGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACACATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGACTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGAGCGTGATAATCGTGGCTATTATTGGGTTGGAGGTCACTGGCGCGACCACGACtggtgg    >  1:478228/1‑250 (MQ=255)
                                                                                                                                                                                                                                    ggAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGaa  <  1:401410/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                            
GAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAA  >  NC_000913/4928‑5406

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