Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,653,135 A→C 7.6% K214N (AAA→AAC sseA → 3‑mercaptopyruvate sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,653,1350AC 7.6% 1308.0 / 19.0 445K214N (AAA→AACsseA3‑mercaptopyruvate sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base A (225/175);  new base C (2/43);  total (227/218)
Fisher's exact test for biased strand distribution p-value = 1.83e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 474 total aligned reads displayed.

CATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGG  >  NC_000913/2652986‑2653284
                                                                                                                                                     |                                                                                                                                                     
 aTGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGGACCTCGCCCAGGTTTACGTCGCGGACATAGTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTaaaaa                                                                                                                                                      <  1:552245/150‑1 (MQ=255)
 aTGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTaaaaa                                                                                                                                                      <  2:3052407/150‑1 (MQ=255)
 aTGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTaaaaa                                                                                                                                                      >  2:2567234/1‑150 (MQ=255)
    aaaaTACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGa                                                                                                                                                   >  1:7469507/1‑150 (MQ=255)
     aaaTACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGAc                                                                                                                                                  <  2:5465994/150‑1 (MQ=255)
        tACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGa                                                                                                                                               >  1:4340781/1‑150 (MQ=255)
            gcgcAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGaa                                                                                                                                           >  2:6521406/1‑150 (MQ=255)
            gcgcAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGaa                                                                                                                                           >  2:1123424/1‑150 (MQ=255)
              gccaaTTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGACCTGGTGCGCGAAGGCGAACTAAAAACGACCGATGaccg                                                                                                                                         >  1:5683513/4‑147 (MQ=255)
                 aaTTTTTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGATTGTTCCGTGTACGTAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGa                                                                                                                                      <  2:282198/150‑1 (MQ=255)
                  attattGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGAGGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGAt                                                                                                                                     <  2:1857339/150‑1 (MQ=255)
                  attattGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGAt                                                                                                                                     >  1:524124/1‑150 (MQ=255)
                   ttattGATGCCCGCCGGGGTGCACGTTTTAACGCAGAAGTTTATGACCCTCGCCCAGTTTTACGTCGCGGAAATATTCCCGGAGAACTTATTGTTCCGTGGACGTAACTGGTGCGCTAAGGCTAACTAAAAACGACCGATGAACTGGATg                                                                                                                                    <  2:7855920/150‑1 (MQ=255)
                    ttttGATGCCCGCCCGGCTGCACTTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATTTTCCCGGAGCACTGAATGTTCCGTGGACGTAAGTGGTGCGCGAAGGCTAACTAAAAACGACCGATGAACTGGATGc                                                                                                                                   <  2:1006826/148‑1 (MQ=255)
                      tttaTGCCCGCCGGGCTGCACTTTTTAACGCAGAAGTTGATGAACCTCGCCCTGGTTTACTTCGCGGACATTTTCCCGGAGCACTGAATTTTCCGTGGACGGAACTGGTGCTCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGa                                                                                                                                 <  2:6774547/147‑1 (MQ=255)
                      ttGAGGCCCGCCCGGCTGCACGTTTTAACGCAGACGTTGATGAACCTCGCCCAGGTTTACGTGGCGGACATATTCCCGGAGCACTTAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGa                                                                                                                                 <  2:1165654/150‑1 (MQ=255)
                             ccgcccgGCTGCACTTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGTACATATTCCCGGAGCACTGAATGTTCCGTGGACTTAACTGGTGCGCGAAGGCTAACTAAAAACGACCGATGAACTGGATGCGATAttttt                                                                                                                          <  2:7327999/150‑1 (MQ=255)
                             ccgcccgGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGTACATATTCCCGGAGCACTGAATGTTCCGTGGACGTAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATAttttt                                                                                                                          <  2:3704682/150‑1 (MQ=255)
                                     cTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTcgcg                                                                                                                  >  1:1288841/1‑150 (MQ=255)
                                      ggCACGTTTTACCGCTGATGTGGATGCCCGTCGCCCTGGTTTCGGTCGCGGACATTTGCCGGGTGCAGTGTTGGTTCGGTGGAGGGACCTGGTGGGCGAAGGCGAACTAAACAGGACCGATGAACTGGATGCGATATTTTTGGGTCgcgg                                                                                                                 <  2:5077173/149‑1 (MQ=255)
                                        cACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCgcggcg                                                                                                               >  2:5214365/1‑150 (MQ=255)
                                         aCGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGt                                                                                                              >  1:3758425/1‑150 (MQ=255)
                                             tttAACGCAGATGTTGATGATCCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGCGCCCTGAATGTTGCGTGGACGGAAGTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGc                                                                                                          <  2:6262732/150‑1 (MQ=255)
                                             tttAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGc                                                                                                          <  2:4840287/150‑1 (MQ=255)
                                             tttAACGAATAAGTTCATGACCCTCGCCCTGGTTTCCGTCGCGTAAATATCCCGGTAGAACTTTATGTTCGGTGTACGTAAGTGGTGCGCGAAGCCTAAATAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGc                                                                                                          <  1:5939782/150‑1 (MQ=255)
                                                aaCGCTGATGTGTATGAACCTCGCCCTGTTTTAGGTCGCGTACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAAGTGGTGCGCTAAGGCGAACTAAACACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTAc                                                                                                       <  2:3585430/150‑1 (MQ=255)
                                                  cGCAGATTTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATTTTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAAATAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGa                                                                                                     <  2:1460620/150‑1 (MQ=255)
                                                  cGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGa                                                                                                     <  2:711444/150‑1 (MQ=255)
                                                  cGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAAGGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAACACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGa                                                                                                     <  2:6895667/150‑1 (MQ=255)
                                                          ttGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATCTTTTTTGGTCGCGGCGTCAGCTACGACAAACCac                                                                                             >  1:7373597/1‑149 (MQ=255)
                                                          ttGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCaa                                                                                             >  2:5518640/1‑150 (MQ=255)
                                                          ttGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCaa                                                                                             <  2:7759233/150‑1 (MQ=255)
                                                           tgctgaACCTCGCCCTGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGTAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAAt                                                                                            <  2:1541737/147‑1 (MQ=255)
                                                            gatgaACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAAtt                                                                                           >  2:4799476/1‑150 (MQ=255)
                                                              tgaACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTAt                                                                                         >  2:2598379/1‑150 (MQ=255)
                                                               gaACCTCGCCCGGGTTTACGTGGGGGACTTTTTCCGGGTGCACTTATTGTTCGGTGGAGGGAACTGGGGGGGGACGGCTAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATc                                                                                        <  1:8418149/150‑1 (MQ=255)
                                                               gaACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGTAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATc                                                                                        <  1:5418322/150‑1 (MQ=255)
                                                                  ccTCGCCCTGGTTTCCGCCGCGTAAATATGCCCGGCGCACTTAATGTCCGGTGGCCGGAACTGGTGCGCGAAGGCGAAATAAACACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTc                                                                                     <  1:4031693/150‑1 (MQ=255)
                                                                  ccTCGCCCCGGTTTACGTCGCGTACATTTTCCCGGAGCACTGAATGTTCGGTGGACGGAAGTGGTGCGCGAAGGCTAACAAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTc                                                                                     <  2:3191310/150‑1 (MQ=255)
                                                                  ccTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAGGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTc                                                                                     >  1:8199006/1‑150 (MQ=255)
                                                                  ccTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTc                                                                                     <  2:8412006/150‑1 (MQ=255)
                                                                  ccTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTc                                                                                     >  1:251153/1‑150 (MQ=255)
                                                                   cTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAg                                                                                   >  1:1934422/1‑150 (MQ=255)
                                                                   cTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAg                                                                                   >  1:6167440/1‑150 (MQ=255)
                                                                       ccctgttttACGTCGCGCACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAAGTGGTGCGCTAAGGCTAAATAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTg                                                                                <  1:2427188/144‑1 (MQ=255)
                                                                       cccAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTg                                                                                >  2:827136/1‑150 (MQ=255)
                                                                       cccAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTg                                                                                >  2:4630123/1‑150 (MQ=255)
                                                                            gTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGct                                                                           >  1:4543657/1‑150 (MQ=255)
                                                                              ttACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAACACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGctct                                                                         <  2:1939668/150‑1 (MQ=255)
                                                                                aCGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTgg                                                                       >  1:711192/1‑150 (MQ=255)
                                                                                 cGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGgt                                                                      <  1:8069527/150‑1 (MQ=255)
                                                                                   tCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGgtgt                                                                    >  2:3999413/1‑150 (MQ=255)
                                                                                   tCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGgtgt                                                                    >  1:6646190/1‑150 (MQ=255)
                                                                                   tCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGAAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGgtgt                                                                    >  2:4691374/1‑150 (MQ=255)
                                                                                      cgGACATATTCCCGGAGCACTTATTGTTCCGTGGACGGAACTGGTGCGCGACGGCGACCAAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAAc                                                                 <  1:4204848/150‑1 (MQ=255)
                                                                                      cgGACATAGCCCGGGAGCATGGATGGTTCGGGGGCCGGAAGTGGTGCGCGACGGCTAACTAAACACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAAc                                                                 <  1:7913082/150‑1 (MQ=255)
                                                                                         aCATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGc                                                              >  1:3496500/1‑150 (MQ=255)
                                                                                              ttCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCg                                                         >  2:7642773/1‑150 (MQ=255)
                                                                                              ttCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATGGTTTGGTCGCGGCGTCAGCTACGACACACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCg                                                         >  2:6371345/1‑150 (MQ=255)
                                                                                                  cGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGt                                                     >  1:7722483/1‑150 (MQ=255)
                                                                                                   ggAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGtt                                                    >  1:3255189/1‑150 (MQ=255)
                                                                                                       cACTGATTGTTCGGTGGACGTAAGTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGt                                                <  1:7322573/150‑1 (MQ=255)
                                                                                                       cACTGAATGTTCCGTGGACGGAACTGGTGCGCTAAGGCTAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGt                                                <  1:6790450/150‑1 (MQ=255)
                                                                                                       cACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGt                                                <  2:1288841/150‑1 (MQ=255)
                                                                                                       cACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGt                                                >  1:3644963/1‑150 (MQ=255)
                                                                                                         cTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTa                                              <  2:813114/150‑1 (MQ=255)
                                                                                                         cTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTa                                              <  2:2963836/150‑1 (MQ=255)
                                                                                                          tGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAg                                             >  1:1421849/1‑150 (MQ=255)
                                                                                                          tGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAg                                             >  2:2300520/1‑150 (MQ=255)
                                                                                                          tGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAg                                             >  1:5878763/1‑150 (MQ=255)
                                                                                                          tGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAg                                             <  2:4541205/150‑1 (MQ=255)
                                                                                                          tGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAg                                             >  1:2113253/1‑150 (MQ=255)
                                                                                                             cagTTCCGTCCACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCAc                                          >  1:1502098/3‑150 (MQ=255)
                                                                                                                 tCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTcgc                                      >  2:6571637/1‑150 (MQ=255)
                                                                                                                   cGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGa                                    >  2:903495/1‑150 (MQ=255)
                                                                                                                    gTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGAc                                   >  1:1966532/1‑150 (MQ=255)
                                                                                                                     tGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACg                                  >  1:2907808/1‑150 (MQ=255)
                                                                                                                     tGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGAGGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACGCGCGACg                                  >  2:1371245/1‑150 (MQ=255)
                                                                                                                      gcaCGGAAGTGTGGCGCGACGGCTAAATAAACACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGc                                 <  1:2279966/148‑1 (MQ=255)
                                                                                                                       gACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCt                                <  2:3853204/150‑1 (MQ=255)
                                                                                                                        aCGGACCTGGTGCGCGAAGGCGAACTAAACACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTg                               <  2:4198266/150‑1 (MQ=255)
                                                                                                                        aCGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTg                               <  2:7527032/150‑1 (MQ=255)
                                                                                                                          ggAACGGGTGCGCGAAGGCGAACTAAACACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGa                             <  1:1086919/150‑1 (MQ=255)
                                                                                                                            aaCTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGAtg                           >  1:8031633/1‑150 (MQ=255)
                                                                                                                               tGGTGCGCGAAGGCTAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGc                        <  1:5015706/150‑1 (MQ=255)
                                                                                                                               tGGTGCGCGAAGGCGAACTAAACACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGc                        <  2:6074783/150‑1 (MQ=255)
                                                                                                                               tGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGc                        <  2:5884130/150‑1 (MQ=255)
                                                                                                                               tGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGc                        <  1:4594372/150‑1 (MQ=255)
                                                                                                                                   gcgcGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCaaa                    <  1:2567234/150‑1 (MQ=255)
                                                                                                                                      cgAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGt                 >  2:6988871/1‑150 (MQ=255)
                                                                                                                                          ggCGAACTAAACACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGaaa             <  1:183789/150‑1 (MQ=255)
                                                                                                                                          ggCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGaaa             <  1:6378550/150‑1 (MQ=255)
                                                                                                                                          ggCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGaaa             <  2:954566/150‑1 (MQ=255)
                                                                                                                                          ggCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGaaa             <  1:2667853/150‑1 (MQ=255)
                                                                                                                                          ggCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGaaa             <  1:4109619/150‑1 (MQ=255)
                                                                                                                                           gCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAAc            >  2:3404326/1‑150 (MQ=255)
                                                                                                                                             gAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTg          >  2:5647886/1‑150 (MQ=255)
                                                                                                                                                 tAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTacg      >  1:4832366/1‑150 (MQ=255)
                                                                                                                                                 tAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTGTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTacg      >  2:5805607/1‑150 (MQ=255)
                                                                                                                                                  aaacaCGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTacga     <  2:4444226/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGG  >  NC_000913/2652986‑2653284

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: