Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,776,173 A→G Y196C (TAC→TGC)  ydiF → putative acetyl‑CoA:acetoacetyl‑CoA transferase: alpha subunit/beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,776,1730AG99.0% 385.1 / ‑4.4 102Y196C (TAC→TGC) ydiFputative acetyl‑CoA:acetoacetyl‑CoA transferase: alpha subunit/beta subunit
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (51/50);  total (52/50)
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.18e-01

Only 100 of 102 total aligned reads displayed.

CCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGG  >  NC_000913/1776025‑1776317
                                                                                                                                                    |                                                                                                                                                
ccACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCg                                                                                                                                                >  2:926547/1‑151 (MQ=255)
   cGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCa                                                                                                                                             <  1:1187469/151‑1 (MQ=255)
   cGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCa                                                                                                                                             <  2:1385959/151‑1 (MQ=255)
   cGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCa                                                                                                                                             <  1:139368/151‑1 (MQ=255)
         cAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTg                                                                                                                                       >  2:564429/1‑151 (MQ=255)
           aGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATGACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGa                                                                                                                                      >  2:1243201/1‑150 (MQ=255)
             gcggcAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGa                                                                                                                                      >  1:1194633/1‑148 (MQ=255)
                gcAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATg                                                                                                                                 >  2:3520/1‑150 (MQ=255)
                 cAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATg                                                                                                                                 >  2:2497/1‑149 (MQ=255)
                      tgaatgaaGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGAt                                                                                                                          >  1:538086/1‑151 (MQ=255)
                       gaatgaaGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGa                                                                                                                           <  1:290298/149‑1 (MQ=255)
                       gaatgaaGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGa                                                                                                                           <  2:1650815/149‑1 (MQ=255)
                       gaatgaaGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATg                                                                                                                         >  2:1141055/1‑151 (MQ=255)
                           gaaGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTAt                                                                                                                      <  1:761680/150‑1 (MQ=255)
                             aGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTAtctc                                                                                                                   >  2:818993/1‑151 (MQ=255)
                               tCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCg                                                                                                                  >  1:1466650/1‑150 (MQ=255)
                               tCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCg                                                                                                                  >  1:1154964/1‑150 (MQ=255)
                                cACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGa                                                                                                                 <  2:254737/150‑1 (MQ=255)
                                  cTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGc                                                                                                              >  1:1596379/1‑151 (MQ=255)
                                    aaaGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCa                                                                                                             <  1:373258/150‑1 (MQ=255)
                                    aaaGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCa                                                                                                             <  2:666993/150‑1 (MQ=255)
                                    aaaGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCAt                                                                                                            >  2:253137/1‑151 (MQ=255)
                                     aagaagACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCAtt                                                                                                           >  1:1141221/1‑151 (MQ=255)
                                      agaagaCCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTg                                                                                                          <  2:1417576/151‑1 (MQ=255)
                                        aagaCCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGt                                                                                                        >  2:1116328/1‑151 (MQ=255)
                                            ccTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTAtt                                                                                                    <  2:143133/151‑1 (MQ=255)
                                                aTTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGccc                                                                                                <  1:422387/151‑1 (MQ=255)
                                                aTTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGccc                                                                                                <  2:1202639/151‑1 (MQ=255)
                                                 ttAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGccc                                                                                                >  2:1705003/1‑150 (MQ=255)
                                                  tAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGccc                                                                                                >  2:1107459/1‑149 (MQ=255)
                                                  tAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGccc                                                                                                >  2:1543350/1‑149 (MQ=255)
                                                  tAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCa                                                                                               <  1:668350/150‑1 (MQ=255)
                                                   aaaCTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAg                                                                                              <  1:1067236/150‑1 (MQ=255)
                                                   aaaCTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAg                                                                                              <  1:1167907/150‑1 (MQ=255)
                                                    aaCTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGGTATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGc                                                                                            <  2:726108/151‑1 (MQ=255)
                                                    aaCTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGc                                                                                            >  2:1491507/1‑151 (MQ=255)
                                                     aCTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCg                                                                                           >  1:1218503/1‑151 (MQ=255)
                                                        ggTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCgg                                                                                          <  1:463218/149‑1 (MQ=255)
                                                        ggTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTg                                                                                        <  2:1707007/151‑1 (MQ=255)
                                                           cGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCCCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGcaca                                                                                    <  1:113903/151‑1 (MQ=255)
                                                             aGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGcaca                                                                                    >  1:1505644/1‑150 (MQ=255)
                                                                 tGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATaa                                                                                >  2:745495/1‑150 (MQ=255)
                                                                  gATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATaa                                                                                <  1:845643/149‑1 (MQ=255)
                                                                  gATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAAcg                                                                              <  1:201474/151‑1 (MQ=255)
                                                                  gATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAAcg                                                                              >  1:520718/1‑151 (MQ=255)
                                                                    tAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAAcggc                                                                            >  1:1604940/1‑151 (MQ=255)
                                                                     aaCAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAAcggcg                                                                           <  2:1212918/151‑1 (MQ=255)
                                                                     aaCAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAAcggc                                                                            <  2:821453/150‑1 (MQ=255)
                                                                     aaCAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCCCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAAcggc                                                                            >  1:1428967/1‑150 (MQ=255)
                                                                      aCAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAAcggcgg                                                                          >  1:1508838/1‑151 (MQ=255)
                                                                        aaaGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTa                                                                        <  2:743406/151‑1 (MQ=255)
                                                                         aaGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTAt                                                                       <  2:402630/151‑1 (MQ=255)
                                                                            aaTATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATtgt                                                                    <  1:1546976/151‑1 (MQ=255)
                                                                             atatCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATtgt                                                                    >  1:800977/1‑150 (MQ=255)
                                                                             atatCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGAGCACAATAACGGCGGGAtt                                                                      >  1:775956/1‑148 (MQ=255)
                                                                               atCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGtga                                                                  <  2:55136/150‑1 (MQ=255)
                                                                                  tctATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTgtggtg                                                                >  1:63236/1‑149 (MQ=255)
                                                                                   ctATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGtgat                                                                 <  1:665076/147‑1 (MQ=255)
                                                                                    tATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGGATTGtgatggttc                                                            >  2:474281/1‑147 (MQ=255)
                                                                                     aTTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGc                                                            <  2:1169923/150‑1 (MQ=255)
                                                                                      ttACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAAcggcggattttt                                                                    >  1:928579/1‑135 (MQ=255)
                                                                                      ttACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTGTtttgtt                                                                 >  1:1304271/1‑142 (MQ=255)
                                                                                       tACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTGGTGTGATTATGCggg                                                         >  1:1240949/1‑151 (MQ=255)
                                                                                       tACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATtgtggatgagga                                                           >  2:839610/1‑141 (MQ=255)
                                                                                       tACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAg                                                          <  2:212049/150‑1 (MQ=255)
                                                                                       tACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGGATTGTGATGATGCAgg                                                         >  2:1298258/1‑151 (MQ=255)
                                                                                        aCAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAgg                                                         <  1:1491407/150‑1 (MQ=255)
                                                                                           aaGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTTTGATGATGCAGGTg                                                       >  1:1501711/1‑149 (MQ=255)
                                                                                           aaGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCCGGGGcc                                                     >  1:191608/1‑150 (MQ=255)
                                                                                              cGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTTTGGTGATGCAGGGGgaga                                                   >  2:1677791/1‑150 (MQ=255)
                                                                                                   gcgcCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATg                                              <  1:1388749/150‑1 (MQ=255)
                                                                                                    cgcCAGATATTGCCTTCATTCGCGCTACCGCCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGt                                            <  2:1505741/151‑1 (MQ=255)
                                                                                                    cgcCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGt                                            <  2:923183/151‑1 (MQ=255)
                                                                                                    cgcCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGGTGTAGGGGGAGaaaa                                                >  2:1287135/1‑147 (MQ=255)
                                                                                                        aGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGGATTGTGATTATGGAGGGGGAGAAAAAgggt                                           >  1:48445/1‑146 (MQ=255)
                                                                                                        aGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGCATTGTTATGATGTAGGGGCGGGAAATGGGTa                                          >  1:56126/1‑149 (MQ=255)
                                                                                                            atTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGaaa                                     <  2:444206/150‑1 (MQ=255)
                                                                                                                 cTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGGTGATGCAGGGGGAGgaaaaaggtaaaaaa                                     >  1:1060058/1‑134 (MQ=255)
                                                                                                                   tCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACg                               <  1:926480/149‑1 (MQ=255)
                                                                                                                   tCATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGACATGTGGTGGTGTAGGAGGGGAAAATGGGTAAGaaactcaagc                              >  1:1513718/1‑143 (MQ=255)
                                                                                                                   gcATTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTTTGATGATGCAGGGGCAGAAAAAGGTTaaaaaa                                     >  2:424508/2‑143 (MQ=255)
                                                                                                                     aTTCGCGNTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAANNNGNNGTNATGNNNNNNGACGCATNGGNTATTGCCCNGGCNNNGNANAATAACGGCGGGATTGTGAGGATGCGGGGGGAGAAAATGGGGAAGAAAACCACGccg                            >  2:739231/1‑148 (MQ=255)
                                                                                                                     aTTCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGc                           <  1:290317/151‑1 (MQ=255)
                                                                                                                       tCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCa                          <  2:369453/150‑1 (MQ=255)
                                                                                                                       tCGCGCTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGGGGAGAAAATGGGTAAGAAAAACACGCt                             >  1:66856/1‑147 (MQ=255)
                                                                                                                          cgcTACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCt                      <  2:1169018/151‑1 (MQ=255)
                                                                                                                             tACCACCTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGAGTATGAAGGGGGAGAAAATGGGTTAGAAAACCCCGCGGCa                          >  2:229770/1‑144 (MQ=255)
                                                                                                                              accaccTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTaaa                   <  2:1591504/150‑1 (MQ=255)
                                                                                                                              accaccTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTaaa                   <  1:533702/150‑1 (MQ=255)
                                                                                                                              accaccTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTaaa                   <  2:1596476/150‑1 (MQ=255)
                                                                                                                               ccaccTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATc                 <  2:1021669/151‑1 (MQ=255)
                                                                                                                               ccaccTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGGGCAGAAAAAGGTTAAGAAAACCACGCTGCATCCTaaa                   >  2:536229/1‑149 (MQ=255)
                                                                                                                                caccTGCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCt                <  1:179755/151‑1 (MQ=255)
                                                                                                                                     gCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTc             <  1:1087018/149‑1 (MQ=255)
                                                                                                                                     gCGACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCg           <  1:761303/151‑1 (MQ=255)
                                                                                                                                       gACAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTTTGATGATGCAGGGGCAGAAAAAGGGTAAGAAAACCACGCCGCATCCTAAATCTGTcc            >  2:1173035/1‑148 (MQ=255)
                                                                                                                                        aCAGTGAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGGGCAGAAAAAGGGTAAGAAAAACACGCCGCATCCTAAAACTGttc            >  1:1350008/1‑145 (MQ=255)
                                                                                                                                             gAAGGCTGCGCCTCTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAGAGCCACGCTGCATCCTAAATCTGTCCGTATTCCgg   <  2:1368258/151‑1 (MQ=255)
                                                                                                                                             gAAGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGGATTGTGATGATGCGGGTGCAGAAAATGGGTAAGAAAACCACGCTCCATCCTATATCTGTCCGTAtttc     >  1:156313/1‑147 (MQ=255)
                                                                                                                                               aGGCTGCGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAAAGGTTAAGAAAACCACGCTGCATCCTAAATCTGTCCCTATTCCggg  >  2:568172/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGG  >  NC_000913/1776025‑1776317

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