Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,068,317 T→G 6.1% intergenic (+63/+194) ymdF → / ← rutG KGG family protein/pyrimidine permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,068,3170TG 6.1% 347.2 / 12.2 99intergenic (+63/+194)ymdF/rutGKGG family protein/pyrimidine permease
Reads supporting (aligned to +/- strand):  ref base T (38/55);  new base G (2/4);  total (40/59)
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.54e-01

Only 100 of 101 total aligned reads displayed.

TAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGT  >  NC_000913/1068173‑1068464
                                                                                                                                                |                                                                                                                                                   
tAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTAc                                                                                                                                                  >  1:387268/1‑148 (MQ=255)
  aTTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTAccc                                                                                                                                                >  1:166104/1‑148 (MQ=255)
     tCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCggt                                                                                                                                             <  2:223656/148‑1 (MQ=255)
       aaGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCACCGGGTAACCggcggt                                                                                                                                          <  4:571773/149‑1 (MQ=255)
        aGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCggtggt                                                                                                                                          <  2:81084/148‑1 (MQ=255)
          aaTGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCACCGGGTAACCggcggtgat                                                                                                                                       <  4:82236/149‑4 (MQ=255)
            tGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGAt                                                                                                                                      >  2:482890/1‑148 (MQ=255)
            tGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATa                                                                                                                                     <  3:369202/149‑1 (MQ=255)
            tGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATa                                                                                                                                     <  4:219444/149‑1 (MQ=255)
                ccGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACtga                                                                                                                                 <  4:674482/149‑1 (MQ=255)
                ccGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACtg                                                                                                                                  >  2:271921/1‑148 (MQ=255)
                ccGCAGCGCGAATCTGAAGCAGGTAAAAAATGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTTTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACtg                                                                                                                                  <  2:877546/148‑1 (MQ=255)
                 cGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACtga                                                                                                                                 <  2:387268/148‑1 (MQ=255)
                        cgcATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAAt                                                                                                                         <  4:561070/149‑1 (MQ=255)
                         gcATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATg                                                                                                                        >  3:662559/1‑149 (MQ=255)
                            tCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTa                                                                                                                      <  1:271921/148‑1 (MQ=255)
                            tCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTa                                                                                                                      <  2:321523/148‑1 (MQ=255)
                            tCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTa                                                                                                                      >  2:893921/1‑148 (MQ=255)
                             cTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACt                                                                                                                                   >  2:771878/1‑134 (MQ=255)
                             cTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACt                                                                                                                                   <  1:771878/134‑1 (MQ=255)
                              tGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGc                                                                                                                   >  3:113269/1‑149 (MQ=255)
                                aaGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCtt                                                                                                                 <  3:81232/149‑1 (MQ=255)
                                  gCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGc                                                                                                               <  3:848012/149‑1 (MQ=255)
                                   cAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGcc                                                                                                              <  4:438123/149‑1 (MQ=255)
                                   cAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGc                                                                                                               <  2:66697/148‑1 (MQ=255)
                                   cAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGc                                                                                                               >  1:541132/1‑148 (MQ=255)
                                     ggTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACtga                                                                                                                                 >  4:304226/1‑128 (MQ=255)
                                     ggTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACtga                                                                                                                                 <  3:304226/128‑1 (MQ=255)
                                              ggtggtAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGAt                                                                                                   <  4:187634/149‑1 (MQ=255)
                                               gtggtAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCTCGATg                                                                                                  <  4:205303/149‑1 (MQ=255)
                                               gtggtAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGAt                                                                                                   <  1:774102/148‑1 (MQ=255)
                                               gtggtAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATg                                                                                                  <  4:590494/149‑1 (MQ=255)
                                               gtggtAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCACCGGGTAACcggcgg                                                                                                                                           <  4:669292/108‑1 (MQ=255)
                                               gtggtAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCACCGGGTAACCggcgg                                                                                                                                           >  3:669292/1‑105 (MQ=255)
                                                tggtAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATg                                                                                                  <  1:893921/148‑1 (MQ=255)
                                                tggtAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATg                                                                                                  <  1:32146/148‑1 (MQ=255)
                                                  gtAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCtt                                                                                               >  3:288847/1‑149 (MQ=255)
                                                  gtAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCtt                                                                                               >  4:172457/1‑149 (MQ=255)
                                                  gtAAAAGAAGTCACGGCAAAAGCGACAACCAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGAAGCTCGTTAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAAATTGCCTGATGCGCGATGCtt                                                                                               >  4:708130/1‑149 (MQ=255)
                                                   tAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTa                                                                                                                      >  1:201011/1‑125 (MQ=255)
                                                   tAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTa                                                                                                                      <  2:201011/125‑1 (MQ=255)
                                                       agCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAg                                                                                          <  4:33005/149‑1 (MQ=255)
                                                       agCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAg                                                                                          <  4:897314/149‑1 (MQ=255)
                                                        gCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAgg                                                                                         <  3:104233/149‑1 (MQ=255)
                                                        gCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAg                                                                                          >  1:282735/1‑148 (MQ=255)
                                                          aGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATg                                                                                                  >  2:525255/1‑138 (MQ=255)
                                                          aGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATg                                                                                                  <  1:525255/138‑1 (MQ=255)
                                                            tCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCt                                                                                      >  1:271298/1‑148 (MQ=255)
                                                             cACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGc                                                                                                               <  1:576719/122‑1 (MQ=255)
                                                             cACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGc                                                                                                               >  2:576719/1‑122 (MQ=255)
                                                                  cAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCaga                                                                                <  2:824109/148‑1 (MQ=255)
                                                                    aaaGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCagaaga                                                                             <  4:662559/149‑1 (MQ=255)
                                                                      aGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGAt                                                                            <  2:282735/148‑1 (MQ=255)
                                                                         gACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCACCGGGTAACCggcggtgatg                                                                                                                                      >  4:379933/1‑83 (MQ=38)
                                                                         gACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCACCGGGTAACCggcggtgatg                                                                                                                                      <  3:379933/87‑5 (MQ=38)
                                                                          aCAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCaa                                                                       >  3:390732/1‑149 (MQ=255)
                                                                          aCAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCa                                                                        <  2:200788/148‑1 (MQ=255)
                                                                           cAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGa                                                                                                           <  2:148114/112‑1 (MQ=255)
                                                                           cAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGa                                                                                                           >  1:148114/1‑112 (MQ=255)
                                                                                aGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGtt                                                                                                                                         >  3:90226/1‑77 (MQ=255)
                                                                                aGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGtt                                                                                                                                         <  4:90226/77‑1 (MQ=255)
                                                                                aGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAatatat                                                                  <  2:175624/148‑1 (MQ=255)
                                                                                aGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATAtt                                                                 <  4:883592/149‑1 (MQ=255)
                                                                                    gggCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGAt                                                                            <  1:142973/134‑1 (MQ=255)
                                                                                    gggCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGAt                                                                            >  2:142973/1‑134 (MQ=255)
                                                                                    gggCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGaa                                                              <  1:105054/148‑1 (MQ=255)
                                                                                          aTCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCAc                                                       >  3:739949/1‑149 (MQ=255)
                                                                                             aaTGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGt                                                    >  3:3020/1‑149 (MQ=255)
                                                                                              aTGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGt                                                    <  2:541132/148‑1 (MQ=255)
                                                                                                gACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGtttt                                                 <  4:58858/149‑1 (MQ=255)
                                                                                                gACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGttt                                                  <  1:714629/148‑1 (MQ=255)
                                                                                                       gCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGc                                           >  2:768023/1‑148 (MQ=255)
                                                                                                         aTTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGAt                                                                            >  3:244814/1‑113 (MQ=255)
                                                                                                         aTTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGAt                                                                            <  4:244814/113‑1 (MQ=255)
                                                                                                           ttttGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGAt                                      >  4:117976/1‑149 (MQ=255)
                                                                                                             ttGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATa                                     <  2:166104/148‑1 (MQ=255)
                                                                                                                 cTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGAtgcggc                                 >  2:306744/1‑143 (MQ=255)
                                                                                                                  tGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGt                               <  4:390732/149‑1 (MQ=255)
                                                                                                                     ctacagacAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGttt                             <  1:810213/141‑1 (MQ=255)
                                                                                                                      gCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTAc                           <  4:113269/149‑1 (MQ=255)
                                                                                                                        tCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTAcgc                         <  4:288847/149‑1 (MQ=255)
                                                                                                                        tCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTAcgc                         >  3:695694/1‑149 (MQ=255)
                                                                                                                         cGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTAcgc                         >  2:261672/1‑148 (MQ=255)
                                                                                                                         cGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTAcgc                         <  1:482890/148‑1 (MQ=255)
                                                                                                                            cAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCa                                                                                  <  2:896685/88‑1 (MQ=255)
                                                                                                                            cAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCa                                                                                  >  1:896685/1‑88 (MQ=255)
                                                                                                                             aaaaaGCCATCACCGCCGGTTACCCGGTGGTTGATAGTGATGACGAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGATTTGTAGGCCGGATAAGGCGTTTACGCCGCAt                    >  3:762823/1‑149 (MQ=255)
                                                                                                                              aaaaGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCAt                    >  1:171772/1‑148 (MQ=255)
                                                                                                                               aaaGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATcc                  <  4:695694/149‑1 (MQ=255)
                                                                                                                                 aGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCgg                <  3:117976/149‑1 (MQ=255)
                                                                                                                                  gCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc               <  3:172457/149‑1 (MQ=255)
                                                                                                                                    cATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAt             >  3:790347/1‑149 (MQ=255)
                                                                                                                                     aTCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATg            >  3:23966/1‑149 (MQ=255)
                                                                                                                                      tCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCgg                >  2:310034/1‑144 (MQ=255)
                                                                                                                                      tCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCgg                <  1:310034/144‑1 (MQ=255)
                                                                                                                                        aCCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAAc         <  4:23966/149‑1 (MQ=255)
                                                                                                                                         ccgccgGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCAGATAAGGCGTTTACGCCGCATCCGGCATGAACa        >  3:28593/1‑149 (MQ=255)
                                                                                                                                             cgGTCACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGTATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGc    <  4:762823/149‑1 (MQ=255)
                                                                                                                                                ttACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGt  >  2:693822/1‑148 (MQ=255)
                                                                                                                                                ttACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGt  <  2:271298/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
TAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGT  >  NC_000913/1068173‑1068464

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