Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,597,121 C→T 5.7% intergenic (+109/+29) 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,1210CT 5.7% 283.0 / 5.1 88intergenic (+109/+29)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 T (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 = 8.21e-01

AACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCCGGCGGCACAATTACCAGTG  >  NC_000913/4596974‑4597268
                                                                                                                                                   |                                                                                                                                                   
aaCAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCAcc                                                                                                                                                    >  4:124413/1‑149 (MQ=255)
        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:248154/149‑1 (MQ=255)
                                tttCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTAcc                                                                                                                    >  5:75064/1‑149 (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                                                                                                                                     <  4:267713/101‑1 (MQ=38)
                                                               aTGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCATAGGGTGCGACCTacac                                                                                                                                     >  3:267713/1‑101 (MQ=38)
                                                               aTGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCAt                                                                                     <  1:157705/149‑1 (MQ=255)
                                                                tGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTTTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATTCGACCTACACATCATTAAGGAATTACCCTTCACGTTGGACCACTTTAATTTCAACCATg                                                                                    >  2:83438/1‑149 (MQ=255)
                                                                tGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTaa                                                                                               <  5:274523/138‑1 (MQ=255)
                                                                tGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTaa                                                                                               >  6:274523/1‑138 (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                                              >  1:272724/1‑143 (MQ=255)
                                                                                                            aataatGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATAtt                                              <  2:272724/143‑1 (MQ=255)
                                                                                                              taatGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCgt                                      <  8:367188/149‑1 (MQ=255)
                                                                                                              taatGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCgt                                      <  8:361597/149‑1 (MQ=255)
                                                                                                              taatGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCgt                                      >  4:75090/1‑149 (MQ=255)
                                                                                                                atGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCgttg                                    >  1:278953/1‑149 (MQ=255)
                                                                                                                atGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCgttg                                    <  1:83438/149‑1 (MQ=255)
                                                                                                                atGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCgttg                                    <  7:41633/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)
                                                                                                                                                   |                                                                                                                                                   
AACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCCGGCGGCACAATTACCAGTG  >  NC_000913/4596974‑4597268

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: ‑