Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,975,282 G→C 5.2% S46T (AGC→ACC) ‡ wecE → TDP‑4‑oxo‑6‑deoxy‑D‑glucose transaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,975,2820GC 5.2% 409.5 / 9.0 116S46T (AGC→ACC) ‡wecETDP‑4‑oxo‑6‑deoxy‑D‑glucose transaminase
Reads supporting (aligned to +/- strand):  ref base G (56/54);  new base C (4/2);  total (60/56)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 116 total aligned reads displayed.

CATGATTCCATTTAACGCACCGCCGGTGGTGGGAACCGAACTCGACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAA  >  NC_000913/3975145‑3975428
                                                                                                                                         |                                                                                                                                                  
caTGATTCCATTTAACGCACCGCCGGTGGTGGGAACCGAACTCGACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAgtgt                                                                                                                                         >  2:770369/1‑149 (MQ=255)
 aTGATTCCATTTAACGCACCGCCGGTGGTGGGAACCGAACTCGACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAgtgtg                                                                                                                                        >  1:1294394/1‑148 (MQ=255)
 aTGATTCCATTTAACGCACCGCCGGTGGTGGGAACCGAACTCGACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGtt                                                                                                                                        >  2:1262705/1‑149 (MQ=255)
 aTGATTCCATTTAACGCACCGCCGGTGGTGGGAACCGAACTCGACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGAGCCAAAGTGtt                                                                                                                                        <  1:770369/149‑1 (MQ=255)
     ttCCATTTAACGCACCGCCGGTGGTGGGAACCGAACTCGACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTg                                                                                                                                    >  2:818448/1‑149 (MQ=255)
          tttAACGCACCGCCGGTGGTGGGAACCGAACTCGACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTa                                                                                                                                       <  4:734080/141‑1 (MQ=255)
          tttAACGCACCGCCGGTGGTGGGAACCGAACTCGACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGcc                                                                                                                               >  1:618246/1‑149 (MQ=255)
          tttAACGCACCGCCGGTGGTGGGAACCGAACTCGACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGcc                                                                                                                               >  1:560942/1‑149 (MQ=255)
              aCGCACCGCCGGTGGTGGGAACCGAACTCGACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTcc                                                                                                                           <  1:818448/149‑1 (MQ=255)
              aCGCACCGCCGGTGGTGGGAACCGAACTCGACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTcc                                                                                                                           <  4:584003/149‑1 (MQ=255)
               cGCACCGCCGGTGGTGGGAACCGAACTCGACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCaaa                                                                                                                                             <  4:1367724/130‑1 (MQ=255)
               cGCACCGCCGGTGGTGGGAACCGAACTCGACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCaaa                                                                                                                                             >  3:1367724/1‑130 (MQ=255)
                gCACCGCCGGTGGTGGGAACCGAACTCGACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCGCTGCCaaaac                                                                                                                                           >  2:928667/1‑129 (MQ=255)
                gCACCGCCGGTGGTGGGAACCGAACTCGACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCGCTGCCaaaa                                                                                                                                            <  2:694980/130‑2 (MQ=255)
                gCACCGCCGGTGGTGGGAACCGAACTCGACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCGCTGCCaaaa                                                                                                                                            >  1:694980/1‑129 (MQ=255)
                gCACCGCCGGTGGTGGGAACCGAACTCGACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTg                                                                                                                         >  2:427832/1‑149 (MQ=255)
                 cACCGCCGGTGGTGGGAACCGAACTCGACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGc                                                                                                                        >  1:1058137/1‑149 (MQ=255)
                        ggtggtGGGAACCGAACTCGACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCtt                                                                                                                 <  1:1313635/149‑1 (MQ=255)
                              gggAACCGAACTCGACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCaaa                                                                                                                                             >  3:1176990/1‑115 (MQ=255)
                              gggAACCGAACTCGACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCaaa                                                                                                                                             <  4:1176990/115‑1 (MQ=255)
                              gggAACCGAACTCGACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTgg                                                                                                           <  1:427832/149‑1 (MQ=255)
                               ggAACCGAACTCGACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGga                                                                                                          <  3:707667/149‑1 (MQ=255)
                               ggAACCGAACTCGACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGga                                                                                                          <  2:1344646/149‑1 (MQ=255)
                                 aaCCGAACTCGACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTc                                                                                                                <  2:1280762/141‑1 (MQ=255)
                                  aCCGAACTCGACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTc                                                                                                                <  2:302179/140‑1 (MQ=255)
                                  aCCGAACTCGACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTc                                                                                                                >  1:302179/1‑140 (MQ=255)
                                  aCCGAACTCGACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGAt                                                                                                       <  2:1294394/149‑1 (MQ=255)
                                    cGAACTCGACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCGCTGCCaaaac                                                                                                                                           >  2:1076566/1‑109 (MQ=255)
                                    cGAACTCGACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCGCTGCCaaaac                                                                                                                                           >  4:245259/1‑109 (MQ=255)
                                    cGAACTCGACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGAGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCGCTGCCaaaac                                                                                                                                           <  1:1076566/111‑3 (MQ=255)
                                     gAACTCGACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATggc                                                                                                    <  3:704542/149‑1 (MQ=255)
                                         tCGACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAgtg                                                                                                                                          <  1:812838/107‑1 (MQ=255)
                                         tCGACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAgtg                                                                                                                                          >  2:812838/1‑107 (MQ=255)
                                           gACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCgct                                                                                              >  1:506338/1‑149 (MQ=255)
                                           gACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCgct                                                                                              >  4:761264/1‑149 (MQ=255)
                                            aCTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCgctg                                                                                             <  2:618246/149‑1 (MQ=255)
                                              tataTGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCgctgct                                                                                           <  1:841481/149‑1 (MQ=255)
                                                tatGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCgctgctgc                                                                                         >  1:355550/1‑149 (MQ=255)
                                                     aGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGat                                                                                    <  3:761264/149‑1 (MQ=255)
                                                      gTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTc                                                                                       >  2:725658/1‑145 (MQ=255)
                                                      gTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTc                                                                                       <  1:725658/145‑1 (MQ=255)
                                                      gTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGata                                                                                   >  2:1337708/1‑149 (MQ=255)
                                                            aaTGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCgctgctgct                                                                                        <  4:866362/138‑1 (MQ=255)
                                                            aaTGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCgctgctgct                                                                                        >  3:866362/1‑138 (MQ=255)
                                                            aaTGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGc                                                                             <  4:1392751/149‑1 (MQ=255)
                                                            aaTGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGc                                                                             >  1:1148204/1‑149 (MQ=255)
                                                              tGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCt                                                                           >  3:449582/1‑149 (MQ=255)
                                                               gggTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTg                                                                          <  4:449582/149‑1 (MQ=255)
                                                                 gTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACt                                                                                                                                     <  3:279062/88‑1 (MQ=255)
                                                                 gTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACt                                                                                                                                     >  4:279062/1‑88 (MQ=255)
                                                                      ggCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGa                                                                                                                                   <  3:1294214/85‑1 (MQ=255)
                                                                       gCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCt                                                                                                                          <  2:560043/93‑1 (MQ=255)
                                                                       gCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGaa                                                                  <  4:842660/149‑1 (MQ=255)
                                                                        cAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTCGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAg                                                                 <  2:1148204/149‑1 (MQ=255)
                                                                            cTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACt                                                                                                                                     <  4:243574/77‑1 (MQ=255)
                                                                            cTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACt                                                                                                                                     >  3:243574/1‑77 (MQ=255)
                                                                            cTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGAt                                                             <  2:506338/149‑1 (MQ=255)
                                                                                  ggCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGcc                                                       >  4:1184095/1‑149 (MQ=255)
                                                                                       tGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCt                                                  <  1:1337708/149‑1 (MQ=255)
                                                                                        ggCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGc                                                        >  1:1095694/1‑142 (MQ=255)
                                                                                        ggCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTa                                                 >  3:473666/1‑149 (MQ=255)
                                                                                         gCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGata                                                                                   <  3:713481/114‑1 (MQ=255)
                                                                                         gCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGata                                                                                   >  4:713481/1‑114 (MQ=255)
                                                                                           ggTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTacac                                              >  2:701960/1‑149 (MQ=255)
                                                                                           ggTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTacac                                              >  3:852454/1‑149 (MQ=255)
                                                                                            gTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGa                                                              <  2:121821/132‑1 (MQ=255)
                                                                                            gTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGa                                                              >  1:121821/1‑132 (MQ=255)
                                                                                             ttttACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCt                                            <  1:1262705/149‑1 (MQ=255)
                                                                                                  cccGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGtc                                       <  2:1058137/149‑1 (MQ=255)
                                                                                                   ccGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCa                                                                                                                       <  3:984591/68‑1 (MQ=255)
                                                                                                   ccGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCt                                                                                                                  >  3:691225/1‑73 (MQ=255)
                                                                                                     gTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGa                                                              >  3:334417/1‑123 (MQ=255)
                                                                                                     gTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGa                                                              <  4:334417/123‑1 (MQ=255)
                                                                                                         cTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCttt                                          >  2:326369/1‑139 (MQ=255)
                                                                                                          tGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGc                               <  1:701960/149‑1 (MQ=255)
                                                                                                              agcagTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAAt                           >  2:1465665/1‑149 (MQ=255)
                                                                                                               gcagTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTaca                                               <  2:372582/128‑1 (MQ=255)
                                                                                                               gcagTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTaca                                               >  1:372582/1‑128 (MQ=255)
                                                                                                                cagTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGcc                              <  4:329117/144‑1 (MQ=255)
                                                                                                                cagTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGcc                              >  3:329117/1‑144 (MQ=255)
                                                                                                                 agTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGAt                                                             <  4:336810/112‑1 (MQ=255)
                                                                                                                    ggCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGAACATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCttt                     <  4:852454/149‑1 (MQ=255)
                                                                                                                     gCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGcc                                                       <  2:868607/114‑1 (MQ=255)
                                                                                                                     gCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGcc                                                       >  1:868607/1‑114 (MQ=255)
                                                                                                                      cTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTtgt                   <  4:473666/149‑1 (MQ=255)
                                                                                                                        ggAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTtgt                   <  3:1151118/147‑1 (MQ=255)
                                                                                                                        ggAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTtgt                   >  4:1151118/1‑147 (MQ=255)
                                                                                                                           gCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGtgctgc              >  1:21366/1‑149 (MQ=255)
                                                                                                                            cAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCACTGCTGCTCGATATCCAGCCTGGCGATGAAGt                                                                <  4:197025/98‑1 (MQ=255)
                                                                                                                            cAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCACTGCTGCTCGATATCCAGCCTGGCGATGAAGt                                                                >  3:197025/1‑98 (MQ=255)
                                                                                                                             aaCGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGAt                                                             >  1:359415/1‑100 (MQ=255)
                                                                                                                             aaCGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGAt                                                             <  2:359415/100‑1 (MQ=255)
                                                                                                                             aaCGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCttt                     <  1:291815/140‑1 (MQ=255)
                                                                                                                             aaCGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCttt                     >  2:291815/1‑140 (MQ=255)
                                                                                                                             aaCGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGt            >  2:406671/1‑149 (MQ=255)
                                                                                                                              aCGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTTGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGAt                                                             <  2:1249744/99‑1 (MQ=255)
                                                                                                                              aCGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGAt                                                             >  1:1249744/1‑99 (MQ=255)
                                                                                                                                 ttttGGCAGCGCCAAAGTGTTACTGCCGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACAACTTTGTCTCCACCGCCAATGCCTTtgt                   <  3:306491/138‑1 (MQ=255)
                                                                                                                                 ttttGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGgcg        >  1:6515/1‑149 (MQ=255)
                                                                                                                                       cAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCaaaaa  >  1:367777/1‑149 (MQ=255)
                                                                                                                                         |                                                                                                                                                  
CATGATTCCATTTAACGCACCGCCGGTGGTGGGAACCGAACTCGACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAA  >  NC_000913/3975145‑3975428

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