Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,906,609 G→C 5.1% intergenic (+26/+34) ygcG → / ← eno TPM domain protein, putative phosphatase/enolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,906,6090GC 5.1% 215.1 / 3.5 78intergenic (+26/+34)ygcG/enoTPM domain protein, putative phosphatase/enolase
Reads supporting (aligned to +/- strand):  ref base G (26/46);  new base C (2/2);  total (29/49)
Fisher's exact test for biased strand distribution p-value = 6.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

AAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATAGAA  >  NC_000913/2906461‑2906757
                                                                                                                                                    |                                                                                                                                                    
aaaaGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTATCTCAACCATAAAAAATGCCAGCCCgg                                                                                                                                                      >  2:52341/1‑149 (MQ=255)
aaaaGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCAGGAGGCTCAGATAGCGGCGTATTTAGCGGCGGGGGCGGATCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCACCCATAAAAAATGCCAGACCgg                                                                                                                                                      >  2:549243/1‑149 (MQ=255)
 aaaGAGCATCTTGGAGAAGTTATCATTCAGACAACAGCGGTTCCGCTTGAGGCTCAGATAGCGGCGGCGTTAGAGGCGGGGGCGGTTCTTCTGGAGGGGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGa                                                                                                                                                     <  1:761605/149‑1 (MQ=255)
          cTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCAt                                                                                                                                            <  1:991751/149‑1 (MQ=255)
            ggtactaGTTATCATTCTGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCTCCGGGCTGGCAttt                                                                                                                                          <  1:233227/143‑1 (MQ=255)
            ggTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCAttt                                                                                                                                          >  2:435231/1‑149 (MQ=255)
            ggTAGAAGTTATCATTCAGACAACAGCGGTTCCGCGGGAGGCTCAGATAGCGGCGGATTTCGCGGCGGGGGCGGTTCTTCTGGAGGGGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGGGGCTGGCAttt                                                                                                                                          <  2:874579/149‑1 (MQ=255)
             gTAGAAGTTATCATTCAGACAACAGCGGGTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCAtttt                                                                                                                                         <  2:986376/149‑1 (MQ=255)
              tcgAAGGTCTCATTCAGACAACAGCCGTTCCGATGGAGGCGCAGATCGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGGGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCTCCGGGCTGGCAttttt                                                                                                                                        <  2:872679/147‑1 (MQ=255)
                gAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGCGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTACCACAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTaa                                                                                                                                      >  1:348333/1‑149 (MQ=255)
                 aaGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTaaa                                                                                                                                     <  2:725915/149‑1 (MQ=255)
                  aggtATCATGCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGCTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAAt                                                                                                                                    <  1:358867/146‑1 (MQ=255)
                      aTCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGTGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCAGGAGGCAGGCATTTTTAAATCAGa                                                                                                                                >  1:1012955/1‑149 (MQ=255)
                      aTCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGa                                                                                                                                >  2:154585/1‑149 (MQ=255)
                         aTTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATaa                                                                                                                             >  1:313174/1‑149 (MQ=255)
                         aTTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCCGCGGGGGCGGTTCTTCAGGAGGTGGCGGCGCATCCCGCCGCTGGTAACTCAACCATAAAAAATGCCCGCCCGGAGGCGGGCATTGTTAAATCAGATaa                                                                                                                             >  1:707941/1‑149 (MQ=255)
                          ttCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGTCATTGTTAAATCAGATaaa                                                                                                                            >  1:954353/1‑149 (MQ=255)
                          ttCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATaaa                                                                                                                            >  1:530543/1‑149 (MQ=255)
                                 tccaGCGGAACCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTTACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAagtcagtc                                                                                                                     >  1:770283/3‑149 (MQ=255)
                                 aacaGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAagtcagtc                                                                                                                     >  1:626114/1‑149 (MQ=255)
                                  acaGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCt                                                                                                                    <  1:149767/149‑1 (MQ=255)
                                    aGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCGTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATATAGTCAGTCttt                                                                                                                  >  1:760298/1‑148 (MQ=255)
                                       ggTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGc                                                                                                               >  1:38211/1‑149 (MQ=255)
                                        gTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGcc                                                                                                              >  1:99285/1‑149 (MQ=255)
                                        gTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGcc                                                                                                              >  1:703208/1‑149 (MQ=255)
                                              cTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGAGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTAAGCCTGGCCt                                                                                                        <  2:745993/149‑1 (MQ=255)
                                                ggAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCttt                                                                                                      <  2:718662/149‑1 (MQ=255)
                                                 gAGGCGCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTg                                                                                                     <  2:262056/149‑1 (MQ=255)
                                                   ggCTCAGTTAGCGGCGGATTTAGCGGCGGCGGCGGGTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGTCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTCGCCTTTGAt                                                                                                   <  2:855141/149‑1 (MQ=255)
                                                   ggCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGAt                                                                                                   <  1:579908/149‑1 (MQ=255)
                                                    gCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGAtc                                                                                                  <  1:370642/149‑1 (MQ=255)
                                                    gCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGAtc                                                                                                  <  1:47746/149‑1 (MQ=255)
                                                     gtctgttAGCGGCGTTTTTAGCGGCGGGGGGGGTTCTTCTGGAGGTGGCGGCGCATGCGGGCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTATATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGTtct                                                                                                 <  1:554840/143‑1 (MQ=255)
                                                       cAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATctct                                                                                               <  2:516593/149‑1 (MQ=255)
                                                       cAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATctct                                                                                               <  1:127109/149‑1 (MQ=255)
                                                       cAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATctct                                                                                               <  2:598809/149‑1 (MQ=255)
                                                        aGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCtt                                                                                              <  1:376917/149‑1 (MQ=255)
                                                         gATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCttt                                                                                             <  2:703208/149‑1 (MQ=255)
                                                         gATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCttt                                                                                             <  2:99285/149‑1 (MQ=255)
                                                          aTAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTa                                                                                            <  1:999419/149‑1 (MQ=255)
                                                          aTAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCATGGCCTTTGATCTCTTTa                                                                                            >  2:832476/1‑149 (MQ=255)
                                                          aTAGCGGCGGATTTAGAGGCGTGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACGATAAAAAATGCCAGCCCGTAGGCTGGCATTTTTAAATCAGATAAAGTCTGTCTTATGCCTGGCCTTTGATCTCTTTa                                                                                            <  2:141542/149‑1 (MQ=255)
                                                              cggcggATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGAc                                                                                        <  2:791779/149‑1 (MQ=255)
                                                                  ggATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGtt                                                                                    <  1:652779/149‑1 (MQ=255)
                                                                  ggATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGtt                                                                                    <  2:776594/149‑1 (MQ=255)
                                                                        aGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCAGTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACgg                                                                              >  1:606522/1‑149 (MQ=255)
                                                                           ggcggGGGCGGATCTTATGAAGGTGGCGGCGCATCCGGCCACTGGTAACTCAACCATAAAAAATGCCAGCCCGTAGTCTGGCATTTTTAAAACAGATAAAGTCATTCTTATGCCTGTCCTTTGATATCTTTACGTCCTTTGTACTGTGc                                                                           >  2:618990/1‑149 (MQ=255)
                                                                            gcggGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGTCCGCTGGTAACTAAACCATAAAAAATTCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTAATGCCTGGCCTTTGATCTCTTTACGACCGTAGTACGGTGCt                                                                          >  2:683991/1‑149 (MQ=255)
                                                                             cggGGGCGGTTCTTCTGGAGGTGGCGGGGCATCCGGCGGCTGGTAAGTCAACGATAAAATATGGCAGCCCGGAGGCTGGCATTTTTAAATCATATAAAGGCAGTGTTATGCCTGTCCTTTGATCTCTTTACGACCGTTGTACGGTGCtt                                                                         <  1:62187/149‑1 (MQ=255)
                                                                                gggCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCttttt                                                                      >  2:289065/1‑149 (MQ=255)
                                                                                gggCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCttttt                                                                      <  1:435231/149‑1 (MQ=255)
                                                                                 ggCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTc                                                                     >  2:31809/1‑149 (MQ=255)
                                                                                  gCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCTTTGTACGGTGCTTTTTcc                                                                    >  2:37069/1‑148 (MQ=255)
                                                                                   cGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGc                                                                   <  1:29281/149‑1 (MQ=255)
                                                                                    ggTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGcc                                                                  <  2:489698/149‑1 (MQ=255)
                                                                                         ttctGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCagag                                                             >  1:104480/1‑149 (MQ=255)
                                                                                          tctGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATTAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGc                                                            >  2:585315/1‑149 (MQ=255)
                                                                                              gAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGcttct                                                        <  1:832355/149‑1 (MQ=255)
                                                                                               aGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGcttcct                                                       <  2:512862/149‑3 (MQ=255)
                                                                                                ggTGGCGGCGCATCCGTCCGCTGGTAACTCAACCATTAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATATAGTCAGTCTTATGCCTGGCCTTTTATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCttcttc                                                      <  1:549243/149‑1 (MQ=255)
                                                                                                ggTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGGCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCttcttc                                                      <  2:370363/149‑1 (MQ=255)
                                                                                                ggTGGCGGCGCATCCGGACGCTGGTAACTAAACCATAAAAAATGCCAGCCCGGATGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCttcttc                                                      <  2:856102/149‑1 (MQ=255)
                                                                                                   ggcggcGCATCCGGCCGCTGGTAACACAACCATAAAAAATGCCAGCACGGAGGCTGGCATTTTTAAATCAGATAAAGTCACTCTTATGACTGGCCTATGATCTCTTTACGACCGTTGTACGGTGCATTATAGCACATAGCTACTTCGaa                                                   >  2:310290/1‑148 (MQ=255)
                                                                                                    gcggcgCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATa                                                  <  2:305787/149‑1 (MQ=255)
                                                                                                    gcggcgCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATa                                                  <  1:482294/149‑1 (MQ=255)
                                                                                                    gcggcgCATCCGGCCGCTGGTAACTCAACCATAAAAAATGACAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATa                                                  <  2:114509/149‑1 (MQ=255)
                                                                                                     cggcgCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATAc                                                 <  2:626114/149‑1 (MQ=255)
                                                                                                         gcATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAAt                                             <  1:98457/149‑1 (MQ=255)
                                                                                                             ccGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGc                                         <  1:499427/149‑1 (MQ=255)
                                                                                                                  cGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGtt                                    <  1:787629/149‑1 (MQ=255)
                                                                                                                            tCAACCATAAAAAATGCCAGCCCGGAGGCTGGCAGTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTGTGATCTCTTTACGACCGTTGGACGGTGCTTTTTCGCCCAGCGCTTCTTCGCTACGAATCAGCGGGATGTATTTAGCa                          >  2:451209/1‑149 (MQ=255)
                                                                                                                                 cATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACg                     <  2:532017/149‑1 (MQ=255)
                                                                                                                                  tttaaaaaTGCCAGCCTCCGGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGATTCAGCTGGTTGTATTTAGCAACACgg                    >  2:260547/4‑149 (MQ=255)
                                                                                                                                    aaaaaaTGCCTGGCCGGTGGCGGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTTATCTCTTTACGACCGATGTACGGTGCTTTTTCGCCCAGAGCTTATTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTc                  <  2:815370/149‑1 (MQ=255)
                                                                                                                                    aaaaaaTGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGGCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTACGATACGAATCAGCTGGTTGTATTTAGCAACACGGTc                  <  1:620225/149‑1 (MQ=255)
                                                                                                                                     aaaaaTGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGGCAGTCTTATGCCTGGCCTTTGAGCTCTTTACGATCGTTGTCCGGGGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACCCGGTCa                 >  2:100888/1‑149 (MQ=255)
                                                                                                                                      aaaaTGCCAGCCCGGCGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACAGTTGTACGGTGCTTTTTCGCCCTGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTGAg                >  1:103127/1‑149 (MQ=255)
                                                                                                                                       aaaTGCCAGCCTCCGGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGa               >  2:105425/1‑149 (MQ=255)
                                                                                                                                            ccAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGc          <  2:131749/149‑1 (MQ=255)
                                                                                                                                             acgCCCGGAGGCTGGCATTTTTAAATCAGAATCAGTCAGGCTTATGCCTGGCCTTTGATCTCTTTCCGACCGTTGTACGGTGCGTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCt         <  1:504500/147‑1 (MQ=255)
                                                                                                                                                 caggAGGCTGGCATTTTTAATTTAGATAAAGTCAGTCTTATGCCTGGCCTTTGTTCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGGAACACGGTCAGAACGGCTCATa     <  2:707941/147‑1 (MQ=255)
                                                                                                                                                  cggcggCTGGCATTATTAAATAAGTTAAAGTCAGTATTATTCCTGGCCTTTGATATCTTTAAGACCGTTGAACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATAg    <  2:954353/149‑1 (MQ=255)
                                                                                                                                                    gAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATAGaa  <  1:289065/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATAGAA  >  NC_000913/2906461‑2906757

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: