Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,173,692 G→C 5.3% pseudogene (587/914 nt) gatC ← pseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,173,6920GC 5.3% 265.9 / 5.0 75pseudogene (587/914 nt)gatCpseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC
Reads supporting (aligned to +/- strand):  ref base G (41/30);  new base C (2/2);  total (43/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.05e-01

TTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATACCTTCCAGCTCAAAGAAATTTCGGGTATCGCGGGCAAACCAGTCGCCGAGCTTATAAACAAACGCCGCGTGAATTACCACACC  >  NC_000913/2173547‑2173834
                                                                                                                                                 |                                                                                                                                              
tttGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGCGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATATTCTCGATGATAGCAt                                                                                                                                             >  3:4018/1‑149 (MQ=255)
 ttGACATCGTAACCCGCGAGGATGCCGATAATCAGTCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATc                                                                                                                                            <  4:450627/149‑1 (MQ=255)
   gACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAACATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGa                                                                                                                                          >  4:176148/1‑149 (MQ=255)
    aCATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGAc                                                                                                                                         <  3:246996/149‑1 (MQ=255)
      aTCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCa                                                                                                                                       >  1:188481/1‑149 (MQ=255)
      aTCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCa                                                                                                                                       >  3:416516/1‑149 (MQ=255)
        cGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGc                                                                                                                                     <  2:149420/149‑1 (MQ=255)
                      aTGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGccc                                                                                                                       <  4:416516/149‑1 (MQ=255)
                          cGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATACCATCGAc                                                                                                                                         <  1:58632/127‑1 (MQ=255)
                          cGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATACCATCGAc                                                                                                                                         >  4:113316/1‑127 (MQ=255)
                          cGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATACCATCGAc                                                                                                                                         <  3:113316/127‑1 (MQ=255)
                          cGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATACCATCGAc                                                                                                                                         >  2:58632/1‑127 (MQ=255)
                                   gCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGAc                                                                                                          <  4:228872/149‑1 (MQ=255)
                                    ccccATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCAtata                                                                                                                 >  4:622031/1‑141 (MQ=255)
                                    ccccATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCAtata                                                                                                                 <  3:622031/141‑1 (MQ=255)
                                     cccATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGt                                                                                                        >  3:226995/1‑149 (MQ=255)
                                      ccATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTa                                                                                                       >  3:682395/1‑149 (MQ=255)
                                      ccATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTa                                                                                                       <  3:440045/149‑1 (MQ=255)
                                      ccATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTa                                                                                                       >  3:421760/1‑149 (MQ=255)
                                       cATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGTAATCGGCCCCATATACGCCGACGTAc                                                                                                      >  2:350406/1‑149 (MQ=255)
                                             aaaaCCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAGTTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCg                                                                                                >  4:47231/1‑149 (MQ=255)
                                                  ccACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTCTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATa                                                                                           >  3:607461/1‑149 (MQ=255)
                                                    aCGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGc                                                                                         >  4:687105/1‑149 (MQ=255)
                                                         gACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATAc                                                                                    >  4:348803/1‑149 (MQ=255)
                                                           cAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATACCt                                                                                  <  4:226995/149‑1 (MQ=255)
                                                                  gCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATACCTTCCAGCt                                                                           <  3:71767/149‑1 (MQ=255)
                                                                          gACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATACCTTCCAGCTCAAAGaaa                                                                   >  1:81382/1‑149 (MQ=255)
                                                                          gACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATACCTTCCAGCTCAAAGaaa                                                                   <  3:47231/149‑1 (MQ=255)
                                                                             cAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATACCTTCCAGCTCAAAGAAAttt                                                                >  2:648121/1‑149 (MQ=255)
                                                                              aaaTTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATACCTTCCAGCTCAAAGAAATTTc                                                               >  2:622040/1‑149 (MQ=255)
                                                                               aaTTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATACCTTCCAGCTCAAAGAAATTTCg                                                              <  2:361423/149‑1 (MQ=255)
                                                                               aaTTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATACCTTCCAGCTCAAAGAAATTTCg                                                              <  2:317734/149‑1 (MQ=255)
                                                                                aTTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATACCTTCCAGCTCAAAGAAATTTCgg                                                             >  2:131465/1‑149 (MQ=255)
                                                                                      gcTGAATATCGTCGGGGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATACCTTCCAGCTCAAAGAAATTTCGGGTATcg                                                       >  4:510393/1‑149 (MQ=255)
                                                                                        tGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATACCTTCCAGCTCAAAGAAATTTCGGGTATcgcg                                                     >  4:312061/1‑149 (MQ=255)
                                                                                         gAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATACCTTCCAGCTCAAAGAAATTTCGGGTATCGCgg                                                    <  2:673737/149‑1 (MQ=255)
                                                                                          aaTATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATACCTTCCAGCTCAAAGAAATTTCGGGTATCGCggg                                                   >  2:434958/1‑149 (MQ=255)
                                                                                            tatCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATACCTTCCAGCTCAAAGAAATTTCGGGTATCGCGGGCa                                                 >  4:388439/1‑149 (MQ=255)
                                                                                            tatCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATACCTTCCAGCTCAAAGAAATTTCGGGTATCGCGGGCa                                                 >  2:352198/1‑149 (MQ=255)
                                                                                              tcgtcgGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATa                                                                                           >  1:459193/1‑105 (MQ=255)
                                                                                              tcgtcgGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATa                                                                                           <  2:459193/105‑1 (MQ=255)
                                                                                              tcgtcgGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATACCTTCCAGCTCAAAGAAATTTCgg                                                             >  4:590570/1‑135 (MQ=255)
                                                                                              tcgtcgGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATACCTTCCAGCTCAAAGAAATTTCgg                                                             <  3:590570/135‑1 (MQ=255)
                                                                                              tcgtcgGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATACCTTCCAGCTCAAAGAAATTTCGGGTATCGCGGGCaaa                                               >  1:128125/1‑149 (MQ=255)
                                                                                               cgtcgGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATACCTTCCAGCTCAAAGAAATTTCGGGTATCGCGGGCAAAc                                              <  1:689645/149‑1 (MQ=255)
                                                                                                 tcgGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATACCTTCCAGCTCAAAGAAATTTCGGGTATCGCGGGCAAACCa                                            <  2:81382/149‑1 (MQ=255)
                                                                                                 tcgGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATACCTTCCAGCTCAAAGAAATTTCGGGTATCGCGGGCAAACCa                                            <  3:687105/149‑1 (MQ=255)
                                                                                                    gcgcTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATACCTTCCAGCTCAAAGAAATTTCGGGTATCGCGGGCAAACCAGTc                                         <  2:188481/149‑1 (MQ=255)
                                                                                                      gcTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATACCTTCCAGCTCAAAGAAATTTCGGGTATCGCGGGCAAACCAGTCGc                                       <  2:49987/149‑1 (MQ=255)
                                                                                                       cTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATACCTTCCAGCTCAAAGAAATTTCGGGTATCGCGGGCaaa                                               >  2:203166/1‑140 (MQ=255)
                                                                                                       cTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATACCTTCCAGCTCAAAGAAATTTCGGGTATCGCGGGCaaa                                               <  1:203166/140‑1 (MQ=255)
                                                                                                       cTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATACCTTCCAGCTCAAAGAAATTTCGGGTATCGCGGGCAAACCAGTCGcc                                      <  2:129484/149‑1 (MQ=255)
                                                                                                        tAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATACCTTCCAGCTCAAAGAAATTTCGGGTATCGCGGGCAAACCAGTCGCCg                                     >  3:500342/1‑149 (MQ=255)
                                                                                                         aaaTTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATACCTTCCAGCTCAAAGAAATTTCGGGTATCGCGGGCAAACCAGTCGCCGa                                    >  1:571733/1‑149 (MQ=255)
                                                                                                           aTTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATACCTTCCAGCTCAAAGAAATTTCGGGTATCGCGGGCAAACCAGTCGCCGGGc                                  <  4:421760/149‑1 (MQ=255)
                                                                                                           aTTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATACCTTCCAGCTCAAAGAAATTTCGGGTATCGCGGGCAAACCAGTCGCCGAGc                                  <  4:682395/149‑1 (MQ=255)
                                                                                                               aaTTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATACCTTCCAGCTCAAAGAAATTTCGGGTATCGCGGGCAAACCAGTCGCCGAGCTtat                              <  4:382480/149‑1 (MQ=255)
                                                                                                                 ttCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATACCTTCCAGCTCAAAGAAATTTCGGGTATCGCGGGCAAACCAGTCGCCGAGCTTATaa                            >  4:479741/1‑149 (MQ=255)
                                                                                                                     gTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATACCTTCCAGCTCaaa                                                                       <  1:90577/102‑1 (MQ=255)
                                                                                                                     gTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATACCTTCCAGCTCaaa                                                                       >  2:90577/1‑102 (MQ=255)
                                                                                                                     gTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATACCTTCCAGCTCAAAGAAATTTCGGGTATCGCGGGCAAACCAGTCGCCGAGCTTATaaacaa                        >  4:521965/1‑149 (MQ=255)
                                                                                                                      ttAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATACCTTCCAGCTCAAAGAAATTTCGGGTATCGCGGGCAAACCAGTCGCCGAGCTTATaaacaaa                       <  3:348803/149‑1 (MQ=255)
                                                                                                                             cTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATACCTTCCAGCTCAAAGAAATTTCGGGTATCGCGGGCAAACCAGTCGCCGAGCTTATAAACAAACGCcgcg                >  4:361486/1‑149 (MQ=255)
                                                                                                                             cTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATACCTTCCAGCTCAAAGAAATCTCGGGTATCGCGGGCAAACCAGTCGCCGAGCGTATAAACAAACGCcgcg                >  1:247801/1‑149 (MQ=255)
                                                                                                                              tGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGAGCGAAATAGCAATACCTTCCAGCTCAAAGAAATTTCGGGTATCGCGGGCAAACCAGTCGCCGAGCTTATAAACAAACGCCGCGt               <  3:479741/149‑1 (MQ=255)
                                                                                                                               gggATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATACCTTCCAGCTCAAAGAAATTTCGGGTATCGCGGGCAAACCAGTCGCCGAGCTTATAAACAAACGCCGCGTg              >  2:314399/1‑149 (MQ=255)
                                                                                                                                  aTTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATACCTTCCAGCTCAAAGAAATTTCGGGTATCGCGGGCAAACCAGTCGCCGAGCTTATAAACAAACGCCGCGTGAAt           >  2:621225/1‑149 (MQ=255)
                                                                                                                                   ttttCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATACCTTCCAGCTCAAAGAAATTTCGGGTATCGCGGGCAAACCAGTCGCCGAGCTTATAAACAAACGCCGCGTGAAtt          >  3:558071/1‑149 (MQ=255)
                                                                                                                                     ttCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATACCTTCCAGCTCAAAGAAATTTCGGGTATCGCGGGCAAACCAGTCGCCGAGCTTATAAACAAACGCCGCGTGAATTAc        >  1:304651/1‑149 (MQ=255)
                                                                                                                                     ttCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATACCTTCCAGCTCAAAGAAATTTCGGGTATCGCGGGCAAACCAGTCGCCGAGCTTATAAACAAACGCCGCGTGAATTAc        <  1:648121/149‑1 (MQ=255)
                                                                                                                                         cGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATACCTTCCAGCGCAAAGAAATTTAGGGTATCGCGGGCAAACCAGTCGCCGAGCTTATAAACAAACGCCGCGTGAATTACcaca    >  2:493933/1‑149 (MQ=255)
                                                                                                                                          gatgatAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATACCTTCCAGCTCAAAGAAATTTCg                                                              <  2:70550/90‑1 (MQ=255)
                                                                                                                                          gatgatAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATACCTTCCAGCTCAAAGAAATTTCg                                                              >  1:70550/1‑90 (MQ=255)
                                                                                                                                          gataaaaGCATCGACCAGCACCGCAAACGGCCCCATATACGCCGACGTACCGTGCGGAATAACAATACCTTCCAGCTCAAACAAATTTCGGGTATCGCGGGCAAACCAGTCGCCGAGCTTATAAACAAACGCCGCGTGAATTACCacac   >  4:589446/7‑149 (MQ=255)
                                                                                                                                           atgatAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATACCTTCCAGCTCAAAGAAATTTCGGGTATCGCGGGCAAACCAGTCGCCGAGCTTATAAACAAACGCCGCGTGAATTACCACAcc  >  3:188222/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
TTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGCATCGACCAGCACCGCAATCGGCCCCATATACGCCGACGTACCGTGCGGAATAGCAATACCTTCCAGCTCAAAGAAATTTCGGGTATCGCGGGCAAACCAGTCGCCGAGCTTATAAACAAACGCCGCGTGAATTACCACACC  >  NC_000913/2173547‑2173834

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