Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,689,821 T→G 5.5% intergenic (+70/‑31) manA → / → ydgA mannose‑6‑phosphate isomerase/DUF945 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,689,8210TG 5.5% 271.1 / 8.3 92intergenic (+70/‑31)manA/ydgAmannose‑6‑phosphate isomerase/DUF945 family protein
Reads supporting (aligned to +/- strand):  ref base T (45/42);  new base G (3/2);  total (48/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.87e-01

AGCGTTTATTGCCGCCAACGAATCACCGGTGACTGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGG  >  NC_000913/1689673‑1689969
                                                                                                                                                    |                                                                                                                                                    
aGCGTTTATTGCCGCCAACGAATCACCGGTGACTGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGt                                                                                                                                                      >  5:194765/1‑149 (MQ=255)
  cGTTTATTGCCGCCAACGAATCACCGGTGACTGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTcc                                                                                                                                                    >  4:244976/1‑149 (MQ=255)
   gTTTATTGCCGCCAACGAATCACCGGTGACTGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGGTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCg                                                                                                                                                   >  7:176871/1‑149 (MQ=255)
           ccgccAACGAATCACCGGTGACTGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGAttt                                                                                                                                           >  7:251233/1‑149 (MQ=255)
               cAACGAATCACCGGTGACTGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCAtt                                                                                                                                       >  5:201331/1‑149 (MQ=255)
                    aaTCACCGGTGACTGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATg                                                                                                                                  >  1:96631/1‑149 (MQ=255)
                        aCCGGTGACTGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGAt                                                                                                                               >  2:112173/1‑148 (MQ=255)
                         ccGGTGACTGTCAAAGCCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGACTGATTTCATTTTATGGATaa                                                                                                                             <  5:130937/149‑1 (MQ=255)
                           ggTGACTGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTTAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATc                                                                                                                           <  2:96631/149‑1 (MQ=255)
                            gTGACTGTCAAAGGCCACGGCCGTTTAGCGCGTGTGTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCa                                                                                                                          <  7:195331/149‑1 (MQ=255)
                               aCTGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTa                                                                                                                       <  3:111927/149‑1 (MQ=255)
                                cTGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTAt                                                                                                                      <  6:347301/149‑1 (MQ=255)
                                cTGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTAt                                                                                                                      <  3:137184/149‑1 (MQ=255)
                                 tGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCAatat                                                                                                                      >  6:229713/1‑148 (MQ=255)
                                 tGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTACGCCTGATTTCATTTTATGGATTATCATTAt                                                                                                                      >  3:192590/1‑148 (MQ=255)
                                 tGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCGGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTAGTAAAGGCTTATAACACCTTCAGGCGGCCAGTCAGCCTGATTTCATTTTATGGATAATCATTAt                                                                                                                      >  1:385166/1‑148 (MQ=255)
                                  gTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATg                                                                                                                     >  5:14510/1‑148 (MQ=255)
                                  gTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGACTGGCCGCCTGa                                                                                                                                              <  1:249128/123‑1 (MQ=255)
                                  gTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGACTGGCCGCCTGa                                                                                                                                              >  2:249128/1‑123 (MQ=255)
                                   tCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAACTTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGaa                                                                                                                   <  8:176871/149‑1 (MQ=255)
                                    cAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGACTGGCCGCCTGaag                                                                                                                                            <  1:106616/123‑3 (MQ=255)
                                    cAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGACTGGCCGCCTGaag                                                                                                                                            >  2:106616/1‑121 (MQ=255)
                                     aaaGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTAt                                                                                                                      >  2:270576/1‑144 (MQ=255)
                                     aaaGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTAt                                                                                                                      <  1:270576/144‑1 (MQ=255)
                                     aaaGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAAt                                                                                                                  >  2:236599/1‑148 (MQ=255)
                                       aGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATaaa                                                                                                               <  8:251233/149‑1 (MQ=255)
                                       aGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATaaa                                                                                                               >  8:79980/1‑149 (MQ=255)
                                          ccACGGGCGTTTAGCGCCGGTTTACAACAATCTGTAAGATCCTTCTGAAAAAATTAACATCTCTTGCTAACCTTATTAAAGGCTTAAAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCg                                                                                                            <  2:265223/149‑1 (MQ=255)
                                           cACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACCTCTCGTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCAGCCTGAGGTCAGTGTATGGATAATCAGTATGAAGAAATCGc                                                                                                           >  6:270150/1‑149 (MQ=255)
                                             cGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGccg                                                                                                         <  6:14510/149‑3 (MQ=255)
                                             cGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTg                                                                                                         <  7:79980/149‑1 (MQ=255)
                                              ggCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTgg                                                                                                        <  1:112173/149‑1 (MQ=255)
                                                       gcgcGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATCGATAATCATTATGAATAAATCGCTGGTAGCGGTa                                                                                                <  4:129314/148‑1 (MQ=255)
                                                         gcgTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTa                                                                                                      >  1:9256/1‑140 (MQ=255)
                                                         gcgTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTa                                                                                                      <  2:9256/140‑1 (MQ=255)
                                                              tttACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCGGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCAt                                                                                        >  8:30237/1‑149 (MQ=255)
                                                               ttACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCAtt                                                                                       >  2:214916/1‑149 (MQ=255)
                                                                 acaacaAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCAttgt                                                                                     >  5:24848/1‑149 (MQ=255)
                                                                 acaacaAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCAttgt                                                                                     >  4:217799/1‑149 (MQ=255)
                                                                  caacaaGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCAttgtt                                                                                    <  3:244976/149‑1 (MQ=255)
                                                                   aacaaGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCgg                                                                                                  <  5:2897/134‑1 (MQ=255)
                                                                   aacaaGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCgg                                                                                                  >  6:2897/1‑134 (MQ=255)
                                                                     caaGCTGTAAGAGCTTACTGAAAAAATTAACATCTTTTGATAATCTTATTAAAGGCTTATAAAACCTTCAGGCGGCCAGTCCTCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCTGTAGGCGTCATTGTTgcg                                                                                 <  8:377458/149‑1 (MQ=255)
                                                                     caaGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATaa                                                                                                                             >  1:93365/1‑105 (MQ=255)
                                                                     caaGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATaa                                                                                                                             <  2:93365/105‑1 (MQ=255)
                                                                     caaGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTAt                                                                                                                      <  6:118241/112‑1 (MQ=255)
                                                                     caaGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTAt                                                                                                                      >  5:118241/1‑112 (MQ=255)
                                                                     caaGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTgcg                                                                                 >  7:411925/1‑149 (MQ=255)
                                                                         cTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAg                                                                             >  3:279340/1‑149 (MQ=255)
                                                                            tAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCg                                                                          >  7:343236/1‑149 (MQ=255)
                                                                            tAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCg                                                                          >  2:395760/1‑149 (MQ=255)
                                                                                  cTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCt                                                                    <  6:194765/149‑1 (MQ=255)
                                                                                     aCTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTg                                                                   <  6:201331/147‑1 (MQ=255)
                                                                                     aCTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGa                                                                 >  7:28695/1‑149 (MQ=255)
                                                                                      cTGAATAAATTAACATCTCTTGCTAAGCTTATTAAATGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGATGGTAGCGGTAGGCGTCATTGTTCCGCTAGGCGTAGTCTgta                                                                 <  6:307818/148‑3 (MQ=255)
                                                                                       tGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATAGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGa                                                                 >  2:293110/1‑147 (MQ=255)
                                                                                        gAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGa                                                                 >  2:200472/1‑146 (MQ=255)
                                                                                        gAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGa                                                                 <  1:200472/146‑1 (MQ=255)
                                                                                         aaaaaaTTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTg                                                                                                         <  3:117272/105‑1 (MQ=255)
                                                                                         aaaaaaTTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTg                                                                                                         >  4:117272/1‑105 (MQ=255)
                                                                                          aaaaaTTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCGGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAggc                                                            <  7:30237/149‑1 (MQ=255)
                                                                                          aaaaaTTAACATCTCTTGCTAAGCTTATTAAAGCCTAATAACACCTTCAGGCGGCCAGTCCGCCTGATATCATTTTATGGATAATCATTATTAATAAATCGCTGGAAGAGGTAGCCTTCAt                                                                                        >  2:396677/1‑121 (MQ=255)
                                                                                             aaTTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAggcggc                                                         <  1:236599/149‑1 (MQ=255)
                                                                                                   cATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATgg                                                   <  6:24848/149‑1 (MQ=255)
                                                                                                      ctctTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGta                                                 >  6:253601/1‑148 (MQ=255)
                                                                                                      ctctTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGta                                                 >  6:360801/1‑148 (MQ=255)
                                                                                                       tctTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGtata                                               >  3:96326/1‑149 (MQ=255)
                                                                                                        catGCTAAGCTTATTAAAGGCATATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATAc                                              <  1:325298/147‑1 (MQ=255)
                                                                                                         ttGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGta                                                 <  4:279340/145‑1 (MQ=255)
                                                                                                            cTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGc                                          <  1:395760/149‑1 (MQ=255)
                                                                                                             tAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGc                                          >  5:252237/1‑148 (MQ=255)
                                                                                                             tAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGc                                          >  8:167485/1‑148 (MQ=255)
                                                                                                                 cTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCaag                                       >  5:114746/1‑147 (MQ=255)
                                                                                                                  ttattaAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATg                                                    <  7:309230/133‑1 (MQ=255)
                                                                                                                  ttattaAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATg                                                    >  8:309230/1‑133 (MQ=255)
                                                                                                                    attaAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATgg                                                   >  6:235252/1‑132 (MQ=255)
                                                                                                                     ttaAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCAt                                                     >  3:132927/1‑129 (MQ=255)
                                                                                                                     ttaAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCAt                                                     <  4:132927/129‑1 (MQ=255)
                                                                                                                     ttaAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGAtt                                 <  4:96326/149‑1 (MQ=255)
                                                                                                                     ttaAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGAt                                  <  5:360801/148‑1 (MQ=255)
                                                                                                                     ttaAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGAt                                  <  1:214916/148‑1 (MQ=255)
                                                                                                                       aaaGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGa                                                                 <  6:333136/115‑1 (MQ=255)
                                                                                                                       aaaGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTg                                >  5:311722/1‑148 (MQ=255)
                                                                                                                        aaGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGa                               >  4:117264/1‑148 (MQ=255)
                                                                                                                         aGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGaaa                             <  8:343236/149‑1 (MQ=255)
                                                                                                                               ataACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCAtc                       >  8:333103/1‑149 (MQ=255)
                                                                                                                                  atgaaatCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCAtctc                     >  7:38930/7‑148 (MQ=255)
                                                                                                                                  acacCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCg                    <  8:28695/149‑1 (MQ=255)
                                                                                                                                   cacCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCg                    >  4:157685/1‑148 (MQ=255)
                                                                                                                                    acCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCg                    >  3:275402/1‑147 (MQ=255)
                                                                                                                                          aGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATg            <  7:333103/149‑1 (MQ=255)
                                                                                                                                           ggcggcCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATgg           <  4:275402/149‑1 (MQ=255)
                                                                                                                                                  aGTCCTCCGGATTTCATTTTATGGTTAATCATTATGAATAACTCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCa    <  5:229713/149‑1 (MQ=255)
                                                                                                                                                    tCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAgg  >  1:195931/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGCGTTTATTGCCGCCAACGAATCACCGGTGACTGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGG  >  NC_000913/1689673‑1689969

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