Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,873,084 T→C S177G (AGC→GGC)  dgoA ← 2‑dehydro‑3‑deoxy‑6‑phosphogalactonate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,873,0840TC99.4% 1108.3 / ‑5.6 331S177G (AGC→GGC) dgoA2‑dehydro‑3‑deoxy‑6‑phosphogalactonate aldolase
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (179/148);  total (179/152)
Fisher's exact test for biased strand distribution p-value = 2.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 353 total aligned reads displayed.

GAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAG  >  NC_000913/3872948‑3873218
                                                                                                                                        |                                                                                                                                      
      ccAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAgcc                                                                                                                               <  2:5354965/140‑1 (MQ=255)
      ccAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAgcc                                                                                                                               <  2:5225852/140‑1 (MQ=255)
          cGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAgcccgcc                                                                                                                           <  1:532023/140‑1 (MQ=255)
          cGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAgcccgcc                                                                                                                           <  1:6854872/140‑1 (MQ=255)
           ggggAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAgcccgccc                                                                                                                          <  1:5076244/140‑1 (MQ=255)
           ggggAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAgcccgccc                                                                                                                          <  2:606188/140‑1 (MQ=255)
              gAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTg                                                                                                                       >  1:7247760/1‑140 (MQ=255)
                ggTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTTTCGATCCGCTTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGca                                                                                                                     <  2:5417547/140‑1 (MQ=255)
                ggTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGca                                                                                                                     <  2:478434/140‑1 (MQ=255)
                    aaCGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAAc                                                                                                                 >  1:5919488/1‑140 (MQ=255)
                    aaCGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAAc                                                                                                                 <  1:6617278/140‑1 (MQ=255)
                       gATACGTGTTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTg                                                                                                              <  2:1398636/140‑1 (MQ=255)
                       gATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTg                                                                                                              <  2:7116416/140‑1 (MQ=255)
                       gATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTg                                                                                                              <  1:4933266/140‑1 (MQ=255)
                       gATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTg                                                                                                              <  2:6343569/140‑1 (MQ=255)
                       gATACGTGGTAATTTTGCTGATCTTCATTGCCCTGTCTCCCGATACGCATTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTg                                                                                                              <  2:2239655/140‑1 (MQ=255)
                         tACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGTTAGCGATGGCCGCCTAAGCCCTCCCCTGCACAACCTGCg                                                                                                            >  2:3023910/1‑140 (MQ=255)
                         tACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCg                                                                                                            >  1:5737283/1‑140 (MQ=255)
                           cGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTc                                                                                                          >  2:6218698/1‑140 (MQ=255)
                            gTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCt                                                                                                         >  2:4964864/1‑140 (MQ=255)
                              ggTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTAt                                                                                                       >  2:6017306/1‑140 (MQ=255)
                                 aaTTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCa                                                                                                    >  2:1603146/1‑140 (MQ=255)
                                  aTTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCAc                                                                                                   <  1:1208834/140‑1 (MQ=255)
                                   ttttGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAAcc                                                                                                                >  1:649720/1‑126 (MQ=255)
                                   ttttGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACt                                                                                                  >  2:3771447/1‑140 (MQ=255)
                                       ggTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGACATCGCCGCTTAAGCCCGCCCCTGCACAACCAGCGTCTATCCACTgcgc                                                                                              >  2:2479231/1‑140 (MQ=255)
                                       ggTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGACATCGCCGCGTAAGCCCGCCCCTGCACAACCCGCGTCTATCCACTgcgc                                                                                              >  2:3099238/1‑140 (MQ=255)
                                       ggTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAAATCGCCGCCTAACCCCGCCCCTGCACAACCCGCGTCTATCCACTgcgc                                                                                              >  2:3253728/1‑140 (MQ=255)
                                        gTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGCTCGCCGCCCAAGGCCGTCCCTGCACGACCTGCGTCTATCCACTGCGcc                                                                                             >  2:1990389/1‑140 (MQ=255)
                                        gTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGcc                                                                                             >  2:1126457/1‑140 (MQ=255)
                                          gATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAg                                                                                           >  1:4389449/1‑140 (MQ=255)
                                             tttCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGtt                                                                                        >  2:4066872/1‑140 (MQ=255)
                                                          ccTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACGACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCAcg                                                                           >  2:4881052/1‑140 (MQ=255)
                                                          ccTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACTTGCGTCTTTCCACTACGCCAGGTTTTCTGGCGTCAcg                                                                           >  1:2879443/1‑140 (MQ=255)
                                                          ccTCTCGATACGCCTTAACAAATGCTGCTGCCTCCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCAcg                                                                           >  1:6767370/1‑140 (MQ=255)
                                                            tctcGATACGCCTTAACAAATGCTGCTGCGTGCTGCGTGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCAcgcc                                                                         <  1:3068889/140‑1 (MQ=255)
                                                                gATAGGCGTTAACAAATGCTGCTGCCTGCGGCGCGTTGCGCTTTACGGATTCCCCGGGGCGATAGGGACCTCCGCCTAAGCCCGCCCCTGCACAACCGTCGTCTAGCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCa                                                                     <  2:1847685/140‑1 (MQ=255)
                                                                gATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCa                                                                     <  2:5849225/140‑1 (MQ=255)
                                                                gATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAAGTTTTCTGGCGTCACGCCGCCa                                                                     <  2:68529/140‑1 (MQ=255)
                                                                agacagCCTTAACAAGTGCTGCGGCCTGCCGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCa                                                                     <  1:3654527/135‑1 (MQ=255)
                                                                  tACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCTGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAAc                                                                   <  2:1803408/140‑1 (MQ=255)
                                                                    cGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATGCACTGCGTCAGGTTGTCTGGCGTCACGCCGCCAACgg                                                                 >  1:2073410/1‑140 (MQ=255)
                                                                    cGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACgg                                                                 >  2:4585863/1‑140 (MQ=255)
                                                                     gCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGTGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGc                                                                <  1:1041447/140‑1 (MQ=255)
                                                                     gCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGc                                                                <  1:7437476/140‑1 (MQ=255)
                                                                     gCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGc                                                                <  2:7245510/140‑1 (MQ=255)
                                                                     gCCTTAACAAATGCTGCTGCCTCCTGCGCCGTGCGCTGTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGc                                                                <  1:2721427/140‑1 (MQ=255)
                                                                     gCCTTAACAAATGCTGCGGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGc                                                                <  1:1882481/140‑1 (MQ=255)
                                                                     gCCTTAACAAATGCTGCGGCCTCCTGCGCGGAGCGCTCTACGGATTGCTCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACggt                                                                <  2:2879443/140‑2 (MQ=255)
                                                                     cccTTAACGAATGCTGCCGCCTTCTGCGCGGCGCGGTCTGCGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGc                                                                <  2:3445970/139‑1 (MQ=255)
                                                                     cccTTAACACATGCTGCTGCATGCTGCGGGGTGCGCTCTGCGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGc                                                                <  2:826810/139‑1 (MQ=255)
                                                                       cTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCaa                                                              >  2:4342469/1‑140 (MQ=255)
                                                                        ttAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGc                                                                >  2:6789563/1‑137 (MQ=255)
                                                                          aaCAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGc                                                                >  1:1839428/1‑135 (MQ=255)
                                                                               aTGCTGTTGCTGGCTGCGGGGTGGGCGCTAGGGTTCGCCGGGCGCGATGGAGATCACCGCCTAACCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCg                                                      <  2:5316332/140‑1 (MQ=255)
                                                                               aTGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCg                                                      <  1:1603146/140‑1 (MQ=255)
                                                                                tgctgctgcCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAAc                                                                   >  1:4855414/1‑126 (MQ=255)
                                                                                tgctgctgcCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGa                                                     <  2:5919488/140‑1 (MQ=255)
                                                                                 gctgctgcCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCCCTGCGCCAGGTTATCTGGCGTCACGGCGCCAGCGGCAAAGAcctgtcc                                                    >  2:3654527/1‑134 (MQ=255)
                                                                                 gctgctgcCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGAt                                                    >  2:6449954/1‑140 (MQ=255)
                                                                                 gctgctgcCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGAt                                                    >  1:914126/1‑140 (MQ=255)
                                                                                   tgctgcctgctgcGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGt                                                  <  1:1409415/140‑1 (MQ=255)
                                                                                     ctgcctgctgcGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTcc                                                <  1:3245716/140‑1 (MQ=255)
                                                                                     ctgcctgctgcGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTcc                                                >  1:5312357/1‑140 (MQ=255)
                                                                                     ctgcctgctgcGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTcc                                                <  1:5618550/140‑1 (MQ=255)
                                                                                     ctgcctgctgcGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTcc                                                <  1:4402913/140‑1 (MQ=255)
                                                                                     ctgcctgctgcGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTATGGCGTCCCGCGGCAAACGGCAAAGACTGCGATGTcc                                                >  1:2180367/1‑140 (MQ=255)
                                                                                     ctgcctgctgcGCGGTCCGCTCGACGGATTGACCGTCGCGATAGAGATGGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTcc                                                <  2:3411143/140‑1 (MQ=255)
                                                                                       gcctgctgcGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGa                                              >  1:6998873/1‑140 (MQ=255)
                                                                                        cctgctgcGCGGTGCGCTCTACGGATTGCTCGGCGCGCTAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGAt                                             <  1:2561111/140‑1 (MQ=255)
                                                                                            ctgcGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCa                                         >  1:5958127/1‑140 (MQ=255)
                                                                                            ctgcGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCa                                         >  1:3282870/1‑140 (MQ=255)
                                                                                               cgcgGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATa                                      >  1:4701225/1‑140 (MQ=255)
                                                                                               cgcgGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATa                                      >  2:4433826/1‑140 (MQ=255)
                                                                                                gcgGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATAc                                     <  1:3835786/140‑1 (MQ=255)
                                                                                                gcgGTGCGCTCTACGAATTGCCTGGCGCGATAGAGATGGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATAc                                     <  1:2948691/140‑1 (MQ=255)
                                                                                                      cgcTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCAAACAGCAAAGACTGCGTTGTCCGATGGCAATACCGCttt                               >  2:2681590/1‑140 (MQ=255)
                                                                                                       gcTCTACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAAGGGCAAAGACTGAGACGTCCGACGGCAATACCGCtttt                              >  1:3423262/1‑140 (MQ=255)
                                                                                                        ctctACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTa                             <  2:7353068/140‑1 (MQ=255)
                                                                                                        ctctACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTa                             <  1:791920/140‑1 (MQ=255)
                                                                                                         tctcCGGATTGCCCGCCGCGATCGAGATCGCCGCCTACGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTaa                            <  1:3099238/140‑1 (MQ=255)
                                                                                                         tctACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTaa                            <  1:4066872/140‑1 (MQ=255)
                                                                                                         tctACGGATGGCCCGGCGCGATTGAGATCGCGGCCAAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCGATACCGCTTTTaa                            <  2:2532378/140‑1 (MQ=255)
                                                                                                           tACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACg                          <  2:1414439/140‑1 (MQ=255)
                                                                                                           tACGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACg                          <  2:5737283/140‑1 (MQ=255)
                                                                                                             cGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCt                        <  2:4471880/140‑1 (MQ=255)
                                                                                                             cGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCt                        >  1:7498063/1‑140 (MQ=255)
                                                                                                             cGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCt                        <  2:6664778/140‑1 (MQ=255)
                                                                                                             cGGATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCt                        <  2:4180846/140‑1 (MQ=255)
                                                                                                              ggATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGACAACGGCATAGAGTGCAATGTCCGAAGGCAGTACCGCTTTTAACGCtt                       >  2:1944152/1‑140 (MQ=255)
                                                                                                               gATTGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCttt                      <  1:4208398/140‑1 (MQ=255)
                                                                                                                  tGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTa                             <  1:5493453/130‑1 (MQ=255)
                                                                                                                  tGCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTa                             >  2:5493453/1‑130 (MQ=255)
                                                                                                                   gCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATg                  >  2:3294/1‑140 (MQ=255)
                                                                                                                   gCCCGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACAGCAAAGACGGCGCTGTCCGATGGCAATACCACTTTTAACGCTTTGATg                  >  2:2307300/1‑140 (MQ=255)
                                                                                                                     ccGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTa                >  2:7531674/1‑140 (MQ=255)
                                                                                                                     ccGGCGCGATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTa                >  2:5179162/1‑140 (MQ=255)
                                                                                                                          gcgATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCGGATGGCAATACGGCGTTTAACGCTTTGATGTATTGCg           >  1:1955585/1‑140 (MQ=255)
                                                                                                                          gcgATAGAGATCGCCGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCg           >  1:6359416/1‑140 (MQ=255)
                                                                                                                          gcgATAGAGATCGCCCCCTAAGACCCCCCCTGCACAAACTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCg           <  1:2331050/140‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                      
GAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAG  >  NC_000913/3872948‑3873218

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: