Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,370,155 A→G intergenic (+152/‑61) pspE → / → ycjM thiosulfate:cyanide sulfurtransferase (rhodanese)/alpha amylase catalytic domain family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,370,1550AG100.0% 655.4 / NA 168intergenic (+152/‑61)pspE/ycjMthiosulfate:cyanide sulfurtransferase (rhodanese)/alpha amylase catalytic domain family protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (91/77);  total (91/77)

Only 100 of 168 total aligned reads displayed.

ATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGAAACTTGTAACGCGTCTT  >  NC_000913/1370008‑1370302
                                                                                                                                                   |                                                                                                                                                   
aTTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAg                                                                                                                                                  <  2:930314/151‑1 (MQ=255)
  tCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGt                                                                                                                                                 <  2:1376172/150‑1 (MQ=255)
   cACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGttt                                                                                                                                               <  3:573773/151‑1 (MQ=255)
     ccATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTAt                                                                                                                                             <  2:533079/151‑1 (MQ=255)
      cATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTAtt                                                                                                                                            >  2:253337/1‑151 (MQ=255)
      cATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTAtt                                                                                                                                            <  2:1148603/151‑1 (MQ=255)
      cATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTAtt                                                                                                                                            <  1:222397/151‑1 (MQ=255)
        tGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTAtttt                                                                                                                                          >  4:561004/1‑151 (MQ=255)
         gAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTc                                                                                                                                         <  4:293900/151‑1 (MQ=255)
             ggcgCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCg                                                                                                                                     >  2:993429/1‑151 (MQ=255)
               cgcTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCgg                                                                                                                                    >  1:634764/1‑150 (MQ=255)
               cgcTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGa                                                                                                                                   <  2:1139184/151‑1 (MQ=255)
                gcTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGag                                                                                                                                  <  1:1083414/151‑1 (MQ=255)
                  ttATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGagt                                                                                                                                 >  2:120017/1‑150 (MQ=255)
                   tATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGagt                                                                                                                                 <  1:932317/149‑1 (MQ=255)
                     tGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGtc                                                                                                                             >  3:66411/1‑151 (MQ=255)
                       cgcCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGtct                                                                                                                            <  2:1311368/150‑1 (MQ=255)
                         ccGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGtctctc                                                                                                                         >  3:95187/1‑151 (MQ=255)
                          cGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCg                                                                                                                        >  2:49441/1‑151 (MQ=255)
                           ttttttttttCTGTCACACCTTATTTACATCCCCATAGATTATTTTCTTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGt                                                                                                                       >  1:448375/2‑151 (MQ=255)
                           gTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGt                                                                                                                       >  3:291300/1‑151 (MQ=255)
                             ttttttttCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGtt                                                                                                                      <  4:3936/150‑1 (MQ=255)
                               ttttttCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAATAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTc                                                                                                                    >  2:371179/1‑150 (MQ=255)
                               ttttttCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTc                                                                                                                    <  1:793791/150‑1 (MQ=255)
                               ttttttCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTc                                                                                                                    <  2:1451643/150‑1 (MQ=255)
                               ttttttCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTc                                                                                                                    >  2:23432/1‑150 (MQ=255)
                                  tttCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCTCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATg                                                                                                                 >  2:257537/1‑150 (MQ=255)
                                   ttCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATggg                                                                                                               >  1:1598351/1‑151 (MQ=255)
                                        tCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCg                                                                                                           >  1:1162009/1‑150 (MQ=255)
                                         cacaCCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCt                                                                                                         <  2:1738817/151‑1 (MQ=255)
                                             ccTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTAcc                                                                                                      <  2:1246584/150‑1 (MQ=255)
                                                tATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACgg                                                                                                  >  4:466522/1‑151 (MQ=255)
                                                tATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACg                                                                                                   <  2:72455/150‑1 (MQ=255)
                                                 aTTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGa                                                                                                 <  4:587170/151‑1 (MQ=255)
                                                  tttACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGa                                                                                                 >  3:266756/1‑150 (MQ=255)
                                                     aCATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAgg                                                                                             <  2:993797/151‑1 (MQ=255)
                                                     aCATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAg                                                                                              <  1:877888/150‑1 (MQ=255)
                                                      cacccccATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGaaa                                                                                               <  2:28885/145‑1 (MQ=255)
                                                      cATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGc                                                                                            >  2:926008/1‑151 (MQ=255)
                                                       aTCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGc                                                                                            >  4:429316/1‑150 (MQ=255)
                                                       aTCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGc                                                                                            >  2:680117/1‑150 (MQ=255)
                                                           ccATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACAt                                                                                       <  4:264926/151‑1 (MQ=255)
                                                            cATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATg                                                                                      <  1:1333341/151‑1 (MQ=255)
                                                            cATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATg                                                                                      >  3:232026/1‑151 (MQ=255)
                                                               aGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGAAACATgca                                                                                    >  1:400098/1‑148 (MQ=255)
                                                                gATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAAc                                                                                  >  1:1320912/1‑151 (MQ=255)
                                                                   tattTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAg                                                                                <  1:253294/150‑1 (MQ=255)
                                                                     tttGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGaa                                                                              >  4:230993/1‑150 (MQ=255)
                                                                        gCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAAtt                                                                          >  2:1335128/1‑151 (MQ=255)
                                                                        gCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAAt                                                                           <  1:73278/150‑1 (MQ=255)
                                                                         cGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAAtt                                                                          >  2:943738/1‑150 (MQ=255)
                                                                         cGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTa                                                                         <  1:346905/151‑1 (MQ=255)
                                                                          gTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTAc                                                                        <  1:1069865/151‑1 (MQ=255)
                                                                            cAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTAc                                                                        <  1:893539/149‑1 (MQ=255)
                                                                            cAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACgg                                                                      <  3:545687/151‑1 (MQ=255)
                                                                             aGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACgg                                                                      >  2:1293699/1‑150 (MQ=255)
                                                                                   caAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCt                                                               >  2:1015822/1‑151 (MQ=255)
                                                                                     aaTACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACTGAAAGGCAACATGAAACAGAAAATTACGGATTACCTgg                                                             <  2:1051110/151‑1 (MQ=255)
                                                                                      aTACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTgg                                                             >  2:181465/1‑150 (MQ=255)
                                                                                        aCGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACg                                                          >  2:504746/1‑151 (MQ=255)
                                                                                             ttttCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAAt                                                      <  1:120573/150‑1 (MQ=255)
                                                                                                tCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTa                                                   <  2:116708/150‑1 (MQ=255)
                                                                                                tCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTa                                                   <  1:32122/150‑1 (MQ=255)
                                                                                                tCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTa                                                   <  3:295149/150‑1 (MQ=255)
                                                                                                tCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTAc                                                  >  2:1591220/1‑151 (MQ=255)
                                                                                                 cccTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTAc                                                  <  2:1592921/150‑1 (MQ=255)
                                                                                                 cccTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACg                                                 >  1:1513156/1‑151 (MQ=255)
                                                                                                     ggTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTgg                                             >  2:1373358/1‑151 (MQ=255)
                                                                                                      gTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGa                                            <  3:129956/151‑1 (MQ=255)
                                                                                                        aaacactGATTTCCTGCGTGACTACAACCTCTGAGCTCTCTTGGCATTTTAGTTTCTCCGGCGTAGTTTCACGTTTCATGGGACCGCTACCACGGTAAGTCACTATGAAACTGAAAATACCGGATTACATGAACGAACTCAACGGCGGAAc                                          >  1:1262978/7‑151 (MQ=255)
                                                                                                        aaaaaaTGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAAc                                          >  1:660390/1‑151 (MQ=255)
                                                                                                          aaatTGATTTCCTGCGTGGCTAAAACCCATGGGCGCAATTGGCAGTTTATTTTCTGCGGAGCAGTCTCTCGTTTAATGAGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACAt                                        >  1:709129/1‑151 (MQ=255)
                                                                                                          aaaaTGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACAt                                        >  1:1446424/1‑151 (MQ=255)
                                                                                                           aaaTGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACAtt                                       <  1:940006/151‑1 (MQ=255)
                                                                                                           aaaTGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACAtt                                       >  2:842052/1‑151 (MQ=255)
                                                                                                            aaTGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACAttt                                      >  2:591467/1‑151 (MQ=255)
                                                                                                               gATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTAcc                                   <  1:634123/151‑1 (MQ=255)
                                                                                                                    ccTGCGTGACTAAAACCCTTGTGCTCACTTGGCAGTTTATTTTCTGCGGAGCAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGTCGACATCTACGGTGGAACATTTACCGCAAc                              >  1:181093/1‑151 (MQ=255)
                                                                                                                    ccTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCaa                               <  2:868553/150‑1 (MQ=255)
                                                                                                                    ccTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCaa                               <  2:814654/150‑1 (MQ=255)
                                                                                                                      tGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTc                            >  2:1704278/1‑151 (MQ=255)
                                                                                                                       gCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTc                            >  2:521101/1‑150 (MQ=255)
                                                                                                                       gCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTc                            >  1:416950/1‑150 (MQ=255)
                                                                                                                       gCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACCAAATCTACGGTGGAACATTTACCGCAACt                             <  4:474498/149‑1 (MQ=255)
                                                                                                                           gACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCAttt                        <  2:810099/150‑1 (MQ=255)
                                                                                                                           gACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTa                       <  2:742837/151‑1 (MQ=255)
                                                                                                                             cTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTAc                      >  1:1691952/1‑150 (MQ=255)
                                                                                                                              tAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTAc                      >  1:857158/1‑149 (MQ=255)
                                                                                                                               aaaaCCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTAc                      >  1:1339626/1‑148 (MQ=255)
                                                                                                                               aaaaCCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTAc                      >  3:203491/1‑148 (MQ=255)
                                                                                                                               aaaaCCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGa                   <  1:752442/151‑1 (MQ=255)
                                                                                                                                   cccTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGAAACt               >  2:375291/1‑151 (MQ=255)
                                                                                                                                     cTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGAAACTTg             >  1:1662367/1‑151 (MQ=255)
                                                                                                                                       tgtgCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTGCCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGAAACTTg             <  4:7937/149‑1 (MQ=255)
                                                                                                                                       tgtgCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGAAACTTg             <  1:778103/149‑1 (MQ=255)
                                                                                                                                       tgtgCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGAAACTTGt            <  2:784537/150‑1 (MQ=255)
                                                                                                                                       tgtgCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGAAACTTGt            <  2:446869/150‑1 (MQ=255)
                                                                                                                                        gtgCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGAAACTTGTaa          <  1:800820/151‑1 (MQ=255)
                                                                                                                                             cAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGAAACTTGTAAcgcg      <  2:1513316/150‑1 (MQ=255)
                                                                                                                                                ttGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGAAACTTGTAACGCGTCtt  >  1:98664/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGAAACTTGTAACGCGTCTT  >  NC_000913/1370008‑1370302

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