Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,466,135 G→C 6.0% L16V (CTG→GTG)  treB ← trehalose‑specific PTS enzyme: IIB and IIC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,466,1350GC 6.0% 574.8 / 17.6 161L16V (CTG→GTG) treBtrehalose‑specific PTS enzyme: IIB and IIC component
Reads supporting (aligned to +/- strand):  ref base G (79/72);  new base C (5/5);  total (84/77)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 162 total aligned reads displayed.

CTGAAATTGCCCGGCATTGGTGAAACAGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGAAGATCACAATTTA  >  NC_000913/4465986‑4466283
                                                                                                                                                     |                                                                                                                                                    
cTGAAATTGCCCGGCATTGGTGAAACAGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAg                                                                                                                                                      >  2:1464424/1‑150 (MQ=255)
  gAAATTGCCCGGCATTGGTGAAACAGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGtt                                                                                                                                                    >  1:2770336/1‑150 (MQ=255)
  gAAATTGCCCGGCATTGGTGAAACAGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGtt                                                                                                                                                    >  2:2506026/1‑150 (MQ=255)
    aaTTGCCCGGCATTGGTGAAACAGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTtca                                                                                                                                                  >  1:2760794/1‑150 (MQ=255)
     aTTGCCCGGCATTGGTGAAACAGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTtcaa                                                                                                                                                 <  1:1215810/150‑1 (MQ=255)
       tGCCCGGCATTGGTGAAACAGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTtcaatc                                                                                                                                               <  2:2710147/150‑1 (MQ=255)
         cccGGCATTGGTGAAACCGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTtcaatcaa                                                                                                                                             <  2:768146/150‑1 (MQ=255)
          ccGGCATTGGTGAAACAGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAAc                                                                                                                                            >  1:2751271/1‑150 (MQ=255)
           cGGCATTGGTGAAACAGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAAcc                                                                                                                                           >  2:2142486/1‑150 (MQ=255)
           cGGCATTGGTGAAACAGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAAcc                                                                                                                                           >  2:2014372/1‑150 (MQ=255)
            ggCATTGGTGAAACAGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACcg                                                                                                                                          >  1:703068/1‑150 (MQ=255)
            ggCATTGGTGAAACAGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACcg                                                                                                                                          >  2:2187332/1‑150 (MQ=255)
                      aaaCAGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATcc                                                                                                                                <  1:2078304/150‑1 (MQ=255)
                       aaCAGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCg                                                                                                                               >  2:1165814/1‑150 (MQ=255)
                       aaCAGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCg                                                                                                                               >  1:2781627/1‑150 (MQ=255)
                       aaCAGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCg                                                                                                                               >  2:2508878/1‑150 (MQ=255)
                          aGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTTAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGttt                                                                                                                            >  1:1849909/1‑150 (MQ=255)
                          aGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGttt                                                                                                                            <  2:1022576/150‑1 (MQ=255)
                          aGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGttt                                                                                                                            >  1:2369944/1‑150 (MQ=255)
                           gCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTg                                                                                                                           <  2:2429808/150‑1 (MQ=255)
                                  aCCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTAtt                                                                                                                    <  2:1556576/150‑1 (MQ=255)
                                   ccATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTAttt                                                                                                                   <  1:2296595/150‑1 (MQ=255)
                                       aGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCt                                                                                                               <  1:2889803/150‑1 (MQ=255)
                                       aGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCt                                                                                                               <  1:2555258/150‑1 (MQ=255)
                                         ggAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCtca                                                                                                             <  1:656825/150‑1 (MQ=255)
                                         ggAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCtca                                                                                                             <  1:3083833/150‑1 (MQ=255)
                                         ggAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCtca                                                                                                             <  1:384602/150‑1 (MQ=255)
                                         ggAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCtca                                                                                                             <  1:2418919/150‑1 (MQ=255)
                                         ggAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCtca                                                                                                             <  1:1818536/150‑1 (MQ=255)
                                          gAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCtcat                                                                                                            >  2:2650245/1‑150 (MQ=255)
                                          gAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCtcat                                                                                                            <  1:2145288/150‑1 (MQ=255)
                                               gCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATaa                                                                                                       <  1:2133197/150‑1 (MQ=255)
                                                cTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATaaa                                                                                                      <  2:917857/150‑1 (MQ=255)
                                                 tCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAg                                                                                                     <  1:383049/150‑1 (MQ=255)
                                                   aaTTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCCGTTCAATCAACCGATCGATATCCGTTTGGGTTATTTTGCTCCTCATAAAGcc                                                                                                   >  1:346057/1‑150 (MQ=255)
                                                   aaTTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGcc                                                                                                   >  2:220143/1‑150 (MQ=255)
                                                   aaTTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGcc                                                                                                   >  2:1138949/1‑150 (MQ=255)
                                                   aaTTTCTTTCGGTCTGGCATTGGCCGGGTGGTTGAGGACAAAGCGTAGGCGCGGAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGGTTGGTTTATGTTGCTCATCATAAAGcc                                                                                                   >  1:346106/1‑150 (MQ=255)
                                                      ttctttcGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCa                                                                                                <  2:2661834/150‑1 (MQ=255)
                                                      ttctttcGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCa                                                                                                <  2:1119618/150‑1 (MQ=255)
                                                         tttcGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATgg                                                                                             >  2:1926087/1‑150 (MQ=255)
                                                         tttcGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATgg                                                                                             >  2:2947782/1‑150 (MQ=255)
                                                         tttcGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATgg                                                                                             <  2:1708232/150‑1 (MQ=255)
                                                         tttcGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCACTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATgg                                                                                             >  2:3039301/1‑150 (MQ=255)
                                                         tttcGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCACTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATgg                                                                                             >  2:3039249/1‑150 (MQ=255)
                                                          ttcGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGc                                                                                            <  1:1442379/150‑1 (MQ=255)
                                                           tcGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCa                                                                                           <  2:833019/150‑1 (MQ=255)
                                                           tcGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCa                                                                                           <  2:833232/150‑1 (MQ=255)
                                                             ggTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGa                                                                                         >  2:254691/1‑150 (MQ=255)
                                                              gTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGAt                                                                                        <  1:2387427/150‑1 (MQ=255)
                                                                       tGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACAttttt                                                                               <  1:589325/150‑1 (MQ=255)
                                                                        ggCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTttg                                                                              >  2:2139780/1‑150 (MQ=255)
                                                                        ggCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTttg                                                                              >  2:1108919/1‑150 (MQ=255)
                                                                             gttggttgAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCTTAAAGCCCCATGGCAGATGACATTTttggttgg                                                                         <  2:1761472/150‑1 (MQ=255)
                                                                               tggttgAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCt                                                                       <  2:2858277/150‑1 (MQ=255)
                                                                                ggttgAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTg                                                                      <  1:2508878/150‑1 (MQ=255)
                                                                                  ttgAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCa                                                                    <  2:1907023/150‑1 (MQ=255)
                                                                                  ttgAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCa                                                                    >  1:1293263/1‑150 (MQ=255)
                                                                                  ttgAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCa                                                                    >  1:307086/1‑150 (MQ=255)
                                                                                  ttgAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCa                                                                    >  1:307143/1‑150 (MQ=255)
                                                                                  ttgAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCa                                                                    >  1:368049/1‑150 (MQ=255)
                                                                                    gAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGa                                                                  <  1:2713790/150‑1 (MQ=255)
                                                                                     aGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGac                                                                 >  1:1909352/1‑149 (MQ=255)
                                                                                     aGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGaa                                                                 >  2:1272498/1‑150 (MQ=255)
                                                                                      ggACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAat                                                                >  2:3088513/1‑150 (MQ=255)
                                                                                       gACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAata                                                               <  2:3071238/150‑1 (MQ=255)
                                                                                          aaaGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAatatat                                                            >  2:2505793/1‑150 (MQ=255)
                                                                                                   gCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGaa                                                   >  2:542014/1‑150 (MQ=255)
                                                                                                       gTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAAt                                               <  2:975823/150‑1 (MQ=255)
                                                                                                        tAATACAGTGGCTCACCGTCGCACTATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAAtt                                              <  2:472793/150‑1 (MQ=255)
                                                                                                           tACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTaa                                           >  2:2387202/1‑150 (MQ=255)
                                                                                                           tACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTaa                                           <  1:11137/150‑1 (MQ=255)
                                                                                                           tACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCAAGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGACAAATTTaa                                           >  1:1045818/1‑150 (MQ=255)
                                                                                                              aGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTaaaaa                                        <  2:231028/150‑1 (MQ=255)
                                                                                                               gTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAAc                                       <  1:908879/150‑1 (MQ=255)
                                                                                                               gTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAAc                                       <  1:908836/150‑1 (MQ=255)
                                                                                                                 ggCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACgg                                     <  2:321432/150‑1 (MQ=255)
                                                                                                                   cTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGa                                   >  2:2903243/1‑150 (MQ=255)
                                                                                                                    tCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAATATTTAAAAACGGGaa                                  >  2:2472913/1‑150 (MQ=255)
                                                                                                                    tCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGaa                                  <  1:2947782/150‑1 (MQ=255)
                                                                                                                    tCACCGTCGCAATATTGCCGCGCCCGCCGACCACTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGaa                                  <  1:3039313/150‑1 (MQ=255)
                                                                                                                    tCACCGTCGCAATATTGCCGCGCCCGCCGACCACTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGaa                                  <  1:3039301/150‑1 (MQ=255)
                                                                                                                    tCACCGTCGCAATATTGCCGCGCCCGCCGACCACTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGaa                                  <  1:3039249/150‑1 (MQ=255)
                                                                                                                     cACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGATGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAAc                                 <  1:235938/150‑1 (MQ=255)
                                                                                                                     cACCGTCGCAATATTGCCGCGCCCGCCGACCACTTCAATCAACCGATCGATATCCGATTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAAc                                 <  2:141096/150‑1 (MQ=255)
                                                                                                                             cAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCTAGCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGTCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGa                         >  2:1281969/1‑150 (MQ=255)
                                                                                                                             cAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTGGGTTTATTTTGCTCATCAGAAAGCCCCATGGCAGATGACATTTTTGGGTGGCTGCAGACGCGAGTGCGCGGACAAATTTAAAAACGGGACCGTTCCCGa                         >  2:2003008/1‑150 (MQ=255)
                                                                                                                                tatTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAAc                      >  1:234564/1‑150 (MQ=255)
                                                                                                                                tatTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAAc                      >  2:1315712/1‑150 (MQ=255)
                                                                                                                                      cgcgCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCg                >  2:63006/1‑150 (MQ=255)
                                                                                                                                       gcgcCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATGTAAAAACGGGAACGTTCCCGAAACGCAGCGa               >  2:2298347/1‑150 (MQ=255)
                                                                                                                                           ccgccgACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGAAGAt           <  2:2781627/150‑1 (MQ=255)
                                                                                                                                            cgccgACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGAAGATc          >  1:2740379/1‑150 (MQ=255)
                                                                                                                                            cgccgACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGAAGATc          >  1:1422468/1‑150 (MQ=255)
                                                                                                                                             gccgACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGGAGATca         >  2:1794023/1‑150 (MQ=255)
                                                                                                                                                gACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGAAGATCACaa      >  1:2179065/1‑150 (MQ=255)
                                                                                                                                                gACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGAAGATCACaa      >  1:1041361/1‑150 (MQ=255)
                                                                                                                                                 aCCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGAAGATCACAAt     >  1:187144/1‑150 (MQ=255)
                                                                                                                                                   cAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGAAGATCACAAttt   >  1:1607336/1‑150 (MQ=255)
                                                                                                                                                    aGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCTGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGAAGATCACAATTTa  <  1:168383/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTGAAATTGCCCGGCATTGGTGAAACAGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTGCTCATCATAAAGCCCCATGGCAGATGACATTTTTGGTTGGCTGCAGAATATATTGCGCGGAAAAATTTAAAAACGGGAACGTTCCCGAAACGCAGCGAAGATCACAATTTA  >  NC_000913/4465986‑4466283

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