Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,633,634 C→T 24.0% E38K (GAG→AAG)  yjjX ← non‑canonical purine NTP phosphatase, ITPase/XTPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,633,6340CT24.0% 163.5 / 49.5 79E38K (GAG→AAG) yjjXnon‑canonical purine NTP phosphatase, ITPase/XTPase
Reads supporting (aligned to +/- strand):  ref base C (25/35);  new base T (7/12);  total (32/47)
Fisher's exact test for biased strand distribution p-value = 7.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGCTGTCACCATCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAA  >  NC_000913/4633495‑4633769
                                                                                                                                           |                                                                                                                                       
gTGCTGTCACCATCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCtctc                                                                                                                                         >  2:754873/1‑140 (MQ=255)
  gCTGTCACCATCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGa                                                                                                                                       <  2:1488371/140‑1 (MQ=255)
    tgacACCATCGATGCCAGCTTCAATCGCCACCCAAAACTCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACg                                                                                                                                     <  2:615732/137‑1 (MQ=255)
    tGTCACCATCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACg                                                                                                                                     >  2:1304864/1‑140 (MQ=255)
    tGTCACCATCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACg                                                                                                                                     >  2:764081/1‑140 (MQ=255)
       cACCATCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCa                                                                                                                                  <  1:1602824/140‑1 (MQ=255)
       cACCATCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCa                                                                                                                                  <  1:1605840/140‑1 (MQ=255)
        aCCATCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCaa                                                                                                                                 <  2:5035/140‑1 (MQ=255)
            tCGATGCCAGCTTCAATCGCGACCCAAAAATCAGCCTCTGGAAGTAACCGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCTGTACACCGCTCTCGATGGCAACGGa                                                                                                                             <  2:254946/140‑1 (MQ=255)
            tCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGa                                                                                                                             <  2:1426132/140‑1 (MQ=255)
            tCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGa                                                                                                                             <  2:517807/140‑1 (MQ=255)
            tCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGa                                                                                                                             <  2:275712/140‑1 (MQ=255)
            tCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGa                                                                                                                             <  2:274666/140‑1 (MQ=255)
            tCGATGCCAGATTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGa                                                                                                                             <  1:1029952/140‑1 (MQ=255)
              gATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAATCGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATg                                                                                                                           <  2:840909/140‑1 (MQ=255)
              gATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATg                                                                                                                           >  2:835301/1‑140 (MQ=255)
              gATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATg                                                                                                                           <  1:661623/140‑1 (MQ=255)
              gATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATg                                                                                                                           <  2:888066/140‑1 (MQ=255)
                tGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCa                                                                                                                         >  1:204022/1‑140 (MQ=255)
                tGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCa                                                                                                                         <  2:1188894/140‑1 (MQ=255)
                tGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCa                                                                                                                         <  2:730276/140‑1 (MQ=255)
                tGCCAGCTTCAATCGCCACCCAAAAATCAGCCGCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGcc                                                                                                                         >  2:1512061/1‑139 (MQ=255)
                   cAGCTTCAATCGCCACCCCAAAATCAGCGTCTGGACGTAATCGGCGGGCATTGGCTCCCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCGAAAGGGCTGTTCCGGTACACCGCTCGCGACGGCAACGGATGCActa                                                                                                                      <  2:344788/140‑1 (MQ=255)
                   cAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACCGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCGCTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAata                                                                                                                      <  2:658969/140‑1 (MQ=255)
                      cTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTACACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCACTATgg                                                                                                                   >  2:1144241/1‑140 (MQ=255)
                      cTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCACCGGATGCAATATgg                                                                                                                   >  2:900774/1‑140 (MQ=255)
                      cTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCACTCGACGGCAACGGATGCAATATgg                                                                                                                   >  2:1528890/1‑140 (MQ=255)
                        tCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCa                                                                                                                 <  1:238459/140‑1 (MQ=255)
                        tCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCa                                                                                                                 <  1:1219248/140‑1 (MQ=255)
                               ccACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCt                                                                                                          >  2:1141895/1‑140 (MQ=255)
                                cACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTTCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCtt                                                                                                         >  2:1347574/1‑140 (MQ=255)
                                cACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCtt                                                                                                         >  2:317807/1‑140 (MQ=255)
                                cACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCtt                                                                                                         >  2:1023324/1‑140 (MQ=255)
                                           cAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGAtc                                                                                              <  1:122670/140‑1 (MQ=255)
                                             gCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGAtctc                                                                                            >  2:346576/1‑140 (MQ=255)
                                             gCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGAtctc                                                                                            >  1:213084/1‑140 (MQ=255)
                                             gCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGAtctc                                                                                            >  1:213140/1‑140 (MQ=255)
                                                 ctcccaGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCCTTTCCTCACTGCCAAAGGGCAGTTCCGGTACACCGCACTCGACGGCAACGGATGCAATATGGTAGGATCCGTCGCCGAAGATATCGTGa                                                                                        <  2:1624426/135‑1 (MQ=255)
                                                 ctGGAAGTAAACGGCGGGCATTGGCTTCCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGa                                                                                        <  2:1649494/140‑1 (MQ=255)
                                                 ctGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGa                                                                                        <  2:783187/140‑1 (MQ=255)
                                                 ctGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGa                                                                                        <  2:851467/140‑1 (MQ=255)
                                                        tAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCt                                                                                 >  1:1532635/1‑140 (MQ=255)
                                                                             cagATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAg                                                            >  1:265500/3‑140 (MQ=255)
                                                                                  ttCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGa                                                       <  1:1134770/140‑1 (MQ=255)
                                                                                  ttCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGa                                                       <  1:589042/140‑1 (MQ=255)
                                                                                  ttCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGa                                                       <  1:1549354/140‑1 (MQ=255)
                                                                                      tGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATtgg                                                   >  2:1455971/1‑140 (MQ=255)
                                                                                               gcgcgTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGcac                                          <  2:513614/140‑1 (MQ=255)
                                                                                                cgcgTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGcaca                                         <  2:46722/140‑1 (MQ=255)
                                                                                                cgcgTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGcaca                                         <  2:46388/140‑1 (MQ=255)
                                                                                                  cgTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGa                                       <  2:1057583/140‑1 (MQ=255)
                                                                                                  cgTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGa                                       <  1:150756/140‑1 (MQ=255)
                                                                                                  cgTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGa                                       <  1:149227/140‑1 (MQ=255)
                                                                                                  cgTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGa                                       <  1:120656/140‑1 (MQ=255)
                                                                                                       ccTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACt                                  <  2:1467307/140‑1 (MQ=255)
                                                                                                        cTCACTGCCATAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTAGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACtt                                 <  2:989529/140‑1 (MQ=255)
                                                                                                        cTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCTATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACtt                                 <  2:1216205/140‑1 (MQ=255)
                                                                                                             ttccAAAGGGCAGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTAAGCGGGAATGGTGGTCGCACAGACCACTTGGTGc                            >  2:441378/3‑140 (MQ=255)
                                                                                                             tGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGc                            >  2:1175110/1‑140 (MQ=255)
                                                                                                                cAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATa                         <  2:1085466/140‑1 (MQ=255)
                                                                                                                  aaGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAAt                       <  2:1271729/140‑1 (MQ=255)
                                                                                                                   aGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTACGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCGGTATTTTAGCGGGATTGGTGGTCGCTCAGACGACTGGGTGCATAATc                      >  2:291806/1‑140 (MQ=255)
                                                                                                                   aGGGCTGTTCCGGTACACAGCTCTTCAAGGGAACGGATGCACTAAGGCAGGATCCTTCGCGGTAGATCTCGTGACATGCCTGCAGGATGGCCGGAGTGTTAGCGGGATGGGTGGTCGCTCAGACAACTTGGTGcctgagc                      >  2:1429678/1‑134 (MQ=255)
                                                                                                                     ggCTGTTCCGGTACACCGCACTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAg                    <  1:685662/140‑1 (MQ=255)
                                                                                                                       cTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCa                  >  1:734787/1‑140 (MQ=255)
                                                                                                                       cTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCa                  >  1:401209/1‑140 (MQ=255)
                                                                                                                         gTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCAtt                >  1:322796/1‑140 (MQ=255)
                                                                                                                         gTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCAtt                >  1:817327/1‑140 (MQ=255)
                                                                                                                              ggTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCa           <  2:1556519/140‑1 (MQ=255)
                                                                                                                                tACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGa         >  2:621433/1‑140 (MQ=255)
                                                                                                                                    ccGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAAt     >  1:1517702/1‑140 (MQ=255)
                                                                                                                                       ctctTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTaa  <  2:366552/140‑1 (MQ=255)
                                                                                                                                       ctctTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTaa  <  2:541167/140‑1 (MQ=255)
                                                                                                                                       ctctTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTaa  >  2:1525674/1‑140 (MQ=255)
                                                                                                                                       ctctTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTaa  <  2:366776/140‑1 (MQ=255)
                                                                                                                                       ctctCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGGATTTTAGCGGGAGTGGGGGTCGCACTGCCAACTTGGTGCATAATCAGCATTCCTCAGAAAATTaa  >  1:1482011/1‑140 (MQ=255)
                                                                                                                                       ctctCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTg                             <  2:1176797/113‑1 (MQ=255)
                                                                                                                                       ctctCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTg                             >  1:1176797/1‑113 (MQ=255)
                                                                                                                                       ctctCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATa                         >  2:590538/1‑117 (MQ=255)
                                                                                                                                       ctctCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATa                         <  1:590538/117‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                       
GTGCTGTCACCATCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAA  >  NC_000913/4633495‑4633769

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: