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

CTGCAACAGCAAGCCAGCGCCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGATCTTGCGTGGA  >  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)
   cAACAGCAAGCCAGCGCCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑t                                                                                                                                                <  2:241876/149‑1 (MQ=255)
         cAAGCCAGCGCCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGaa                                                                                                                                          <  2:323142/149‑1 (MQ=255)
           aGCCAGCGCCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAAtt                                                                                                                                        >  1:853159/1‑149 (MQ=255)
              cAGCGCCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCAt                                                                                                                                     >  1:859586/1‑149 (MQ=255)
              cAGCGCCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCACTGC‑TGCTGAATTCAt                                                                                                                                     <  4:621811/149‑1 (MQ=255)
               aGCGCCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATc                                                                                                                                    <  3:248423/149‑1 (MQ=255)
                gcgcCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATcc                                                                                                                                   >  1:394620/1‑149 (MQ=255)
                 cgcCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATccc                                                                                                                                  <  4:527674/149‑1 (MQ=255)
                  gcCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATcccc                                                                                                                                 <  1:19555/149‑1 (MQ=255)
                           ggCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGc                                                                                                                        <  3:367393/149‑1 (MQ=255)
                              cTGGATGTGTGCGAAATGGTATTGCTGGAGGTCCAGGTGCAGCTTGAGCAGATGCTGCTGGAAGTCCAAAGTAACAGATAGCGCAGATAACGACTATCGATGCCAATGAATTCAGCAGCAG‑TGCGGAATTCATCCCCATGCGCGGCttt                                                                                                                     <  4:819916/149‑1 (MQ=255)
                                 gAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGc                                                                                                                  <  4:214417/149‑1 (MQ=255)
                                    gagTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAAt                                                                                                               <  2:394620/149‑1 (MQ=255)
                                     agTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAAtt                                                                                                              <  3:118483/149‑1 (MQ=255)
                                     agTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAAtt                                                                                                              <  3:237141/149‑1 (MQ=255)
                                         cGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCa                                                                                                          >  2:262540/1‑149 (MQ=255)
                                          gAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCAt                                                                                                         >  2:792675/1‑149 (MQ=255)
                                            aaTGGCATGGCTGGAAGCCCAGGAGAAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCt                                                                                                       <  3:576957/149‑1 (MQ=255)
                                                gcatggcTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCaa                                                                                                   >  2:490235/1‑149 (MQ=255)
                                                  atggcTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGc                                                                                                                        >  1:231990/1‑126 (MQ=255)
                                                  atggcTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGc                                                                                                                        <  2:231990/126‑1 (MQ=255)
                                                  atggcTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAAc                                                                                                 >  3:324675/1‑149 (MQ=255)
                                                               cAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGt                                                                                    >  4:764113/1‑149 (MQ=255)
                                                               cAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGt                                                                                    >  1:666521/1‑149 (MQ=255)
                                                               cAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGt                                                                                    >  1:673934/1‑149 (MQ=255)
                                                                 ggAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCt                                                                                  >  1:883009/1‑149 (MQ=255)
                                                                 ggAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATACTGCCGCGTCt                                                                                  <  3:692545/149‑1 (MQ=255)
                                                                   agcagcTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGt                                                                                <  2:883009/149‑1 (MQ=255)
                                                                     cagcTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGttt                                                                              >  2:223151/1‑149 (MQ=255)
                                                                     cagcTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCACTGC‑TGCTGAATTCAttg                                                                                                                                   >  2:40970/1‑94 (MQ=38)
                                                                     cagcTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCA‑CTGCTGCTGAATTCAttg                                                                                                                                   <  1:40970/96‑3 (MQ=38)
                                                                      agcTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTtc                                                                             <  2:260081/149‑1 (MQ=255)
                                                                      agcTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTtc                                                                             <  3:41594/149‑1 (MQ=255)
                                                                       gcTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTtcg                                                                            >  4:418874/1‑149 (MQ=255)
                                                                        cTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTc                                                                                   >  2:766497/1‑141 (MQ=255)
                                                                        cTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTc                                                                                   <  1:766497/141‑1 (MQ=255)
                                                                        cTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTtcgt                                                                           >  3:426547/1‑149 (MQ=255)
                                                                         ttGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTtcgtc                                                                          <  2:859586/149‑1 (MQ=255)
                                                                         ttGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTtcgtc                                                                          <  3:518569/149‑1 (MQ=255)
                                                                           gAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTcgtcgt                                                                        <  4:415752/149‑1 (MQ=255)
                                                                             gCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCa                                                                      >  1:290828/1‑149 (MQ=255)
                                                                                gATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGt                                                                   >  1:860715/1‑149 (MQ=255)
                                                                                    ctgctgGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGaaa                                                               <  3:683796/149‑1 (MQ=255)
                                                                                    ctgctgGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGaaa                                                               <  4:696014/149‑1 (MQ=255)
                                                                                     tgctgGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAAt                                                              >  1:191681/1‑149 (MQ=255)
                                                                                        tgGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCa                                                           <  4:39678/149‑1 (MQ=255)
                                                                                         ggAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCAc                                                          <  2:666521/149‑1 (MQ=255)
                                                                                         ggAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCAc                                                          <  3:398050/149‑1 (MQ=255)
                                                                                         ggAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCAc                                                          <  1:31774/149‑1 (MQ=255)
                                                                                           aaGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCAccc                                                        >  4:219709/1‑149 (MQ=255)
                                                                                               ccAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTa                                                    >  2:383915/1‑149 (MQ=255)
                                                                                                cAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTAc                                                   >  3:335377/1‑149 (MQ=255)
                                                                                                cAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCAGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTAc                                                   >  2:683716/1‑149 (MQ=255)
                                                                                                cAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCACTGC‑TGCTGAATTCAtt                                                                                                                                    >  3:781818/1‑67 (MQ=21)
                                                                                                cAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCACTGC‑TGCTGAATTCAtt                                                                                                                                    <  4:781818/68‑2 (MQ=21)
                                                                                                 aaaGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACt                                                  >  3:61509/1‑149 (MQ=255)
                                                                                                      aaCTTATGGCTCAGAGAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGTTGAAATTCACCCCGTACTGGCAg                                             >  2:140001/1‑149 (MQ=255)
                                                                                                       aCTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGc                                            <  3:686632/149‑1 (MQ=255)
                                                                                                         tGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCgg                                          >  1:894457/1‑149 (MQ=255)
                                                                                                           aTGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCt                                        >  3:263107/1‑149 (MQ=255)
                                                                                                            tGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTg                                       <  1:241533/149‑1 (MQ=255)
                                                                                                             ggCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTgg                                      >  1:13746/1‑149 (MQ=255)
                                                                                                              gcgcAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCaaa                                                                                                  <  4:182340/88‑1 (MQ=255)
                                                                                                              gcgcAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCaaa                                                                                                  >  3:182340/1‑88 (MQ=255)
                                                                                                               cgcAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAg                                    <  4:206719/149‑1 (MQ=255)
                                                                                                                 cAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGtt                                  <  3:354439/149‑1 (MQ=255)
                                                                                                                   gATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGc                                <  2:290828/149‑1 (MQ=255)
                                                                                                                     tAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGccc                              >  1:159947/1‑149 (MQ=255)
                                                                                                                      aaCGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGc                                <  4:33669/146‑1 (MQ=255)
                                                                                                                      aaCGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGc                                >  3:33669/1‑146 (MQ=255)
                                                                                                                       acgacgACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGa                            <  1:262540/149‑1 (MQ=255)
                                                                                                                         gacgacCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACg                          >  3:672691/1‑149 (MQ=255)
                                                                                                                                   gCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGt                >  1:143534/1‑149 (MQ=255)
                                                                                                                                     cAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGCCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCg              >  1:305436/1‑149 (MQ=255)
                                                                                                                                     cAATGAATTCAGCAGCAG‑TGCTGAATGCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCTTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCg              >  2:420061/1‑149 (MQ=255)
                                                                                                                                      aaTGAATTCAGCAGCAG‑TGCTGAATTCGTCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGa             >  3:690059/1‑149 (MQ=255)
                                                                                                                                       aTGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGAt            <  4:314386/149‑1 (MQ=255)
                                                                                                                                       aTGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGAt            >  4:163484/1‑149 (MQ=255)
                                                                                                                                          aaTTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGATCtt         >  3:594302/1‑149 (MQ=255)
                                                                                                                                           aTTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGATCTTg        >  3:226012/1‑149 (MQ=255)
                                                                                                                                              cagcagcag‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCg                                           <  3:490794/111‑1 (MQ=255)
                                                                                                                                              cagcagcag‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCg                                           >  4:490794/1‑111 (MQ=255)
                                                                                                                                              cagcagcag‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGATCTTGCGt     >  4:496158/1‑149 (MQ=255)
                                                                                                                                              cagcagcag‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGATCTTGCGt     >  2:364746/1‑149 (MQ=255)
                                                                                                                                               agcagcag‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGATCTTGCGTg    <  2:71243/149‑1 (MQ=255)
                                                                                                                                                 cagcag‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGATCTTGCGTGGa  <  4:226012/149‑1 (MQ=255)
                                                                                                                                                 cagcag‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGATCTTGCGTGGa  <  3:163484/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CTGCAACAGCAAGCCAGCGCCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAG‑TGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGATCTTGCGTGGA  >  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: ‑