Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,297,196 A→C 5.7% intergenic (‑71/+23) yjcO ← / ← fdhF Sel1 repeat‑containing protein YjcO/formate dehydrogenase H

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,297,1960AC 5.7% 429.7 / 15.0 141intergenic (‑71/+23)yjcO/fdhFSel1 repeat‑containing protein YjcO/formate dehydrogenase H
Reads supporting (aligned to +/- strand):  ref base A (68/64);  new base C (6/2);  total (74/67)
Fisher's exact test for biased strand distribution p-value = 2.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.82e-01

Only 100 of 149 total aligned reads displayed.

AATACTGGCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGG  >  NC_000913/4297047‑4297341
                                                                                                                                                     |                                                                                                                                                 
 aTACTGGCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAag                                                                                                                                                  <  2:1000410/150‑1 (MQ=255)
   aCTGGCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAagga                                                                                                                                                >  1:849621/1‑150 (MQ=255)
    cTGGCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGTAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAaggag                                                                                                                                               <  1:1306543/150‑1 (MQ=255)
    cTGGCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAaggag                                                                                                                                               >  2:703485/1‑150 (MQ=255)
      ggCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGc                                                                                                                                             >  1:918334/1‑150 (MQ=255)
      ggCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGc                                                                                                                                             >  1:1313696/1‑150 (MQ=255)
      ggCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGc                                                                                                                                             <  2:1404225/150‑1 (MQ=255)
       gCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCt                                                                                                                                            <  2:1008829/150‑1 (MQ=255)
        cTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTg                                                                                                                                           >  2:1883163/1‑150 (MQ=255)
         tGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCTTTCGGAGGCTGt                                                                                                                                          <  1:1670073/150‑1 (MQ=255)
         tGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGt                                                                                                                                          >  2:2122254/1‑150 (MQ=255)
            cAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGa                                                                                                                                       >  2:436178/1‑150 (MQ=255)
             aGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGaa                                                                                                                                      >  2:1675371/1‑150 (MQ=255)
              ggCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGaaa                                                                                                                                     >  1:2103667/1‑150 (MQ=255)
               gCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAg                                                                                                                                    <  2:2189870/150‑1 (MQ=255)
                     aGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGt                                                                                                                              <  1:222352/150‑1 (MQ=255)
                        cGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTAtt                                                                                                                           <  1:1012280/150‑1 (MQ=255)
                         gTTGGCGTGGGCAAAGAATGTCAAAAACAAAATCAAAGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTa                                                                                                                          >  2:611634/1‑150 (MQ=255)
                            ggCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGc                                                                                                                       >  2:2383825/1‑150 (MQ=255)
                              cGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCa                                                                                                                     <  2:2091854/150‑1 (MQ=255)
                               gTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAg                                                                                                                    <  1:2048022/150‑1 (MQ=255)
                               gTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAg                                                                                                                    >  1:1651593/1‑150 (MQ=255)
                               gTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAg                                                                                                                    <  1:1457848/150‑1 (MQ=255)
                                tGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGt                                                                                                                   <  1:726668/150‑1 (MQ=255)
                                 gggCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTg                                                                                                                  <  1:2220641/150‑1 (MQ=255)
                                  ggCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGc                                                                                                                 >  2:511048/1‑150 (MQ=255)
                                   gCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGcc                                                                                                                <  1:770683/150‑1 (MQ=255)
                                            ttCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCgcgc                                                                                                       <  1:465515/149‑1 (MQ=255)
                                              cAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAg                                                                                                     >  1:315498/1‑150 (MQ=255)
                                               aaaaaCAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAgg                                                                                                    <  1:511048/150‑1 (MQ=255)
                                                      acaTCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGttt                                                                                             >  2:2004737/1‑150 (MQ=255)
                                                        aTCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTc                                                                                           >  2:2474345/1‑150 (MQ=255)
                                                        aTCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTc                                                                                           <  2:918334/150‑1 (MQ=255)
                                                         tCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCgcgc                                                                                                       <  2:1849880/137‑1 (MQ=255)
                                                          cAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCaa                                                                                         <  2:1313696/150‑1 (MQ=255)
                                                           taTGCGTTAATTTTTTTCATTTTGATAATACGATTGGTGAGGGAGTAACCAGTATAAAGAGTTGGAATGGATAATAAAACAGCCTCCGAATGGAGGCTGTAGAAAGGACGTTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAAc                                                                                        <  1:69424/149‑1 (MQ=255)
                                                              gCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACttg                                                                                     <  2:622427/150‑1 (MQ=255)
                                                               cGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACttgt                                                                                    <  2:2141034/150‑1 (MQ=255)
                                                                gATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACttgtt                                                                                   <  2:561389/150‑1 (MQ=255)
                                                                 aTAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACttgttg                                                                                  <  1:1764167/150‑1 (MQ=255)
                                                                    aTTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTAc                                                                               >  2:920264/1‑150 (MQ=255)
                                                                         tttCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTc                                                                          <  2:1941180/150‑1 (MQ=255)
                                                                         tttCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTc                                                                          <  2:1351843/150‑1 (MQ=255)
                                                                           tCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCaa                                                                        >  2:853236/1‑150 (MQ=255)
                                                                           tCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCACAGGCGAGTTTTCAACTTGTTGTACTCGTCaa                                                                        <  1:1848923/150‑1 (MQ=255)
                                                                             aTTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGATTGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGAAGGTATTACGCCAGTGCAGCTTCGCGCAGGCGAGTTTTCAACTTGTAGTACTCGTCAAta                                                                      >  1:1772388/1‑149 (MQ=255)
                                                                             aTTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATc                                                                      >  2:581512/1‑150 (MQ=255)
                                                                             aTTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATc                                                                      >  1:750378/1‑150 (MQ=255)
                                                                             aTTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATc                                                                      >  2:1439824/1‑150 (MQ=255)
                                                                              ttttGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCa                                                                     <  1:436178/150‑1 (MQ=255)
                                                                              ttttGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCa                                                                     >  1:192660/1‑150 (MQ=255)
                                                                              ttttGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCa                                                                     <  1:1439824/150‑1 (MQ=255)
                                                                              ttttGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCa                                                                     <  2:849621/150‑1 (MQ=255)
                                                                                ttGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACg                                                                   >  2:1575032/1‑150 (MQ=255)
                                                                                   ataataCGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTAc                                                                >  2:2001308/1‑150 (MQ=255)
                                                                                   ataataCGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTAc                                                                >  1:1426577/1‑150 (MQ=255)
                                                                                    taataCGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCGATCACGTACt                                                               <  2:1241582/150‑1 (MQ=255)
                                                                                    taataCGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACt                                                               <  2:362143/150‑1 (MQ=255)
                                                                                    taataCGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACt                                                               >  1:278806/1‑150 (MQ=255)
                                                                                          gATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTcgg                                                         >  1:2436723/1‑150 (MQ=255)
                                                                                             tGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTcggcgg                                                      <  1:2004737/150‑1 (MQ=255)
                                                                                              gggTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCggcggc                                                     >  1:1387764/1‑150 (MQ=255)
                                                                                               ggTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGc                                                             >  1:876239/1‑141 (MQ=255)
                                                                                               ggTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGgcggcg                                                    >  2:795985/1‑150 (MQ=255)
                                                                                                   gggAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCt                                                <  1:237176/150‑1 (MQ=255)
                                                                                                   gggAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCt                                                >  1:2093639/1‑150 (MQ=255)
                                                                                                     gAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTga                                              >  2:1484858/1‑150 (MQ=255)
                                                                                                      aGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTgat                                             >  1:1042352/1‑150 (MQ=255)
                                                                                                        tAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTgatcg                                           >  2:1396380/1‑150 (MQ=255)
                                                                                                         aaCCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGTTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCt                                                            >  2:1584096/1‑132 (MQ=255)
                                                                                                         aaCCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGTTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCt                                                            <  1:1584096/132‑1 (MQ=255)
                                                                                                         aaCCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTgatcgg                                          >  2:291883/1‑150 (MQ=255)
                                                                                                          aCCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTa                                                                                                                          <  2:1432913/69‑1 (MQ=255)
                                                                                                           ccAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTgatcggcg                                        <  1:2122254/150‑1 (MQ=255)
                                                                                                           ccAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTgatcggcg                                        <  2:359242/150‑1 (MQ=255)
                                                                                                            cAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTgatcggcga                                       <  2:780516/150‑1 (MQ=255)
                                                                                                              gTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTgatcggcgatc                                     <  2:738088/150‑1 (MQ=255)
                                                                                                                ataAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTgatcggcgatcgg                                   >  2:1645894/1‑150 (MQ=255)
                                                                                                                 taAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTgatcggcgatcggc                                  >  2:2218323/1‑150 (MQ=255)
                                                                                                                 taAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTgatcggcgatcggc                                  >  1:977298/1‑150 (MQ=255)
                                                                                                                  aaaGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCt                                 >  1:1085363/1‑150 (MQ=255)
                                                                                                                      agTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGAc                             <  2:1485850/150‑1 (MQ=255)
                                                                                                                      agTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGAc                             <  1:603469/150‑1 (MQ=255)
                                                                                                                               ggATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCCCGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGgc                    >  1:1431636/1‑150 (MQ=255)
                                                                                                                                gATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGgcg                   >  1:512908/1‑150 (MQ=255)
                                                                                                                                  tAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCa                 >  2:2448139/1‑150 (MQ=255)
                                                                                                                                  tAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCa                 >  1:189421/1‑150 (MQ=255)
                                                                                                                                   ctttctaCAGCCTCCTTTCGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAg                >  2:1006045/7‑150 (MQ=255)
                                                                                                                                   aaaaaaaCAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAg                >  1:217064/1‑150 (MQ=255)
                                                                                                                                   aaaaaaaCAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAg                <  2:1211263/150‑1 (MQ=255)
                                                                                                                                      tctaCAGCCTCCTTTCGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTAt             >  1:585700/4‑150 (MQ=255)
                                                                                                                                       ctaCAGCCTCCTTTCGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTAtt            >  2:2296100/3‑150 (MQ=255)
                                                                                                                                       ctaCAGCCTCCTTTCGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTAtt            >  1:1007928/3‑150 (MQ=255)
                                                                                                                                        taCAGCCTCCTTTCGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTAttt           >  2:2453593/2‑150 (MQ=255)
                                                                                                                                        taCAGCCTCCTTTCGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTAttt           >  2:1882777/2‑150 (MQ=255)
                                                                                                                                          cAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGt         <  1:136352/150‑1 (MQ=255)
                                                                                                                                            gCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTAc       >  1:1834895/1‑150 (MQ=255)
                                                                                                                                               tCCGAAAGGAGGCTGTGGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTcc    >  2:1819856/1‑150 (MQ=255)
                                                                                                                                               tCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTcc    >  1:1595130/1‑150 (MQ=255)
                                                                                                                                                ccGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCg   >  1:1117451/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
AATACTGGCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGG  >  NC_000913/4297047‑4297341

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