Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,483,242 G→C 5.0% S14T (AGC→ACC)  yheT → UPF0017 family putative hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,483,2420GC 5.0% 278.7 / 5.3 80S14T (AGC→ACC) yheTUPF0017 family putative hydrolase
Reads supporting (aligned to +/- strand):  ref base G (39/36);  new base C (2/2);  total (41/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.38e-01

CTGCAACAGCAAGCCAGCGCCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGATCTTGCGTGGA  >  NC_000913/3483095‑3483388
                                                                                                                                                   |                                                                                                                                                   
ctgcAACAGCAAGCCAGCGCCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTcagcagc                                                                                                                                                    >  4:237141/1‑149 (MQ=255)
 tgcAACAGCAAGCCAGCGCCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTcagcagca                                                                                                                                                   >  2:514938/1‑149 (MQ=255)
 tgcAACAGCAAGCCAGCGCCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTcagcagca                                                                                                                                                   <  3:623832/149‑1 (MQ=255)
   cAACAGCAAGCCAGCGCCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGt                                                                                                                                                <  2:241876/149‑1 (MQ=255)
         cAAGCCAGCGCCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGaa                                                                                                                                          <  2:323142/149‑1 (MQ=255)
           aGCCAGCGCCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAAtt                                                                                                                                        >  1:853159/1‑149 (MQ=255)
              cAGCGCCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCAt                                                                                                                                     >  1:859586/1‑149 (MQ=255)
              cAGCGCCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCACTGCTGCTGAATTCAt                                                                                                                                     <  4:621811/149‑1 (MQ=255)
               aGCGCCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATc                                                                                                                                    <  3:248423/149‑1 (MQ=255)
                gcgcCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATcc                                                                                                                                   >  1:394620/1‑149 (MQ=255)
                 cgcCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATccc                                                                                                                                  <  4:527674/149‑1 (MQ=255)
                  gcCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATcccc                                                                                                                                 <  1:19555/149‑1 (MQ=255)
                           ggCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGc                                                                                                                        <  3:367393/149‑1 (MQ=255)
                              cTGGATGTGTGCGAAATGGTATTGCTGGAGGTCCAGGTGCAGCTTGAGCAGATGCTGCTGGAAGTCCAAAGTAACAGATAGCGCAGATAACGACTATCGATGCCAATGAATTCAGCAGCAGTGCGGAATTCATCCCCATGCGCGGCttt                                                                                                                     <  4:819916/149‑1 (MQ=255)
                                 gAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGc                                                                                                                  <  4:214417/149‑1 (MQ=255)
                                    gagTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAAt                                                                                                               <  2:394620/149‑1 (MQ=255)
                                     agTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAAtt                                                                                                              <  3:118483/149‑1 (MQ=255)
                                     agTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAAtt                                                                                                              <  3:237141/149‑1 (MQ=255)
                                         cGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCa                                                                                                          >  2:262540/1‑149 (MQ=255)
                                          gAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCAt                                                                                                         >  2:792675/1‑149 (MQ=255)
                                            aaTGGCATGGCTGGAAGCCCAGGAGAAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCt                                                                                                       <  3:576957/149‑1 (MQ=255)
                                                gcatggcTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCaa                                                                                                   >  2:490235/1‑149 (MQ=255)
                                                  atggcTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGc                                                                                                                        >  1:231990/1‑126 (MQ=255)
                                                  atggcTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGc                                                                                                                        <  2:231990/126‑1 (MQ=255)
                                                  atggcTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAAc                                                                                                 >  3:324675/1‑149 (MQ=255)
                                                               cAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGt                                                                                    >  4:764113/1‑149 (MQ=255)
                                                               cAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGt                                                                                    >  1:666521/1‑149 (MQ=255)
                                                               cAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGt                                                                                    >  1:673934/1‑149 (MQ=255)
                                                                 ggAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCt                                                                                  >  1:883009/1‑149 (MQ=255)
                                                                 ggAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATACTGCCGCGTCt                                                                                  <  3:692545/149‑1 (MQ=255)
                                                                   agcagcTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGt                                                                                <  2:883009/149‑1 (MQ=255)
                                                                     cagcTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGttt                                                                              >  2:223151/1‑149 (MQ=255)
                                                                     cagcTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCACTGCTGCTGAATTCAttg                                                                                                                                   >  2:40970/1‑94 (MQ=38)
                                                                     cagcTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCACTGCTGCTGAATTCAttg                                                                                                                                   <  1:40970/96‑3 (MQ=38)
                                                                      agcTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTtc                                                                             <  2:260081/149‑1 (MQ=255)
                                                                      agcTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTtc                                                                             <  3:41594/149‑1 (MQ=255)
                                                                       gcTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTtcg                                                                            >  4:418874/1‑149 (MQ=255)
                                                                        cTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTc                                                                                   >  2:766497/1‑141 (MQ=255)
                                                                        cTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTc                                                                                   <  1:766497/141‑1 (MQ=255)
                                                                        cTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTtcgt                                                                           >  3:426547/1‑149 (MQ=255)
                                                                         ttGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTtcgtc                                                                          <  2:859586/149‑1 (MQ=255)
                                                                         ttGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTtcgtc                                                                          <  3:518569/149‑1 (MQ=255)
                                                                           gAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTcgtcgt                                                                        <  4:415752/149‑1 (MQ=255)
                                                                             gCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCa                                                                      >  1:290828/1‑149 (MQ=255)
                                                                                gATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGt                                                                   >  1:860715/1‑149 (MQ=255)
                                                                                    ctgctgGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGaaa                                                               <  3:683796/149‑1 (MQ=255)
                                                                                    ctgctgGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGaaa                                                               <  4:696014/149‑1 (MQ=255)
                                                                                     tgctgGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAAt                                                              >  1:191681/1‑149 (MQ=255)
                                                                                        tgGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCa                                                           <  4:39678/149‑1 (MQ=255)
                                                                                         ggAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCAc                                                          <  2:666521/149‑1 (MQ=255)
                                                                                         ggAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCAc                                                          <  3:398050/149‑1 (MQ=255)
                                                                                         ggAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCAc                                                          <  1:31774/149‑1 (MQ=255)
                                                                                           aaGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCAccc                                                        >  4:219709/1‑149 (MQ=255)
                                                                                               ccAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTa                                                    >  2:383915/1‑149 (MQ=255)
                                                                                                cAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTAc                                                   >  3:335377/1‑149 (MQ=255)
                                                                                                cAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCAGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTAc                                                   >  2:683716/1‑149 (MQ=255)
                                                                                                cAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCACTGCTGCTGAATTCAtt                                                                                                                                    >  3:781818/1‑67 (MQ=21)
                                                                                                cAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCACTGCTGCTGAATTCAtt                                                                                                                                    <  4:781818/68‑2 (MQ=21)
                                                                                                 aaaGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACt                                                  >  3:61509/1‑149 (MQ=255)
                                                                                                      aaCTTATGGCTCAGAGAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGTTGAAATTCACCCCGTACTGGCAg                                             >  2:140001/1‑149 (MQ=255)
                                                                                                       aCTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGc                                            <  3:686632/149‑1 (MQ=255)
                                                                                                         tGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCgg                                          >  1:894457/1‑149 (MQ=255)
                                                                                                           aTGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCt                                        >  3:263107/1‑149 (MQ=255)
                                                                                                            tGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTg                                       <  1:241533/149‑1 (MQ=255)
                                                                                                             ggCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTgg                                      >  1:13746/1‑149 (MQ=255)
                                                                                                              gcgcAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCaaa                                                                                                  <  4:182340/88‑1 (MQ=255)
                                                                                                              gcgcAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCaaa                                                                                                  >  3:182340/1‑88 (MQ=255)
                                                                                                               cgcAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAg                                    <  4:206719/149‑1 (MQ=255)
                                                                                                                 cAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGtt                                  <  3:354439/149‑1 (MQ=255)
                                                                                                                   gATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGc                                <  2:290828/149‑1 (MQ=255)
                                                                                                                     tAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGccc                              >  1:159947/1‑149 (MQ=255)
                                                                                                                      aaCGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGc                                <  4:33669/146‑1 (MQ=255)
                                                                                                                      aaCGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGc                                >  3:33669/1‑146 (MQ=255)
                                                                                                                       acgacgACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGa                            <  1:262540/149‑1 (MQ=255)
                                                                                                                         gacgacCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACg                          >  3:672691/1‑149 (MQ=255)
                                                                                                                                   gCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGt                >  1:143534/1‑149 (MQ=255)
                                                                                                                                     cAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGCCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCg              >  1:305436/1‑149 (MQ=255)
                                                                                                                                     cAATGAATTCAGCAGCAGTGCTGAATGCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCTTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCg              >  2:420061/1‑149 (MQ=255)
                                                                                                                                      aaTGAATTCAGCAGCAGTGCTGAATTCGTCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGa             >  3:690059/1‑149 (MQ=255)
                                                                                                                                       aTGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGAt            <  4:314386/149‑1 (MQ=255)
                                                                                                                                       aTGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGAt            >  4:163484/1‑149 (MQ=255)
                                                                                                                                          aaTTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGATCtt         >  3:594302/1‑149 (MQ=255)
                                                                                                                                           aTTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGATCTTg        >  3:226012/1‑149 (MQ=255)
                                                                                                                                              cagcagcagTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCg                                           <  3:490794/111‑1 (MQ=255)
                                                                                                                                              cagcagcagTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCg                                           >  4:490794/1‑111 (MQ=255)
                                                                                                                                              cagcagcagTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGATCTTGCGt     >  4:496158/1‑149 (MQ=255)
                                                                                                                                              cagcagcagTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGATCTTGCGt     >  2:364746/1‑149 (MQ=255)
                                                                                                                                               agcagcagTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGATCTTGCGTg    <  2:71243/149‑1 (MQ=255)
                                                                                                                                                 cagcagTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGATCTTGCGTGGa  <  4:226012/149‑1 (MQ=255)
                                                                                                                                                 cagcagTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGATCTTGCGTGGa  <  3:163484/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CTGCAACAGCAAGCCAGCGCCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGATCTTGCGTGGA  >  NC_000913/3483095‑3483388

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: