Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,597,125 A→G 5.3% intergenic (+113/+25) lgoD → / ← opgB L‑galactonate oxidoreductase; L‑gulonate oxidoreductase/OPG periplasmic biosynthetic phosphoglycerol transferases I (membrane‑bound) and II (soluble)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,597,1250AG 5.3% 302.9 / 5.0 94intergenic (+113/+25)lgoD/opgBL‑galactonate oxidoreductase; L‑gulonate oxidoreductase/OPG periplasmic biosynthetic phosphoglycerol transferases I (membrane‑bound) and II (soluble)
Reads supporting (aligned to +/- strand):  ref base A (45/44);  new base G (3/2);  total (48/46)
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.08e-01

TGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCCGGCGGCACAATTACCAGTGTGTC  >  NC_000913/4596982‑4597272
                                                                                                                                               |                                                                                                                                                   
tGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGa                                                                                                                                                >  2:178403/1‑149 (MQ=255)
  aGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGCTTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGGTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGAcc                                                                                                                                              <  1:70918/149‑1 (MQ=255)
    ttaattaaAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGc                                                                                                                                                  <  8:271877/143‑1 (MQ=255)
    ttaattaaAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGc                                                                                                                                                  >  7:271877/1‑143 (MQ=255)
    ttaattaaAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTa                                                                                                                                            >  7:156232/1‑149 (MQ=255)
      aattaaAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTac                                                                                                                                           <  6:168940/148‑1 (MQ=255)
      aattaaAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTac                                                                                                                                           <  5:255756/148‑1 (MQ=255)
      aattaaAGGCGTAATTACATTCTGATAAGGCGAGATTATTAAAGTAGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGGTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCATAGGGTGCGACCTac                                                                                                                                           <  1:114700/148‑1 (MQ=255)
       attaaAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGCCCTaca                                                                                                                                          >  3:278024/1‑148 (MQ=255)
        ttaaAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTAcaca                                                                                                                                        >  6:248154/1‑149 (MQ=255)
               cGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACAtcagtca                                                                                                                                 <  8:57480/149‑1 (MQ=255)
               cGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACAtcagtca                                                                                                                                 <  4:15608/149‑1 (MQ=255)
                    ttACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAAt                                                                                                                            >  8:264141/1‑149 (MQ=255)
                      aCTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAAtt                                                                                                                           <  7:13287/148‑1 (MQ=255)
                        tttCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTAcc                                                                                                                        >  5:75064/1‑149 (MQ=255)
                        tttCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTAcc                                                                                                                        <  5:248154/149‑1 (MQ=255)
                          tCTGATAAGGCGAGATTATTAAAGTTGCCATTCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTAcc                                                                                                                        <  4:218930/147‑1 (MQ=255)
                          tCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCt                                                                                                                      >  8:54515/1‑149 (MQ=255)
                             gATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATg                                                                                                                                                   >  4:266114/1‑117 (MQ=255)
                             gATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATg                                                                                                                                                   <  3:266114/117‑1 (MQ=255)
                              aTAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTc                                                                                                                    <  2:322566/147‑1 (MQ=255)
                              aTAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCa                                                                                                                   <  4:345102/148‑1 (MQ=255)
                                   gCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGc                                                                                                                                                  >  3:244553/1‑112 (MQ=255)
                                   gCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGc                                                                                                                                                  <  4:244553/112‑1 (MQ=255)
                                        attattAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCAc                                                                                                        <  3:184389/149‑1 (MQ=255)
                                          tattaAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACtt                                                                                                      >  1:103341/1‑149 (MQ=255)
                                              aaaGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGCACTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAAt                                                                                                  >  6:55262/1‑149 (MQ=255)
                                              aaaGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAAt                                                                                                  >  8:294962/1‑149 (MQ=255)
                                                  ttGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTc                                                                                               >  8:41633/1‑148 (MQ=255)
                                                     ccATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTccac                                                                                            <  8:282546/148‑1 (MQ=255)
                                                       aTGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCATAGGGTGCGACCTacac                                                                                                                                         >  3:267713/1‑101 (MQ=38)
                                                       aTGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCATAGGGTGCGACCTacac                                                                                                                                         <  4:267713/101‑1 (MQ=38)
                                                       aTGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCAt                                                                                         <  1:157705/149‑1 (MQ=255)
                                                        tGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTTTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATTCGACCTACACATCATTAAGGAATTACCCTTCACGTTGGACCACTTTAATTTCAACCATg                                                                                        >  2:83438/1‑149 (MQ=255)
                                                        tGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTaa                                                                                                   >  6:274523/1‑138 (MQ=255)
                                                        tGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTaa                                                                                                   <  5:274523/138‑1 (MQ=255)
                                                                 cGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATccc                                                                                >  7:247346/1‑148 (MQ=255)
                                                                  ggggAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATccc                                                                                >  4:127138/1‑147 (MQ=255)
                                                                  ggggAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGa                                                                              <  7:54515/149‑1 (MQ=255)
                                                                    ggAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGc                                                                            <  7:264141/149‑1 (MQ=255)
                                                                    ggAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCAGATCCCGAGc                                                                            >  6:23067/1‑149 (MQ=255)
                                                                          gtTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAATCCGTAGAAGGCGTGTAGGTGGCACCCTATTAGACCTACACATCAGTCTGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACAATGCCGATCCCGAGCTTAcgc                                                                      >  3:234013/1‑149 (MQ=255)
                                                                            tGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATAAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACgcgg                                                                    >  6:370001/1‑149 (MQ=255)
                                                                                 gcTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAgt                                                               <  8:363812/149‑1 (MQ=255)
                                                                                             gTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGatat                                                   >  7:367188/1‑149 (MQ=255)
                                                                                             gTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGatat                                                   >  7:361597/1‑149 (MQ=255)
                                                                                              tttCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTc                                                            >  8:260520/1‑139 (MQ=255)
                                                                                              tttCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTc                                                            <  7:260520/139‑1 (MQ=255)
                                                                                                    aataatGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTTCCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTccc                                               <  1:178403/146‑1 (MQ=255)
                                                                                                    aataatGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATAtt                                                  <  2:272724/143‑1 (MQ=255)
                                                                                                    aataatGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATAtt                                                  >  1:272724/1‑143 (MQ=255)
                                                                                                      taatGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCgt                                          <  8:367188/149‑1 (MQ=255)
                                                                                                      taatGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCgt                                          <  8:361597/149‑1 (MQ=255)
                                                                                                      taatGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCgt                                          >  4:75090/1‑149 (MQ=255)
                                                                                                        atGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCgttg                                        <  7:41633/149‑1 (MQ=255)
                                                                                                        atGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCgttg                                        >  1:278953/1‑149 (MQ=255)
                                                                                                        atGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCgttg                                        <  1:83438/149‑1 (MQ=255)
                                                                                                           tgttgtAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAg                                                                             >  4:128071/1‑109 (MQ=255)
                                                                                                           tgttgtAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAg                                                                             <  3:128071/109‑1 (MQ=255)
                                                                                                           tgttgtAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCgttggtt                                     <  7:294962/149‑1 (MQ=255)
                                                                                                           tgttgtAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCgttggtt                                     <  3:127138/149‑1 (MQ=255)
                                                                                                             ttgtAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCgttggttg                                    <  8:247346/148‑1 (MQ=255)
                                                                                                               gtAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGa                                                                 <  6:210683/117‑1 (MQ=255)
                                                                                                               gtAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGa                                                                 >  5:210683/1‑117 (MQ=255)
                                                                                                               gtAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATTTTCCCTTCGTTGGTTgaga                                 >  6:224494/1‑149 (MQ=255)
                                                                                                                tAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTgaga                                 <  6:75064/148‑1 (MQ=255)
                                                                                                                tAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTgaga                                 >  1:67985/1‑148 (MQ=255)
                                                                                                                    ccGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGt                            >  2:224056/1‑149 (MQ=255)
                                                                                                                       tAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTcc                         <  3:124413/149‑1 (MQ=255)
                                                                                                                        gactgatgtgtAGGTCGCATAGGGTGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCGACCACTTTAATTTCCACCATGCCGATCCAGAGCTTACGCCGCGAGTGTCCGAGGATATTCCCTTCGATGGTTGAGACAGGTTCcg                        >  6:112951/8‑149 (MQ=37)
                                                                                                                         gAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCcg                        <  5:23067/148‑1 (MQ=255)
                                                                                                                          aaGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCcggc                      >  7:104371/1‑149 (MQ=255)
                                                                                                                           tgatgtgtAGGTCGCATAGGGTGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCcggcg                     >  1:43399/5‑149 (MQ=255)
                                                                                                                            ggCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCcggcgg                    <  2:278953/149‑1 (MQ=255)
                                                                                                                                tgtAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCCGGCGGcac                 >  8:254375/1‑148 (MQ=255)
                                                                                                                                 gtAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTgaga                                 >  2:37664/1‑131 (MQ=255)
                                                                                                                                 gtAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTgaga                                 <  1:37664/131‑1 (MQ=255)
                                                                                                                                  tAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCCGGCGGCACaa               >  7:296796/1‑148 (MQ=255)
                                                                                                                                  tAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCCGGCGGCACaa               <  8:296796/148‑1 (MQ=255)
                                                                                                                                  tAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCCGGCGGCACAAt              >  2:369740/1‑149 (MQ=255)
                                                                                                                                   aGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGa                                                                 <  8:47888/97‑1 (MQ=255)
                                                                                                                                   aGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGa                                                                 >  7:47888/1‑97 (MQ=255)
                                                                                                                                   aGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCCGGCGGCACAAt              >  6:146051/1‑148 (MQ=255)
                                                                                                                                     gTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCCGGCGGCACAATTAc           >  7:84421/1‑149 (MQ=255)
                                                                                                                                         cACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCCGGCGGCACAATTACCAgt       <  4:278024/149‑1 (MQ=255)
                                                                                                                                         cACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCCGGCGGCACAATTACCAgt       <  7:254375/149‑1 (MQ=255)
                                                                                                                                         cACCCTATGCGACCTACACATCAGTCAGGAAGTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCCGGCGGCACAATTACCAgt       <  5:55262/149‑1 (MQ=255)
                                                                                                                                          aCCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCCGGCGGCACAATTACCAgtg      >  1:149753/1‑149 (MQ=255)
                                                                                                                                            ccTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCCGGCGGCACAATTACCAgtgtg    >  6:111747/1‑149 (MQ=255)
                                                                                                                                             cTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCCGGCGGCACAATTACCa         <  5:145886/143‑1 (MQ=255)
                                                                                                                                             cTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCCGGCGGCACAATTACCa         >  6:145886/1‑143 (MQ=255)
                                                                                                                                              tATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCCGGCGGCACAATTACCAgt       >  5:200358/1‑144 (MQ=255)
                                                                                                                                              tATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCCGGCGGCACAATTACCAgt       <  6:200358/144‑1 (MQ=255)
                                                                                                                                              tATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCCGGCGGCACAATTACCAGTGTGTc  >  4:316098/1‑149 (MQ=255)
                                                                                                                                               aTGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCgttg                                        >  8:137702/1‑110 (MQ=255)
                                                                                                                                               aTGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCgttg                                        <  7:137702/110‑1 (MQ=255)
                                                                                                                                               aTGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCgtt                                         >  8:108498/1‑109 (MQ=255)
                                                                                                                                               aTGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCgtt                                         <  7:108498/109‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
TGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCCGGCGGCACAATTACCAGTGTGTC  >  NC_000913/4596982‑4597272

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