Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,538,557 T→C 5.8% intergenic (+74/+115) yfeS → / ← cysM WGR domain protein/cysteine synthase B (O‑acetylserine sulfhydrolase B)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,538,5570TC 5.8% 360.5 / 8.0 119intergenic (+74/+115)yfeS/cysMWGR domain protein/cysteine synthase B (O‑acetylserine sulfhydrolase B)
Reads supporting (aligned to +/- strand):  ref base T (49/63);  new base C (4/3);  total (53/66)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.64e-01

Only 100 of 121 total aligned reads displayed.

ATGGGGTGACGGTCAGGACTCACCAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCC  >  NC_000913/2538408‑2538706
                                                                                                                                                     |                                                                                                                                                     
aTGGGGTGACGGTCAGGACTCACCAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCtt                                                                                                                                                       >  1:511757/1‑150 (MQ=255)
      tGACGGTCAGGACTCACCAATTCTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCTCGAATATTAAATAGTGCATTAGTCGGCAGAAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTTACGCCGCATCCgg                                                                                                                                                 <  2:1059594/150‑1 (MQ=38)
        aCGGTCAGGACTCACCAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCg                                                                                                                                                  >  1:265967/1‑147 (MQ=255)
        aCGGTCAGGACTCACCAATACATCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                                                                                               <  2:1193291/150‑1 (MQ=255)
          ggTCAGGACTCACCAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                                                                                                             >  2:1099988/1‑150 (MQ=255)
          ggTCAGGACTCACCAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGAGGCGTAAACGCCTTATCCGGCCTa                                                                                                                                             >  2:1328208/1‑150 (MQ=255)
           gTCACGACTCACCAATACATCAAAAACTGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCCGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTTACGCCGCATCCGGCATaa                                                                                                                                            <  1:305066/150‑2 (MQ=255)
                gACTCACCAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAt                                                                                                                                       <  1:572643/150‑1 (MQ=255)
                   tCACCAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGt                                                                                                                                    <  1:1473097/150‑1 (MQ=255)
                    cACCAATACTTCAAAAAATGATTGCTGACCTTACGGCGTTTCCTCACGAAAATTAATTACAGCATTTGTCGGCAGCAACACCTGTTAAAAATGGGCCCTTTGTATATTCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTg                                                                                                                                   >  1:78097/1‑150 (MQ=37)
                     aCCAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATTCAGCATTTGTCGGCAGCATCAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGc                                                                                                                                  <  2:182535/150‑1 (MQ=255)
                     aCCAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTACGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGc                                                                                                                                  <  1:788208/150‑1 (MQ=255)
                      ccAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCt                                                                                                                                 >  1:356954/1‑150 (MQ=255)
                      ccAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCt                                                                                                                                 <  2:603778/150‑1 (MQ=255)
                       cAATACTTCAAAAACTGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTa                                                                                                                                <  1:282179/150‑1 (MQ=255)
                           aCTTCAAAAAATGATTGATGACCTTACGTCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTc                                                                                                                            <  2:481101/150‑1 (MQ=255)
                           aCTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTc                                                                                                                            <  1:279264/150‑1 (MQ=255)
                             ttCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGAAGCGGCGTAATCGCCTTATCCGGCCTACATAATCGTGCTAATTCaa                                                                                                                          >  1:88189/1‑150 (MQ=255)
                              tCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAat                                                                                                                         <  2:70521/150‑1 (MQ=255)
                                    aaTGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTg                                                                                                                   <  2:1529203/150‑1 (MQ=255)
                                     aTGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGc                                                                                                                  <  2:616815/150‑1 (MQ=255)
                                     aTGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGc                                                                                                                  >  1:393165/1‑150 (MQ=255)
                                       gattgatGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAg                                                                                                                <  2:356954/150‑1 (MQ=255)
                                           gatgaCCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAc                                                                                                            >  1:156061/1‑150 (MQ=255)
                                           gatgaCCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAc                                                                                                            >  1:621712/1‑150 (MQ=255)
                                            atgaCCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAcc                                                                                                           <  1:801435/150‑1 (MQ=255)
                                             tgaCCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCt                                                                                                          >  1:397647/1‑150 (MQ=255)
                                                 cTTACGGCGTTTCCTCACGAAAATTAAATACAGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTa                                                                                                      <  2:1441027/150‑1 (MQ=255)
                                                     cGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACCACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGcc                                                                                                  >  1:230492/1‑150 (MQ=255)
                                                       gCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGCAACCTTGTAGGCCTg                                                                                                >  2:887400/1‑150 (MQ=255)
                                                          tttCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATa                                                                                             <  2:1495378/150‑1 (MQ=255)
                                                            tCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAg                                                                                           <  2:1679356/150‑1 (MQ=255)
                                                              cTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCg                                                                                         >  2:1304416/1‑150 (MQ=255)
                                                                 aCGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAATAAAGTTCGCGTTGTTTATGCCGTATGCGGAGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAg                                                                                      <  2:1178285/150‑1 (MQ=25)
                                                                 aCGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAg                                                                                      >  1:1046364/1‑150 (MQ=255)
                                                                 aCGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCCGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAg                                                                                      >  1:556451/1‑150 (MQ=255)
                                                                  cGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAgc                                                                                     >  1:101590/1‑150 (MQ=255)
                                                                  cGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAgc                                                                                     <  1:1063317/150‑1 (MQ=255)
                                                                   gAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGACGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAgcg                                                                                    <  1:1016235/150‑1 (MQ=255)
                                                                              aCTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGt                                                                         >  1:571611/1‑150 (MQ=255)
                                                                               cTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGtt                                                                        >  1:630857/1‑150 (MQ=255)
                                                                                tGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGttt                                                                       <  1:727356/150‑1 (MQ=255)
                                                                                 gCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGtttt                                                                      >  1:1494265/1‑150 (MQ=255)
                                                                                   aTTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGc                                                                    <  2:45046/150‑1 (MQ=255)
                                                                                    tttGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCg                                                                   <  1:1620941/150‑1 (MQ=255)
                                                                                    tttGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCg                                                                   <  1:1614530/150‑1 (MQ=255)
                                                                                     ttGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGt                                                                  <  1:1366203/150‑1 (MQ=255)
                                                                                       gTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGttt                                                                >  2:1143659/1‑150 (MQ=255)
                                                                                       gTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGttt                                                                >  1:672908/1‑150 (MQ=255)
                                                                                        gcGGCAGCAACAACTGTTAAAAACGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTg                                                               <  2:1307579/149‑1 (MQ=255)
                                                                                         gggCATCAACTACTGTTAAAAAAGTGCGCTTTGTTTATGCGGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTTGCGTTTGt                                                              <  1:237225/149‑1 (MQ=255)
                                                                                          ggCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGtc                                                             >  2:1105299/1‑150 (MQ=255)
                                                                                          ggCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGAGGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGtc                                                             <  1:1296811/150‑1 (MQ=255)
                                                                                           gcagcaACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGtca                                                            <  1:1246686/150‑1 (MQ=255)
                                                                                            cagcaACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGtcat                                                           <  2:1576132/150‑1 (MQ=255)
                                                                                               caacaaCTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAg                                                        <  2:397647/150‑1 (MQ=255)
                                                                                                 acaacTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGtc                                                      <  2:11805/150‑1 (MQ=255)
                                                                                                 acaacTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGta                                                      >  2:91804/1‑149 (MQ=255)
                                                                                                    acTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAACCGCCTTATCCGGCCTACAAAATCGTGCTAATTCAAAATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCTTTTGTCATCAGTCTcc                                                   >  1:382670/1‑150 (MQ=255)
                                                                                                    acTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTATTTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTcc                                                   >  2:205790/1‑150 (MQ=255)
                                                                                                      tGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGa                                                 <  1:1105299/150‑1 (MQ=255)
                                                                                                      tGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCCGa                                                 >  2:122955/1‑150 (MQ=255)
                                                                                                       tttAAAAAAGTGCGCTTTGTTTATGCCGGTTGCGGCGTAAACGCCTTATCCGGCCTACAAATTCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGAt                                                <  1:1004556/149‑1 (MQ=255)
                                                                                                           aaaaaGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTa                                            >  2:405290/1‑150 (MQ=255)
                                                                                                           aaaaaGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTa                                            <  2:1610911/150‑1 (MQ=255)
                                                                                                           aaaaaGTGAGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCTATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGcga                                            <  1:592013/150‑3 (MQ=255)
                                                                                                            aaaaGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTAt                                           <  1:1572738/150‑1 (MQ=255)
                                                                                                             aaaGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTAtt                                          <  1:1396378/150‑1 (MQ=255)
                                                                                                                gTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAAt                                       <  1:1255615/150‑1 (MQ=255)
                                                                                                                 tGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATc                                      <  1:1336367/150‑1 (MQ=255)
                                                                                                                 tGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATc                                      >  1:65307/1‑150 (MQ=255)
                                                                                                                     cTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATACAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTa                                  <  1:1337921/150‑1 (MQ=255)
                                                                                                                     cTTTGTTTATGCCAGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGAGTTTGTCATCAGTCTCCGATGCTATTAATCCTTa                                  >  2:547887/1‑150 (MQ=255)
                                                                                                                        tGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAAt                               <  2:1698342/150‑1 (MQ=255)
                                                                                                                           gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATccc                            >  1:997597/5‑150 (MQ=21)
                                                                                                                           gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGAGGCTATTAATCCTTAAATccc                            >  1:747649/5‑150 (MQ=18)
                                                                                                                             atcccGGATGCGGCGTAAACGCCGTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCg                          <  2:350617/147‑1 (MQ=255)
                                                                                                                             aTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCg                          <  2:1471902/150‑1 (MQ=255)
                                                                                                                             aTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCg                          <  2:393165/150‑1 (MQ=255)
                                                                                                                             aTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCg                          >  2:1247798/1‑150 (MQ=255)
                                                                                                                              tGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGc                         <  2:630857/150‑1 (MQ=255)
                                                                                                                              ggCCGGATAAGGCGTTTACGCCGCATCCGGCATACAAAATCGTGCTAATTCAATATATTGCATAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGc                         >  2:338049/2‑150 (MQ=21)
                                                                                                                                     tGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTgg                  <  2:1547678/150‑1 (MQ=255)
                                                                                                                                      gcggcgTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGc                 <  2:1540468/150‑1 (MQ=255)
                                                                                                                                        ggcgTAAACTCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTGATCCTTAATTCCCCGCCCCCTGGCTa               >  2:73390/1‑150 (MQ=255)
                                                                                                                                         gcgTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTaa              <  1:1499350/150‑1 (MQ=255)
                                                                                                                                         gcgTAAACGCCTTATCCGGCCTACAAAATCGGGCTAATTCAATATATTGCAGACACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCtca              >  1:136812/1‑148 (MQ=255)
                                                                                                                                          cgTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTaaa             >  2:15515/1‑150 (MQ=255)
                                                                                                                                          cgTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTaaa             <  2:1223109/150‑1 (MQ=255)
                                                                                                                                            tttaCGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGc                         >  1:1550120/4‑136 (MQ=255)
                                                                                                                                            tttaCGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGc                         <  2:1550120/133‑1 (MQ=255)
                                                                                                                                              aaCGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGc         <  2:195942/150‑1 (MQ=255)
                                                                                                                                              aaCGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGc         <  1:1581985/150‑1 (MQ=255)
                                                                                                                                              aaCGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGc         >  1:1276771/1‑150 (MQ=255)
                                                                                                                                                cGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGctc       <  2:671266/150‑1 (MQ=255)
                                                                                                                                                cGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGctc       >  2:302628/1‑150 (MQ=255)
                                                                                                                                                  ccTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGTCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCACCTGGCTAAAATGATCtt     >  2:180441/1‑150 (MQ=255)
                                                                                                                                                  ccTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCtt     >  2:1677098/1‑150 (MQ=255)
                                                                                                                                                  ccTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGACCCCTGGCTAAAATGCTCtt     >  2:294184/1‑150 (MQ=255)
                                                                                                                                                     tATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTccc  >  2:1199122/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATGGGGTGACGGTCAGGACTCACCAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCC  >  NC_000913/2538408‑2538706

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