Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,597,123 C→G 5.5% intergenic (+111/+27) 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,1230CG 5.5% 304.2 / 6.0 88intergenic (+111/+27)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 C (41/42);  new base G (3/2);  total (44/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 = 9.45e-01

TGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCCGGCGGCACAATTACCAGTGTG  >  NC_000913/4596982‑4597270
                                                                                                                                             |                                                                                                                                                   
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                                                                                                                                                <  4:244553/112‑1 (MQ=255)
                                   gCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGc                                                                                                                                                >  3:244553/1‑112 (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:361597/1‑149 (MQ=255)
                                                                                             gTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGatat                                                 >  7:367188/1‑149 (MQ=255)
                                                                                              tttCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTc                                                          <  7:260520/139‑1 (MQ=255)
                                                                                              tttCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTc                                                          >  8:260520/1‑139 (MQ=255)
                                                                                                    aataatGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTTCCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTccc                                             <  1:178403/146‑1 (MQ=255)
                                                                                                    aataatGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATAtt                                                <  2:272724/143‑1 (MQ=255)
                                                                                                    aataatGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATAtt                                                >  1:272724/1‑143 (MQ=255)
                                                                                                      taatGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCgt                                        >  4:75090/1‑149 (MQ=255)
                                                                                                      taatGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCgt                                        <  8:367188/149‑1 (MQ=255)
                                                                                                      taatGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCgt                                        <  8:361597/149‑1 (MQ=255)
                                                                                                        atGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCgttg                                      <  7:41633/149‑1 (MQ=255)
                                                                                                        atGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCgttg                                      >  1:278953/1‑149 (MQ=255)
                                                                                                        atGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCgttg                                      <  1:83438/149‑1 (MQ=255)
                                                                                                           tgttgtAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAg                                                                           <  3:128071/109‑1 (MQ=255)
                                                                                                           tgttgtAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAg                                                                           >  4:128071/1‑109 (MQ=255)
                                                                                                           tgttgtAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCgttggtt                                   <  3:127138/149‑1 (MQ=255)
                                                                                                           tgttgtAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCgttggtt                                   <  7:294962/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                               <  1:37664/131‑1 (MQ=255)
                                                                                                                                 gtAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTgaga                               >  2:37664/1‑131 (MQ=255)
                                                                                                                                  tAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCCGGCGGCACaa             >  7:296796/1‑148 (MQ=255)
                                                                                                                                  tAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCCGGCGGCACaa             <  8:296796/148‑1 (MQ=255)
                                                                                                                                  tAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCCGGCGGCACAAt            >  2:369740/1‑149 (MQ=255)
                                                                                                                                   aGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGa                                                               >  7:47888/1‑97 (MQ=255)
                                                                                                                                   aGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGa                                                               <  8:47888/97‑1 (MQ=255)
                                                                                                                                   aGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCCGGCGGCACAAt            >  6:146051/1‑148 (MQ=255)
                                                                                                                                     gTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCCGGCGGCACAATTAc         >  7:84421/1‑149 (MQ=255)
                                                                                                                                         cACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCCGGCGGCACAATTACCAgt     <  7:254375/149‑1 (MQ=255)
                                                                                                                                         cACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCCGGCGGCACAATTACCAgt     <  4:278024/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       >  6:145886/1‑143 (MQ=255)
                                                                                                                                             cTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCCGGCGGCACAATTACCa       <  5:145886/143‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
TGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCCGGCGGCACAATTACCAGTGTG  >  NC_000913/4596982‑4597270

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: