Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,906,608 G→C 6.5% intergenic (+25/+35) 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,6080GC 6.5% 246.8 / 5.6 75intergenic (+25/+35)ygcG/enoTPM domain protein, putative phosphatase/enolase
Reads supporting (aligned to +/- strand):  ref base G (26/44);  new base C (2/3);  total (28/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

AAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATAG  >  NC_000913/2906462‑2906755
                                                                                                                                                  |                                                                                                                                                   
aaaGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGa                                                                                                                                                   <  4:235962/149‑1 (MQ=255)
    agCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCt                                                                                                                                               <  4:121668/149‑1 (MQ=255)
     gCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTg                                                                                                                                              <  4:674870/149‑1 (MQ=255)
     gCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTg                                                                                                                                              <  3:7112/149‑1 (MQ=255)
     gCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTg                                                                                                                                              <  3:362388/149‑1 (MQ=255)
       aTCTGGTAGAAGTTATCATTCAGACAACGGCCGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGGGCTTCTGGAGGTGGCGGCGCCTCCGGCCCCTGGTAACTCAACCATAAAAAATGCCAGCCCGGGGGCCGGc                                                                                                                                            >  4:418456/1‑149 (MQ=255)
       aTCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGc                                                                                                                                            <  3:472654/149‑1 (MQ=255)
        tCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCa                                                                                                                                           >  1:275341/1‑149 (MQ=255)
          tGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCTCCGGGCTGGCAtt                                                                                                                                         <  3:637453/149‑1 (MQ=255)
           ggTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCAttt                                                                                                                                        <  3:95608/149‑1 (MQ=255)
           ggTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCAttt                                                                                                                                        <  4:50524/149‑1 (MQ=255)
               gAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTaa                                                                                                                                    <  3:625106/149‑1 (MQ=255)
                aaGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCATAAAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCCGGAGGTCGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCTCCGGGCTGGCAttttttat                                                                                                                                   <  4:592698/149‑4 (MQ=255)
                aaGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCTCCGGGCTGGCAttttttat                                                                                                                                   <  1:55721/149‑4 (MQ=255)
                aaGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTaaa                                                                                                                                   <  3:34167/149‑1 (MQ=255)
                aaGTTATCAATCAGACAACAGCGGTTACGCGGGAGGCTAAGATAGCGGCGGATTTAGCGGCGGGGGCGCTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATCAAAAATGCCAGCCCGGAGGCTGGCATTTTTaaa                                                                                                                                   <  2:536205/149‑1 (MQ=255)
                  gTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATc                                                                                                                                 >  2:397884/1‑149 (MQ=255)
                   ttATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGGTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCa                                                                                                                                >  3:383726/1‑149 (MQ=255)
                            aGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAagtc                                                                                                                       >  1:398743/1‑149 (MQ=255)
                                        ttCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCt                                                                                                           <  1:609510/149‑1 (MQ=255)
                                        ttCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCt                                                                                                           <  1:467907/149‑1 (MQ=255)
                                          ccGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATTCGCCCGCAGGTAACTACAGCATAAAAAATGCCCGCCCGGAGGCTGGCATTTTTAAATCAGGTAAAGTCAGTCTTATGGCttg                                                                                                         >  4:79362/1‑147 (MQ=255)
                                              tGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCtt                                                                                                     <  2:230249/149‑1 (MQ=255)
                                                    cTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGAtct                                                                                               <  3:220969/149‑1 (MQ=255)
                                                       aGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCtt                                                                                            <  2:191955/149‑1 (MQ=255)
                                                       aGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCtt                                                                                            <  4:624795/149‑1 (MQ=255)
                                                       aGATAGCGGCGGAATTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCCGGGGCTGGCATTTTTAAATCAGATAAAGTGAGTCTTATGCCTGGCCTTTGACCGCtt                                                                                            >  4:441935/1‑149 (MQ=255)
                                                        gATAGCGGCGGATTTAGCGGCGGGGGCCGGTCTTTTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGGCGGCATTTTTAAATCAGATAAAGTCCGTCTTATGCCTGGCCTTTGATCTCttt                                                                                           >  1:476030/1‑149 (MQ=255)
                                                          tAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTAc                                                                                         <  4:136174/149‑1 (MQ=255)
                                                              ggcggATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGAcc                                                                                     >  1:173697/1‑149 (MQ=255)
                                                               gcggATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGGTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCg                                                                                    >  3:65363/1‑149 (MQ=255)
                                                                 ggATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGtt                                                                                  <  3:620041/149‑1 (MQ=255)
                                                                      tAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACg                                                                             <  4:161052/149‑1 (MQ=255)
                                                                           gcggGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCt                                                                        >  4:200882/1‑149 (MQ=255)
                                                                              ggggCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCtttt                                                                     <  2:310675/149‑1 (MQ=255)
                                                                                ggaggTTCTTCTGGAGGTGGCGGCGAATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCACGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTc                                                                   <  2:358751/146‑1 (MQ=255)
                                                                                ggCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTc                                                                   <  4:127254/149‑1 (MQ=255)
                                                                                   ggTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGcc                                                                >  3:500850/1‑149 (MQ=255)
                                                                                     ttcttcTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCa                                                              >  2:410718/1‑149 (MQ=255)
                                                                                     ttcttcTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCa                                                              >  1:150117/1‑149 (MQ=255)
                                                                                       cttctGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCaga                                                            <  3:568007/149‑1 (MQ=255)
                                                                                       cttctGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCaga                                                            <  4:284648/149‑1 (MQ=255)
                                                                                         tctGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGc                                                          <  1:235521/149‑1 (MQ=255)
                                                                                            ggAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGcttc                                                       <  2:275341/149‑1 (MQ=255)
                                                                                              aGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGcttctt                                                     <  3:35084/149‑1 (MQ=255)
                                                                                               ggTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCTAGAGCttcttc                                                    <  3:441935/149‑1 (MQ=255)
                                                                                                  ggcggcGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGAt                                                 <  1:323893/149‑1 (MQ=255)
                                                                                                   gcggcgCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATa                                                >  3:302069/1‑149 (MQ=255)
                                                                                                   gcggcgCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATa                                                <  2:150117/149‑1 (MQ=255)
                                                                                                   gcggcgCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATa                                                >  1:378833/1‑149 (MQ=255)
                                                                                                   gcggcgCATCCGGCCGATGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATa                                                <  2:135166/149‑1 (MQ=255)
                                                                                                           tCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAg                                        >  2:291976/1‑149 (MQ=255)
                                                                                                             cGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCt                                      >  3:237035/1‑149 (MQ=255)
                                                                                                              ggCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTg                                     <  1:696148/149‑1 (MQ=255)
                                                                                                                    tGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTa                               <  2:572737/149‑1 (MQ=255)
                                                                                                                    tGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTa                               >  3:184434/1‑149 (MQ=255)
                                                                                                                           tCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAg                          >  2:542156/1‑147 (MQ=255)
                                                                                                                           tCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAg                          <  1:542156/147‑1 (MQ=255)
                                                                                                                           tCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCa                        >  2:15661/1‑149 (MQ=255)
                                                                                                                           tCAACCATAAAAAATGCCACCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCa                        >  3:563506/1‑149 (MQ=255)
                                                                                                                            cAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCaa                       <  1:291976/149‑1 (MQ=255)
                                                                                                                            cAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCaa                       <  1:397884/149‑1 (MQ=255)
                                                                                                                                  tAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGt                 >  1:462582/1‑149 (MQ=255)
                                                                                                                                  tAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGt                 >  1:425935/1‑149 (MQ=255)
                                                                                                                                  tAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGt                 >  3:477923/1‑149 (MQ=255)
                                                                                                                                      aaaTGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGa             <  3:156176/149‑1 (MQ=255)
                                                                                                                                       aaTGCCAGCCTCCGGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGaa            >  3:340437/1‑149 (MQ=255)
                                                                                                                                       aaTGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTAc                                                                                         >  2:607477/1‑72 (MQ=255)
                                                                                                                                       aaTGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTAc                                                                                         <  1:607477/72‑1 (MQ=255)
                                                                                                                                           ccAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGc        <  2:173697/149‑1 (MQ=255)
                                                                                                                                             aGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCt       <  2:646074/148‑1 (MQ=255)
                                                                                                                                             aGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCt       >  1:646074/1‑148 (MQ=255)
                                                                                                                                                ccGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATa   <  3:370220/149‑1 (MQ=255)
                                                                                                                                                ccGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATa   <  4:63666/149‑1 (MQ=255)
                                                                                                                                                 cGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATAg  <  1:452838/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATAG  >  NC_000913/2906462‑2906755

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