Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,131,524 A→G 6.4% S295S (AGT→AGC) ‡ yghQ ← putative inner membrane polysaccharide flippase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,131,5240AG 6.4% 363.3 / 10.7 118S295S (AGT→AGC) ‡yghQputative inner membrane polysaccharide flippase
Reads supporting (aligned to +/- strand):  ref base A (60/50);  new base G (3/5);  total (63/55)
Fisher's exact test for biased strand distribution p-value = 4.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.07e-01

Only 100 of 125 total aligned reads displayed.

AAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAAT  >  NC_000913/3131377‑3131670
                                                                                                                                                   |                                                                                                                                                  
aaaCACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAAc                                                                                                                                                   >  1:133877/1‑149 (MQ=255)
  acacCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACtt                                                                                                                                                 <  3:518936/149‑1 (MQ=255)
     ccAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATCCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTa                                                                                                                                              <  3:527962/149‑1 (MQ=255)
      cAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAGTTTTTAc                                                                                                                                             >  3:876128/1‑149 (MQ=255)
       aGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTa                                                                                                                                              <  3:863672/147‑1 (MQ=255)
       aGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTa                                                                                                                                              >  4:863672/1‑147 (MQ=255)
       aGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTAcc                                                                                                                                            >  4:576788/1‑149 (MQ=255)
         cGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAGTTTTTACCCa                                                                                                                                          <  1:678743/149‑1 (MQ=255)
           aaaTGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAGTTTTTACCCAga                                                                                                                                        <  2:796451/149‑2 (MQ=255)
           aaaTGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAGTTTTTACCCAga                                                                                                                                        <  4:876128/149‑2 (MQ=255)
            aaTGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAgaa                                                                                                                                       <  1:1015990/149‑3 (MQ=255)
              tGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCaaa                                                                                                                                     >  2:932593/1‑149 (MQ=255)
              tGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCaaa                                                                                                                                     <  1:1049140/149‑1 (MQ=255)
                   ggCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGc                                                                                                                                <  1:541170/149‑1 (MQ=255)
                   ggCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGc                                                                                                                                <  3:31989/149‑1 (MQ=255)
                      ttGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTgg                                                                                                                             <  1:750909/149‑1 (MQ=255)
                          cgacgaTCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGtt                                                                                                                         >  3:519187/1‑149 (MQ=255)
                          cgacgaTCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGtt                                                                                                                         >  3:591811/1‑149 (MQ=255)
                          cgacgaTCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGtt                                                                                                                         >  4:510135/1‑149 (MQ=255)
                          cgacgaTCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGAAAAAACTTTTACCCAGCAAACCCGCTGGCGtt                                                                                                                         >  4:647211/1‑149 (MQ=255)
                           gacgaTCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTc                                                                                                                        >  2:826/1‑149 (MQ=255)
                             cgaTCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCcg                                                                                                                      <  3:146175/149‑1 (MQ=255)
                               aTCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCcggc                                                                                                                    >  3:271303/1‑149 (MQ=255)
                                tCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCcggcg                                                                                                                   >  3:52884/1‑149 (MQ=255)
                                tCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCcggcg                                                                                                                   >  4:571500/1‑149 (MQ=255)
                                  aaCACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCggcggc                                                                                                                 <  1:932593/149‑1 (MQ=255)
                                  aaCACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCggcggc                                                                                                                 >  2:784883/1‑149 (MQ=255)
                                  aaCACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCggcggc                                                                                                                 >  1:1036873/1‑149 (MQ=255)
                                   acacCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACtttt                                                                                                                                               >  2:141905/1‑118 (MQ=255)
                                   acacCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACtttt                                                                                                                                               <  1:141905/118‑1 (MQ=255)
                                           aGCGCCACCAGAATCCCGTTTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCAGCGGCGTCgaagaa                                                                                                        <  3:79053/149‑1 (MQ=255)
                                            gcgcCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAAt                                                                                                       <  1:226241/149‑1 (MQ=255)
                                            gcgcCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAAt                                                                                                       <  3:642949/149‑1 (MQ=255)
                                               ccaccaGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTc                                                                                                    <  3:1021418/149‑1 (MQ=255)
                                                caccaGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAGTTTTTACCCAgaggt                                                                                                                                     >  4:125754/1‑111 (MQ=255)
                                                caccaGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAGTTTTTACCCAgaggt                                                                                                                                     <  3:125754/115‑5 (MQ=255)
                                                caccaGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCa                                                                                                   <  3:510135/149‑1 (MQ=255)
                                                 accaGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCAt                                                                                                  >  3:802674/1‑149 (MQ=255)
                                                 accaGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCAt                                                                                                  >  1:78191/1‑149 (MQ=255)
                                                  ccaGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATg                                                                                                 >  1:494105/1‑149 (MQ=255)
                                                        tACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGAtt                                                                                           >  3:107895/1‑149 (MQ=255)
                                                          ccGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGAtttt                                                                                         <  1:784883/149‑1 (MQ=255)
                                                          ccGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGAtttt                                                                                         >  2:941544/1‑149 (MQ=255)
                                                             aTTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTaaa                                                                                      <  1:440648/149‑1 (MQ=255)
                                                              taccgcccgcCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAAt                                                                                     <  3:647211/147‑1 (MQ=255)
                                                                    ccgccAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCa                                                                                                                                          <  3:116612/90‑1 (MQ=255)
                                                                    ccgccAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCa                                                                                                                                          >  4:116612/1‑90 (MQ=255)
                                                                    ccgccAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCcg                                                                               >  4:723425/1‑149 (MQ=255)
                                                                     cgccAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTcc                                                                                                                       <  2:978533/108‑1 (MQ=255)
                                                                     cgccAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTcc                                                                                                                       >  1:978533/1‑108 (MQ=255)
                                                                      gccAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCcggc                                                                             <  1:41993/149‑1 (MQ=255)
                                                                      gccAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCcggc                                                                             <  1:279795/149‑1 (MQ=255)
                                                                              aCCGAACTACACACCTAACTACCACTGTCCGGTGGTGCGCGGATCTAAGCCCACGACCTCTGCGTAAAAACTTGTACTCGCCAAACCCGCTGGCGCGCCGGCGGCGTCGAAGAATGTCATGGCGACTTTAAATGACCCGGCGGCAGCgg                                                                     <  4:55368/149‑1 (MQ=255)
                                                                               ccGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTAGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCggg                                                                    <  3:768611/149‑1 (MQ=255)
                                                                                cGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGt                                                                   >  3:330990/1‑149 (MQ=255)
                                                                                  gACTTCACATCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCggg                                                                    >  4:932087/1‑146 (MQ=255)
                                                                                  gACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCTCCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCggg                                                                    <  3:932087/146‑1 (MQ=255)
                                                                                          acCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGa                                                         >  3:179025/1‑149 (MQ=255)
                                                                                          acCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGa                                                         >  1:534713/1‑149 (MQ=255)
                                                                                           ccTAACAACCACGTTCTGGTGGTGCGCGGATCTAGACGCAGCACCTCTGGGTAAAAGCGTTTACTCAGCAAACCCGCTGGCGTTCCGGCAGCGTCGAAGAATGTCATGGCGATTTTAAATGACCCGGCGGCAGCGGGTCCTaatgcgtg                                                        <  2:802310/149‑7 (MQ=255)
                                                                                            ctcacaacCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATc                                                       >  4:471690/4‑149 (MQ=255)
                                                                                            cTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATc                                                       >  1:1003737/1‑149 (MQ=255)
                                                                                            cTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATc                                                       <  1:682023/149‑1 (MQ=255)
                                                                                            cTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATc                                                       <  4:802674/149‑1 (MQ=255)
                                                                                             tAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATcc                                                      >  2:463935/1‑149 (MQ=255)
                                                                                              aacaacCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATccc                                                     >  3:56964/1‑149 (MQ=255)
                                                                                               acaacCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATcccc                                                    <  4:271303/149‑1 (MQ=255)
                                                                                                 aacCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCccac                                                  >  1:608880/1‑149 (MQ=255)
                                                                                                  acCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCccacc                                                 <  3:821333/149‑1 (MQ=255)
                                                                                                  acCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCccacc                                                 <  2:1058455/149‑1 (MQ=255)
                                                                                                       ggTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTCACACGATCCCCACCAacac                                            <  4:474623/149‑1 (MQ=255)
                                                                                                        gTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACt                                           >  3:400769/1‑149 (MQ=255)
                                                                                                         tCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAGAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACtg                                          <  4:227575/149‑1 (MQ=255)
                                                                                                          cTGGTGGTGCGCGGATCCAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACtgt                                         <  4:667129/149‑1 (MQ=255)
                                                                                                                    gcgGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAg                               >  2:1026608/1‑149 (MQ=255)
                                                                                                                     cgGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAaca                                                            >  1:405973/1‑119 (MQ=255)
                                                                                                                     cgGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGt                              <  4:124936/149‑1 (MQ=255)
                                                                                                                       gATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTa                            >  1:147352/1‑149 (MQ=255)
                                                                                                                        aTCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTAc                           >  2:24333/1‑149 (MQ=255)
                                                                                                                             aaCGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGcc                      >  3:154245/1‑149 (MQ=255)
                                                                                                                             aaCGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGAGcc                      >  4:698895/1‑149 (MQ=255)
                                                                                                                                    cACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGGCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGa               >  2:535732/1‑149 (MQ=255)
                                                                                                                                     tttgctGGGTAAAAGTTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGAt              >  2:192903/5‑149 (MQ=255)
                                                                                                                                     aCCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGAt              >  4:217329/1‑149 (MQ=255)
                                                                                                                                      ccTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATg             >  1:674183/1‑149 (MQ=255)
                                                                                                                                       ctctGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATgg            >  4:450761/1‑149 (MQ=255)
                                                                                                                                        tctGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCTGATGGa           >  3:633696/1‑149 (MQ=255)
                                                                                                                                         ctGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGTGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAg          <  2:46320/149‑1 (MQ=255)
                                                                                                                                           gggTAAAAACTTTTACCCAGCAACCCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTTCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTg        <  2:312387/149‑1 (MQ=255)
                                                                                                                                            ggTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTgg       <  3:217329/149‑1 (MQ=255)
                                                                                                                                            ggTAAAAACTCTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTgg       >  2:945594/1‑149 (MQ=255)
                                                                                                                                             gTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGt         >  1:964498/1‑146 (MQ=255)
                                                                                                                                              tAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACtgt                                         >  4:389108/1‑113 (MQ=255)
                                                                                                                                              tAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGc     <  3:576788/149‑1 (MQ=255)
                                                                                                                                               aaaaaCTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCa    >  2:417700/1‑149 (MQ=255)
                                                                                                                                                aaaaCTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCaa   >  2:873186/1‑149 (MQ=255)
                                                                                                                                                 aaaCTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCggg                                                                    <  1:856535/83‑1 (MQ=255)
                                                                                                                                                 aaaCTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAAt  >  4:117559/1‑149 (MQ=255)
                                                                                                                                                  aaCTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAacac                                            <  4:906692/106‑1 (MQ=255)
                                                                                                                                                  aaCTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAacac                                            >  3:906692/1‑106 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAAT  >  NC_000913/3131377‑3131670

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