Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,533,415 T→A 9.5% 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,4150TA 9.5% 236.3 / 14.2 74intergenic (+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 (34/33);  new base A (4/3);  total (38/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.92e-01

TGAAAGCTTTACCAACAAGAATATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGA  >  NC_000913/2533269‑2533562
                                                                                                                                                  |                                                                                                                                                   
tGAAAGCTTTACCAACAAGAATATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCttt                                                                                                                                                    <  2:209242/148‑1 (MQ=255)
 gAAAGCTTTACCAACAAGAATATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCttttt                                                                                                                                                  <  3:646005/149‑1 (MQ=255)
 gAAAGCTTTACCAACAAGAATATCGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCtttt                                                                                                                                                   <  1:399727/148‑1 (MQ=255)
  aaaGCTTTACCAACAAGAATATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCttttt                                                                                                                                                  >  2:638872/1‑148 (MQ=255)
  aaaGCTTTACCAACAAGAATATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCttttt                                                                                                                                                  >  1:374781/1‑148 (MQ=255)
   aaGCTTTACCAACAAGAATATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAg                                                                                                                                                >  4:335625/1‑149 (MQ=255)
    aGCTTTACCAACAAGAATATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAg                                                                                                                                                >  2:734527/1‑148 (MQ=255)
     gCTTTACCAACAAGAATATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAgg                                                                                                                                               <  1:412577/148‑1 (MQ=255)
           ccAACAAGAATATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGTTGCttt                                                                                                                                         >  1:644610/1‑148 (MQ=255)
           ccAACAAGAATATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTAAAAAGGTGCtttt                                                                                                                                        <  3:327748/149‑1 (MQ=255)
           ccAACAAGAATATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTAAAAAGGTGCttt                                                                                                                                         <  2:786714/148‑1 (MQ=255)
           ccAACAAGAATATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTAAAAAGGTGCttt                                                                                                                                         <  1:645406/148‑1 (MQ=255)
              acaaGAATATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTa                                                                                                                                                 >  3:775061/1‑137 (MQ=255)
              acaaGAATATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTa                                                                                                                                                 <  4:775061/137‑1 (MQ=255)
              acaaGAATATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTtg                                                                                                                                     >  4:366102/1‑149 (MQ=255)
               caaGAATATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTtg                                                                                                                                     <  2:543358/148‑1 (MQ=255)
               caaGAATATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTAACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTtgt                                                                                                                                    <  3:508512/149‑1 (MQ=255)
                 aGAATATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTtgtg                                                                                                                                   >  1:462175/1‑148 (MQ=255)
                     tatTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTg                                                                                                                              >  3:180097/1‑149 (MQ=255)
                            gTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCaa                                                                                                                        <  2:734984/148‑1 (MQ=255)
                            gTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCa                                                                                                                         <  4:374160/147‑1 (MQ=255)
                            gTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCa                                                                                                                         >  3:374160/1‑147 (MQ=255)
                                 tCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACttt                                                                                                                   <  1:638872/148‑1 (MQ=255)
                                             gggTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTgg                                                                                                       <  2:462175/148‑1 (MQ=255)
                                                  aGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCTAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTg                                                                                                  <  2:644610/148‑1 (MQ=255)
                                                     gTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTc                                                                                              <  3:316864/149‑1 (MQ=255)
                                                      ttATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCa                                                                                             <  3:769304/149‑1 (MQ=255)
                                                          ttAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAGTTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCt                                                                                                                                           <  2:92576/99‑1 (MQ=255)
                                                          ttAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCt                                                                                                                                           >  1:92576/1‑99 (MQ=255)
                                                          ttAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTTATTCAGcc                                                                                          >  1:396365/1‑148 (MQ=255)
                                                             aGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGt                                                                                       >  1:352579/1‑148 (MQ=255)
                                                                  cGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAAc                                                                                 <  3:789587/149‑1 (MQ=255)
                                                                  cGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAAc                                                                                 >  3:810102/1‑149 (MQ=255)
                                                                       tGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGt                                                                             <  2:396365/148‑1 (MQ=255)
                                                                       tGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTa                                                                            >  3:341882/1‑149 (MQ=255)
                                                                       tGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTa                                                                            <  3:802625/149‑1 (MQ=255)
                                                                        gTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTa                                                                            <  1:332623/148‑1 (MQ=255)
                                                                            gCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACttt                                                                                                                   <  4:852046/105‑1 (MQ=255)
                                                                            gCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACttt                                                                                                                   >  3:852046/1‑105 (MQ=255)
                                                                             ccGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTAttt                                                                         >  3:484723/1‑146 (MQ=255)
                                                                             ccGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTAttt                                                                         <  4:484723/146‑1 (MQ=255)
                                                                             ccGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTaa                                                                       <  2:374781/148‑1 (MQ=255)
                                                                              cGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAAc                                                                      >  2:99473/1‑148 (MQ=255)
                                                                                aTCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCa                                                                    >  2:561290/1‑148 (MQ=255)
                                                                                      tCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTaa                                                              >  1:540997/1‑148 (MQ=255)
                                                                                       cACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAg                                                                                                                                                <  1:480774/65‑1 (MQ=255)
                                                                                       cACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAg                                                                                                                                                >  2:480774/1‑65 (MQ=255)
                                                                                          tGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAattat                                                          >  1:752816/1‑148 (MQ=255)
                                                                                           gagaAAGAATTGCAACAGTAATGCCAGCTTGTGAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAattatt                                                         >  2:281066/1‑148 (MQ=255)
                                                                                            agaaagaaTTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTAttt                                                        <  1:269451/148‑1 (MQ=255)
                                                                                             gaaagaaTTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTCATTCAGCCGGTCGGAACTGGTATTTAACCAGACTAATTATTTTg                                                      >  3:107236/1‑149 (MQ=255)
                                                                                                 gaaTTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGATCTGGTATTTAACCAGACTAATTATTTTGATg                                                   >  1:419674/1‑148 (MQ=255)
                                                                                                 gaaTTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATg                                                   <  1:561290/148‑1 (MQ=255)
                                                                                                     tGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGcgcg                                               >  1:570151/1‑148 (MQ=255)
                                                                                                       cAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGaa                                             >  1:320396/1‑148 (MQ=255)
                                                                                                       cAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGaa                                             <  2:320396/148‑1 (MQ=255)
                                                                                                        aaCAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGaaa                                            <  1:317092/148‑1 (MQ=255)
                                                                                                         aCAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTAttt                                                        >  1:841518/1‑135 (MQ=255)
                                                                                                         aCAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTAttt                                                        <  2:841518/135‑1 (MQ=255)
                                                                                                            gTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTaa                                        <  1:734527/148‑1 (MQ=255)
                                                                                                                tGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGtt                                   >  4:213405/1‑149 (MQ=255)
                                                                                                                 gCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGtt                                   <  2:352579/148‑1 (MQ=255)
                                                                                                                  ccAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTa                                  >  2:17763/1‑148 (MQ=255)
                                                                                                                  ccAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTAc                                 <  4:810102/149‑1 (MQ=255)
                                                                                                                   cAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACa                                >  3:218105/1‑149 (MQ=255)
                                                                                                                   cAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACa                                <  4:180097/149‑1 (MQ=255)
                                                                                                                   cAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACa                                >  4:244811/1‑149 (MQ=255)
                                                                                                                           tAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAg                         <  1:264031/148‑1 (MQ=255)
                                                                                                                                     caaaaaaGCACCTAAAAAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGa               >  1:865080/2‑148 (MQ=255)
                                                                                                                                       aaaaaGCACCTAAAAAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAAt             >  2:726317/1‑148 (MQ=255)
                                                                                                                                       aaaaaGCACCTAAAAAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATc            >  4:272476/1‑149 (MQ=255)
                                                                                                                                        aaaaGCACCTAAAAAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCg           >  3:816484/1‑149 (MQ=255)
                                                                                                                                           aGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGAt         >  1:132625/1‑148 (MQ=255)
                                                                                                                                            gCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGAtt        >  1:864977/1‑148 (MQ=255)
                                                                                                                                              aCCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGa                                                                  >  2:865010/1‑88 (MQ=255)
                                                                                                                                              aCCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGa                                                                  <  1:865010/88‑1 (MQ=255)
                                                                                                                                              aCCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTAtttt                                                       <  3:444097/99‑1 (MQ=255)
                                                                                                                                              aCCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTAtttt                                                       >  4:444097/1‑99 (MQ=255)
                                                                                                                                               ccTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTa     <  1:17763/148‑1 (MQ=255)
                                                                                                                                                 tttttAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGa          <  2:115605/141‑1 (MQ=255)
                                                                                                                                                 tttttAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGa          >  1:115605/1‑141 (MQ=255)
                                                                                                                                                 tttttAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGa  >  4:856667/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TGAAAGCTTTACCAACAAGAATATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGA  >  NC_000913/2533269‑2533562

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