Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,597,122 C→A 5.7% intergenic (+110/+28) 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,1220CA 5.7% 274.3 / 5.1 88intergenic (+110/+28)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 A (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 = 3.58e-01

AACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCCGGCGGCACAATTACCAGTGTG  >  NC_000913/4596974‑4597270
                                                                                                                                                    |                                                                                                                                                    
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                                                               <  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     <  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)
                                                                                                                                                    |                                                                                                                                                    
AACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCCGGCGGCACAATTACCAGTGTG  >  NC_000913/4596974‑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: ‑