Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,936,605 T→C 5.1% intergenic (‑291/‑47) zwf ← / → yebK glucose‑6‑phosphate 1‑dehydrogenase/putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,936,6050TC 5.1% 251.2 / 3.0 79intergenic (‑291/‑47)zwf/yebKglucose‑6‑phosphate 1‑dehydrogenase/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (36/38);  new base C (2/2);  total (39/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

TGCGGCAAAACTGATAAAAAAATCATCGTTTTCCATTGGGTAAAAATCTGACACTGATCATGTTATGAAAAAAAATAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTACATCATGAATATGCTGGAAAAAATCCAGTCTCAGCTGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATC  >  NC_000913/1936458‑1936750
                                                                                                                                                   |                                                                                                                                                 
tGCGGCAAAACTGATAAAAAAATCATCGTTTTCCATTGGGTAAAAATCTGACACTGATCATGTTATGAAAAAAAATAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGAtt                                                                                                                                                  >  7:187467/1‑149 (MQ=255)
   ggCAAAACTGATAAAAAAATCATCGTTTTCCATTGGGTAAAAATCTGACACTTATCATGTTATGAAAAAAAATAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGAtttct                                                                                                                                               <  8:28572/149‑1 (MQ=255)
    gCAAAACTGATAAAAAAATCATCGTTTTCCATTGTGTAAAAATCTGACACTGATCATGTTATGAAAAAAAATAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGAtttctt                                                                                                                                              <  5:247970/149‑1 (MQ=255)
    gCAAAACTGATAAAAAAATCATCGTTTTCCATTGGGTAAAAATCTGACACTGATCATGTTATGAAAAAAAATAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGAtttctt                                                                                                                                              >  4:33056/1‑149 (MQ=255)
     cAAAACTGATAAAAAAATCATCGTTTTCCATTGGGTAAAAATCTGACACTGATCATGTTATGAAAAAAAATAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGAtttcttt                                                                                                                                             <  6:188467/149‑1 (MQ=255)
       aaaCTGATAAAAAAATCTTCGTTTTCCATTGGGTAAAAATCTGACACTGATCATGTTATGAAAAAAAATAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCa                                                                                                                                           <  2:216315/149‑1 (MQ=255)
       aaaCTGATAAAAAAATCATCGTTTTCCATTGGGTAAAAATCTGACACTGATCATGTTATGAAAAAAAATAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGAtttcttt                                                                                                                                             >  2:128002/1‑147 (MQ=255)
       aaaCTGATAAAAAAATCATCGTTTTCCATTGGGTAAAAATCTGACACTGATCATGTTATGAAAAAAAATAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGAtttcttt                                                                                                                                             <  1:128002/147‑1 (MQ=255)
        aaCTGATAAAAAAATCATCGTTTTCCATTGGGTAAAAATCTGACACTGATCATGTTATGAAAAAAAATAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGAtttcttt                                                                                                                                             >  1:97318/1‑146 (MQ=255)
        aaCTGATAAAAAAATCATCGTTTTCCATTGGGTAAAAATCTGACACTGATCATGTTATGAAAAAAAATAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCa                                                                                                                                           <  2:97318/148‑1 (MQ=255)
                  aaaaTCATCGTTTTCCATTGGGTAAAAATCTGACACTGATCATGTTATGAAAAAAAATAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATc                                                                                                                                <  8:205164/149‑1 (MQ=255)
                     atcatcGTTTTCCATTGGGTAAAAATCTGACACTGATCATGTTATGAAAAAAAATAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTc                                                                                                                             <  2:68165/149‑1 (MQ=255)
                     atcatcGTTTTCCATTGGGTAAAAATCTGACACTGATCATGTTATGAAAAAAAATAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTc                                                                                                                             <  7:66905/149‑1 (MQ=255)
                      tcatcGTTTTCCATTGGGTAAAAATCTGACACTGATCATGTTATGAAAAAAAATAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCa                                                                                                                            <  7:1074/149‑1 (MQ=255)
                      tcatcGTTTTCCATTGGGTAAAAATCTGACACTGATCATGTTATGAAAAAAAATAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCa                                                                                                                            <  2:225644/149‑1 (MQ=255)
                       catcGTTTTCCATTGGGTAAAAATCTGACACTGATCATGTTATGAAAAAAAATAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCAt                                                                                                                           >  1:280591/1‑149 (MQ=255)
                        atcGTTTTCCATTGGGTAAAAATCTGACACTGATCATGTTATGAAAAAAAATAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTTATTTCTGAGAAAACTGCATTTCGATTTCGTTCAGGGGGGAACTCGTCAtt                                                                                                                          >  4:285802/1‑149 (MQ=255)
                         tcGTTTTCCATTGGGTAAAAATCTGACACTGATCATGTTATGAAAAAAAATAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGAtttcttt                                                                                                                                             <  1:345366/129‑1 (MQ=255)
                         tcGTTTTCCATTGGGTAAAAATCTGACACTGATCATGTTATGAAAAAAAATAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGAtttcttt                                                                                                                                             >  2:345366/1‑129 (MQ=255)
                         tcGTTTTCCATTGGGTAAAAATCTGACACTGATCATGTTATCAAAAAAAATAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTa                                                                                                                         >  3:150199/1‑149 (MQ=255)
                            ttttCCATTGGGTAAAAATCTGACACTGATCATGTTATGAAAAAAAATAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTAccc                                                                                                                      >  6:357107/1‑149 (MQ=255)
                            ttttCCATTGGGTAAAAATCTGACACTGATCATGTTATGAAAAAAAATAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTAccc                                                                                                                      <  8:126781/149‑1 (MQ=255)
                              ttCCATTGGGTAAAAATCTGACACTGATCATGTTATGAAAAAAAATAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGt                                                                                                                    >  6:78686/1‑149 (MQ=255)
                                 cATTGGGTAAAAATCTGCCACTGATCATGTTATGAAAAAAAATAACAACTTTTTTATCTGCTTTTGTCATTAACTCGGCTCACAGGCGTAATATCTGACAAAACTGCATTTCGAGGTCTTGCAGTGCGGAAAACCTCATTACCCGTGAg                                                                                                                 >  8:344507/1‑149 (MQ=255)
                                 cATTGGGTAAAAATCTGACACTGATCATGTTATCAAAAAAAATAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGa                                                                                                                  <  4:150199/148‑1 (MQ=255)
                                  aTTGGGTAAAAATCTGACACTGATCATGTTATGAAAAAAAATAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGt                                                                                                                >  1:155793/1‑149 (MQ=255)
                                        tAAAAATCTGACACTGATCATGTTATGAAAAAAAATAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCtt                                                                                                           >  5:268601/1‑148 (MQ=255)
                                           aaaTCTGACACTGATCATGTTATGAAAAAAAATAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTAca                                                                                                       <  4:184296/149‑1 (MQ=255)
                                              tCTGACACTGATCATGTTATGAAAAAAAATAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTAcatca                                                                                                    >  2:210831/1‑149 (MQ=255)
                                                    acTGATCATGTTATGAAAAAAAATAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTACATCATGAata                                                                                              >  7:342666/1‑149 (MQ=255)
                                                      tGATCATGTTATGAAAAAAAATAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTACATCATGAATATg                                                                                            >  7:361964/1‑149 (MQ=255)
                                                       gATCATGTTATGAAAAAAAATAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTACATCATGAat                                                                                               >  5:341060/1‑145 (MQ=255)
                                                       gATCATGTTATGAAAAAAAATAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTACATCATGAat                                                                                               <  6:341060/145‑1 (MQ=255)
                                                        aTCATGTTATGAAAAAAAATAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGt                                                                                                                              >  5:28295/1‑113 (MQ=255)
                                                        aTCATGTTATGAAAAAAAATAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGt                                                                                                                              <  6:28295/113‑1 (MQ=255)
                                                        aTCATGTTATGAAAAAAAATAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGtctc                                                                                                             >  1:314185/1‑130 (MQ=255)
                                                        aTCATGTTATGAAAAAAAATAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTAcatc                                                                                                     <  2:28136/138‑1 (MQ=255)
                                                        aTCATGTTATGAAAAAAAATAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTAcatc                                                                                                     >  1:28136/1‑138 (MQ=255)
                                                        aTCATGTTATGAAAAAAAATAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTACATCATGAATATGCt                                                                                          >  6:220745/1‑149 (MQ=255)
                                                            tGTTATGAAAAAAAATAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTACATCATGAATATGCTGGaa                                                                                      <  2:155793/149‑1 (MQ=255)
                                                             gTTATGAAAAAAAATAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTAAGTGAGGAAATCGt                                                                                                                              >  2:264028/1‑108 (MQ=255)
                                                             gTTATGAAAAAAAATAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAAAATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGt                                                                                                                              <  1:264028/108‑1 (MQ=255)
                                                                 tGAAAAAAAATAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTACATCATGAATATGCTGGaaaaaa                                                                                  <  5:357107/148‑1 (MQ=255)
                                                                   aaaaaaaaTAACAACTTTTTTATCTTCTTTTGTCATTAACACGGGACACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTTAGTCTCTTTACATCATGAATATGCTGGAAAAAATc                                                                                <  5:276453/148‑1 (MQ=255)
                                                                   aaaaaaaaTAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGAGAAAACTGCATTTCTATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTACATCATGAATATGCTGGAAAAAATc                                                                                >  5:93789/1‑148 (MQ=255)
                                                                   aaaaaaaaTAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTACATCATGAATATGCTGGAAAAAATc                                                                                <  5:220745/148‑1 (MQ=255)
                                                                   aaaaaaaaTAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTACATCATGAATATGCTGGAAAAAATc                                                                                <  8:378647/148‑1 (MQ=255)
                                                                    aaaaaaaTAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTACATCATGAATATGCTGGAAAAAATCCa                                                                              <  2:280591/149‑1 (MQ=255)
                                                                     aaaaaaTAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTACATCATGAATATGCTGGAAAAAATCCAg                                                                             <  8:342666/149‑1 (MQ=255)
                                                                      aaaaaTAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTACATCATGAATATGCTGGAAAAAATCCAGt                                                                            <  3:285802/149‑1 (MQ=255)
                                                                       aaaaTAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTACATCATGAATATGCTGGAAAAAATCCAGtc                                                                           >  1:326065/1‑149 (MQ=255)
                                                                          aTAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGACTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTACATCATGAATATGCTGGAAAAAATCCAGtct                                                                          >  5:245380/1‑147 (MQ=255)
                                                                          aTAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGACTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTACATCATGAATATGCTGGAAAAAATCCAGtct                                                                          >  5:250758/1‑147 (MQ=255)
                                                                           tAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGACTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTACATCATGAATATGCTGGAAAAAATCCAGTCTCAg                                                                       <  6:245380/149‑1 (MQ=255)
                                                                           tAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGACTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTACATCATGAATATGCTGGAAAAAATCCAGTCTCAg                                                                       <  6:250758/149‑1 (MQ=255)
                                                                             acaacTTTTTTATGTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGTAAATCGTCATTACCCGTGAGTCTCTTTACATCATGAATATGCTGTAAAAAATCCAGTCTCAGCt                                                                     <  8:187467/149‑1 (MQ=255)
                                                                             acaacTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCgg                                                                                                                                     <  2:230548/85‑1 (MQ=255)
                                                                             acaacTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCgg                                                                                                                                     >  1:230548/1‑85 (MQ=255)
                                                                               aacTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTACATCATGAATATGCTGGAAAAAATCCAGTCTCAGCTg                                                                    >  4:176669/1‑148 (MQ=255)
                                                                                       tATCTGCTTTTGTCATTAACACGGCTCACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTACATCATGAATATGCTGGAAAAAATCCAGTCTCAGCTGGAACAttt                                                            <  4:48314/148‑1 (MQ=255)
                                                                                        aTCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTACATCATGAATATGCTGGAAAAAATCCAGTCTCAGCTGGAACATTTg                                                           <  6:268601/148‑1 (MQ=255)
                                                                                        aTCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTACATCATGAATATGCTGGAAAAAATCCAGTCTCAGCTGGAACATTTg                                                           <  1:337438/148‑1 (MQ=255)
                                                                                         tCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTACATCATGAATATGCTGGAAAAAATCCAGTCTCAGCTGGAACATTTGa                                                          <  3:33056/148‑1 (MQ=255)
                                                                                                  gTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTACATCATGAATATGCTGGAAAAAATCCAGTCTCAGCTGGAACATTTGAGCAAATCa                                                  >  1:191584/1‑147 (MQ=255)
                                                                                                   tCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTACATCATGAATATGCTGGAAAAAATCCAGTCTCAGCTGGAACATTTGAGCAAATCaga                                                <  5:78686/148‑1 (MQ=255)
                                                                                                     aTTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTACATCATGAATATGCTGGAAAAAATCCAGTCTCAGCTGGAACATTTGAGCAAATCAGAgcg                                             >  1:76613/1‑149 (MQ=255)
                                                                                                                 cacaGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTTTTTACATCATGAATATGCTGGAAAAAATCCAGTCTCAGCTGGAACATTTGa                                                          <  5:351905/124‑1 (MQ=255)
                                                                                                                 cacaGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTATTTACATCATGAATATGCTGGAAAAAATCCAGTCTCAGCTGGAACATTTGa                                                          >  6:351905/1‑124 (MQ=255)
                                                                                                                  acaGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTACATCATGAATATGCTGGAAAAAATCCAGTCTCAGCTg                                                                    <  8:50229/113‑1 (MQ=255)
                                                                                                                    aGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTACATCATGAATATGCTGGAAAAAATCCAGTCTCAGCTGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGt                              >  7:28470/1‑149 (MQ=255)
                                                                                                                    aGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTACATCATGAATATGCTGGAAAAAATCCAGTCTCAGCTGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGt                              >  7:131032/1‑149 (MQ=255)
                                                                                                                      gCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTACATCATGAATATGCTGGAAAAAATCCAGTCTCAGCTGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCa                            >  2:393179/1‑149 (MQ=255)
                                                                                                                        gTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTACATCATGAATATGCTGGAAAAAATCCAGTCTCAGCTGGAACATTTGAGCaaa                                                     >  5:210196/1‑122 (MQ=255)
                                                                                                                        gTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTACATCATGAATATGCTGGAAAAAATCCAGTCTCAGCTGGAACATTTGAGCaaa                                                     <  6:210196/122‑1 (MQ=255)
                                                                                                                        gTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTACATCATGAATATGCTGGAAAAAATCCAGTCTCAGCTGGAACATTTGAGCAAATCAGAGCTCAAAGTTGCCGAGGTCAtt                          >  8:163574/1‑149 (MQ=255)
                                                                                                                         tAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTACATCATGAATATGCTGGAAAAAATCCAGTCTCAGCTGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTc                         >  4:97728/1‑149 (MQ=255)
                                                                                                                           atatCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTACATCATGAATATGCTGGAAAAAATCCAGtctc                                                                         >  8:134171/1‑99 (MQ=255)
                                                                                                                           atatCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTACATCATGAATATGCTGGAAAAAATCCAGtctc                                                                         <  7:134171/99‑1 (MQ=255)
                                                                                                                           atatCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTACATCATGAATATGCTGGAAAAAATCCAGTCTCAGCTGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCt                        <  3:176669/148‑1 (MQ=255)
                                                                                                                              tCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAATTCGTCATTACCCGTGAGTCTCTTTACATCATGAATATGCTGGAAAAAATCCAGTCTCAGCTGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCt                    >  1:65097/1‑149 (MQ=255)
                                                                                                                              tCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTACATCATGAATATGCTGGAAAAAATCCAGTCTCAGCTGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGc                     <  2:76613/148‑1 (MQ=255)
                                                                                                                                     aaaCTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTACATCATGAATATGCTGGAAAAAATCCAGTCTCAGCTGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGccc              >  6:42269/1‑148 (MQ=255)
                                                                                                                                         tGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTAcatc                                                                                                     >  8:111118/1‑57 (MQ=255)
                                                                                                                                         tGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTAcatc                                                                                                     <  7:111118/57‑1 (MQ=255)
                                                                                                                                         tGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTACATCATGAATATGCTGGAAAAAATCCAGTCTCAGCTGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATaa         >  1:313324/1‑149 (MQ=255)
                                                                                                                                            aTTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTACATCATGAATATGCTGGAAAAAATCCAGTCTCAGCTGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAAcgc      <  1:393179/149‑1 (MQ=255)
                                                                                                                                             tttCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTACATCATGAATATGCTGGAAAAAATCCAGTCTCAGCTGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGa                                 <  8:251514/121‑1 (MQ=255)
                                                                                                                                             tttCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTACATCATGAATATGCTGGAAAAAATCCAGTCTCAGCTGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGa                                 >  7:251514/1‑121 (MQ=255)
                                                                                                                                                cGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTACATCATGAATATGCTGGAAAAAATCCAGTCTCAGCTGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATc  >  4:132061/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TGCGGCAAAACTGATAAAAAAATCATCGTTTTCCATTGGGTAAAAATCTGACACTGATCATGTTATGAAAAAAAATAACAACTTTTTTATCTGCTTTTGTCATTAACACGGCACACAGGCGTAATATCTGACAAAACTGCATTTCGATTTCTTTCAGTGCGGAAATCGTCATTACCCGTGAGTCTCTTTACATCATGAATATGCTGGAAAAAATCCAGTCTCAGCTGGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATC  >  NC_000913/1936458‑1936750

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: