Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,533,415 T→A 12.6% intergenic (+35/‑349) cysK → / → ptsH cysteine synthase A, O‑acetylserine sulfhydrolase A subunit/phosphohistidinoprotein‑hexose phosphotransferase component of PTS system (Hpr)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,533,4150TA12.6% 195.0 / 15.0 79intergenic (+35/‑349)cysK/ptsHcysteine synthase A, O‑acetylserine sulfhydrolase A subunit/phosphohistidinoprotein‑hexose phosphotransferase component of PTS system (Hpr)
Reads supporting (aligned to +/- strand):  ref base T (28/41);  new base A (4/6);  total (32/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

TGAAAGCTTTACCAACAAGAATATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATG  >  NC_000913/2533269‑2533561
                                                                                                                                                  |                                                                                                                                                  
tGAAATCTTTACCTACAAGAATATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGTTCTCTACACTGAGAAAGAATTGCATCAGTATTGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCttttt                                                                                                                                                 <  2:325601/150‑1 (MQ=255)
tGAAAGCTTTACCAACAAGAATATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCttttt                                                                                                                                                 >  1:929483/1‑150 (MQ=255)
 gAAAGCTTTACCAACAAGAATATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTa                                                                                                                                                <  1:188747/150‑1 (MQ=255)
      cTTTACCAACAAGAATATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGc                                                                                                                                           >  1:682269/1‑150 (MQ=255)
       tttACCAACAAGAATATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCt                                                                                                                                          >  1:718309/1‑150 (MQ=255)
          aCCAACAAGAATATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCtttt                                                                                                                                       <  2:897881/150‑1 (MQ=255)
           ccAACAAGAATATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTAAAAAGGTGCttttt                                                                                                                                      <  2:587626/150‑1 (MQ=255)
            caacaaGAATATTGTGGTTATTCTACCGTCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCtttttt                                                                                                                                     >  2:141348/1‑150 (MQ=255)
              acaaGAATATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGAATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTAAAAAGGTGCttttttac                                                                                                                                   <  1:882775/150‑3 (MQ=255)
               caaGAATTTTGTGGTTATTCTACCATCATCGGGTGAGCGATATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTAAAAAGGTGCttttttacg                                                                                                                                  <  2:694046/150‑4 (MQ=255)
               caaGAATATTGTTGTTATTCTACCATCATCGGGTGATCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTAAAAAGGTGCttttttacg                                                                                                                                  <  2:961450/150‑4 (MQ=255)
               caaGAATATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTtgtg                                                                                                                                  >  2:607205/1‑150 (MQ=255)
               caaGAATATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTAAAAAGGTGCttttttacg                                                                                                                                  <  2:1460015/150‑4 (MQ=255)
                 aGAATATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCAGATATCTTCACTGAGAAAGAATTTCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAATCACCTTTTTAGGTGCTTTTTTGTGGc                                                                                                                                >  1:296143/1‑150 (MQ=255)
                  gAATATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGcc                                                                                                                               <  1:706512/150‑1 (MQ=255)
                  caaTATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGcc                                                                                                                               <  2:970611/149‑1 (MQ=255)
                     tatTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGATGCTTTTTTGTGGCCTGc                                                                                                                            <  1:1248571/150‑1 (MQ=255)
                      atTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCt                                                                                                                           <  2:1708224/150‑1 (MQ=255)
                        tgtgGTTATTCTACCATCATCGGGTGAGCTTGATTAAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATTCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTc                                                                                                                         <  1:438437/150‑1 (MQ=255)
                         gtgGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAGAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCa                                                                                                                        >  2:672954/1‑150 (MQ=255)
                         gtgGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCa                                                                                                                        <  1:1260230/150‑1 (MQ=255)
                          tgGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCaa                                                                                                                       <  2:316775/150‑1 (MQ=255)
                          tgGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCaa                                                                                                                       <  1:545477/150‑1 (MQ=255)
                          tgGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCaa                                                                                                                       <  2:685288/150‑1 (MQ=255)
                          tgGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCaa                                                                                                                       >  2:306759/1‑150 (MQ=255)
                           ggTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATATCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGAGTAAAAAAGCACCTTTTTAGGAGCTTTTTTGTGGCATGCTTCaaa                                                                                                                      >  2:1183086/1‑150 (MQ=255)
                            ggtATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCtt                                                                                                                          >  1:352868/3‑145 (MQ=255)
                            ggtATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCtt                                                                                                                          <  2:352868/143‑1 (MQ=255)
                            gTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAAc                                                                                                                     <  1:306759/150‑1 (MQ=255)
                            gTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAAc                                                                                                                     <  1:718692/150‑1 (MQ=255)
                               aTTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAGAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACttt                                                                                                                  <  1:672954/150‑1 (MQ=255)
                               aTTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATTCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACttt                                                                                                                  <  2:1519698/150‑1 (MQ=255)
                                     ccATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGcccc                                                                                                            <  2:1682144/150‑1 (MQ=255)
                                        tcatcGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCcctcc                                                                                                         <  2:1654085/150‑1 (MQ=255)
                                         catcGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCcctcct                                                                                                        <  2:446543/150‑1 (MQ=255)
                                         catcGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTACAAACTTTCGCCcctcca                                                                                                        >  2:486258/1‑149 (MQ=255)
                                           tcGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTgg                                                                                                      <  1:1187422/150‑1 (MQ=255)
                                                 gAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGa                                                                                                <  1:607205/150‑1 (MQ=255)
                                                  aGCGTTATTTAAGCACCGCATGGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGAt                                                                                               <  2:171194/150‑1 (MQ=255)
                                                    cGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGAtta                                                                                             >  1:134078/1‑149 (MQ=255)
                                                    cGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTc                                                                                             >  2:455749/1‑150 (MQ=255)
                                                      ttATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAg                                                                                           >  2:1647452/1‑150 (MQ=255)
                                                          ttAAGAACCGCATTGTTTGCCGAACTCTTCACAGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCATTTTAGTTGATTTTTTGTGGCCTGCTTCAAACTTTCGCAACTCCTGGCATTTATTAAGCCTg                                                                                       >  2:1390468/1‑150 (MQ=255)
                                                              gCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCgg                                                                                   <  2:407874/150‑1 (MQ=255)
                                                               cACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGGTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGa                                                                                  <  2:1558688/150‑1 (MQ=255)
                                                                     aTTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTAGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGTCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGt                                                                            <  2:307115/150‑1 (MQ=255)
                                                                      ttGTTTGCCGATGTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTa                                                                           <  2:296143/150‑1 (MQ=255)
                                                                      ttGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTa                                                                           <  1:1357967/150‑1 (MQ=255)
                                                                         tttGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTATTTGTGGCCTGCTTCAAACGTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGGAtgt                                                                        >  1:1617067/1‑148 (MQ=255)
                                                                            gCCGATCTCTTCCCGGAGACAGACTTGCCACAGTAATGACCGCTGGTTAAAAATGCGTAAAACAGCACCTTTTTCGGTGCTTTTTTGTGGCCTGCTTCAACCTTTCGCCCCTCCTGGCATTGATTCCGCCTGTCGGACCTGGTATTTAAc                                                                     <  2:1617067/150‑1 (MQ=255)
                                                                            gCCGATCTCTTCACTGAGAAATAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAAc                                                                     <  1:1311941/150‑1 (MQ=255)
                                                                               gATCTCTTCACTTAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAg                                                                  <  2:102881/150‑1 (MQ=255)
                                                                               gATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAg                                                                  <  2:459267/150‑1 (MQ=255)
                                                                                    ctTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCACCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTaa                                                             >  2:936749/1‑150 (MQ=255)
                                                                                      tCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAatt                                                           <  2:718309/150‑1 (MQ=255)
                                                                                           gagaAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTAtttt                                                      >  1:1048717/1‑150 (MQ=255)
                                                                                            tgtaagaaTTGCAACAGTAATGCCAGATTGTTAAAAATGAGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCTTTGATTAAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTg                                                     <  2:1603254/147‑1 (MQ=255)
                                                                                            agaaagaaTTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTg                                                     >  2:1292796/1‑150 (MQ=255)
                                                                                             gaaagaaTTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGa                                                    >  2:52228/1‑150 (MQ=255)
                                                                                              aaagaaTTGCAACTGTTATTCCAGCTTGTTATAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCGCCTCCTGGCTTTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGAt                                                   <  1:1183086/150‑1 (MQ=255)
                                                                                                   aTTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGcgcg                                              <  1:1614052/150‑1 (MQ=255)
                                                                                                    ttGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGa                                             >  1:1394979/1‑150 (MQ=255)
                                                                                                      gCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTGTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGaaa                                           >  1:159987/1‑150 (MQ=255)
                                                                                                       cAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAAt                                          <  1:199571/150‑1 (MQ=255)
                                                                                                       cAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAAt                                          <  1:1122721/150‑1 (MQ=255)
                                                                                                         aCAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCTCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACACGACTAATTATTTTGTTGCGCGAAATTa                                        >  1:872176/1‑150 (MQ=255)
                                                                                                         aCAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTa                                        >  1:1594533/1‑150 (MQ=255)
                                                                                                              aaTGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGt                                   <  1:526234/150‑1 (MQ=255)
                                                                                                              aaTGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGt                                   <  2:1709816/150‑1 (MQ=255)
                                                                                                               aTGCCAGCTTGTTAAAAATGCGTAATAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGtt                                  <  1:1267048/150‑1 (MQ=255)
                                                                                                                      cTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGaa                           >  1:680819/1‑150 (MQ=255)
                                                                                                                       ttGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGaaa                          >  1:1207226/1‑150 (MQ=255)
                                                                                                                          ttAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGc                       >  2:1618375/1‑150 (MQ=255)
                                                                                                                            aaaaaTGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCa                     <  1:905197/150‑1 (MQ=255)
                                                                                                                                     caaaaaaGCACCTAAAAAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAAt            >  1:836338/2‑150 (MQ=255)
                                                                                                                                     caaaaaaGCACCTAAAAAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAAt            >  1:823878/2‑150 (MQ=255)
                                                                                                                                       aaaaaGCACCTAAAAAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATct          >  2:317423/1‑149 (MQ=255)
                                                                                                                                        aaaaGCACCTAAAAAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGa         >  2:685983/1‑150 (MQ=255)
                                                                                                                                         aaaGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGAt        >  1:1706576/1‑150 (MQ=255)
                                                                                                                                         aaaGCACCTAAAAAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGAt        <  2:836338/150‑1 (MQ=255)
                                                                                                                                            gCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGAtttt     >  1:244329/1‑150 (MQ=255)
                                                                                                                                               cgtttttAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGTTTCAGCGTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCATAGCTGAATCGATTTTATg  <  2:1354376/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TGAAAGCTTTACCAACAAGAATATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATG  >  NC_000913/2533269‑2533561

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: