Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,143,377 C→T S283L (TCA→TTA)  frwC → putative PTS enzyme IIC component FrwC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,143,3770CT96.5% 177.4 / ‑3.4 57S283L (TCA→TTA) frwCputative PTS enzyme IIC component FrwC
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (17/38);  total (18/39)
Fisher's exact test for biased strand distribution p-value = 5.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

CTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGG‑TTCAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGC‑AGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCGACTGATGAAAAAGAAGACGACCTGGATTTGGATTTTGAAATTAATTAATTGAGGA  >  NC_000913/4143135‑4143616
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                
cTTCATGCTGCTTTGCGTTGCTCAGGGTGGTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGG‑TTTCGTTtgtg                                                                                                                                                                                                                                          >  2:183083/1‑251 (MQ=255)
 ttCAAGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGATGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGAGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGGAACCGGCAAAGCGGCACTGGTGATGGGCTGAGGTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGG‑TTTAGTTTgtgg                                                                                                                                                                                                                                         >  2:373252/1‑251 (MQ=255)
   cATGCTGCTTTGCGTTGCTCAGGGTGTTTATCCCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGCAGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGCAACCGGCAAAGCGGCACTGGTGATGGGCTGCGTTGGCGTTACTGAAGCGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGG‑TTTAGTTTgtggtg                                                                                                                                                                                                                                       >  1:263184/1‑251 (MQ=255)
   cATGCTGATTTGCGTTGCTCAGTGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGTTATTTGTATCCCGCCGCTGGGTATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGTGCGCGAAACCGGCAAAGCGGCACTGGTGATGTGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTTCTTGTCATTCCTTCCATCATGGTCGG‑TTTAGTTTgtggt                                                                                                                                                                                                                                        <  1:433352/250‑1 (MQ=255)
      gCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGG‑TTTAGTTTGTGGTGCAg                                                                                                                                                                                                                                    <  1:109012/251‑1 (MQ=255)
        tGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGG‑TTTAGTTTGTGGTGCAGt                                                                                                                                                                                                                                   >  1:113470/1‑250 (MQ=255)
         gATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCTATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGTTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTGTCGCCGCGGCCGATCCGCTGCGTGTCATGCCTTCCATCATGGTCGG‑TTTAGTTTGTGGTGCAGTaa                                                                                                                                                                                                                                 <  1:52736/251‑1 (MQ=255)
           tttGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGACAATCGTCTCCCGCCGCTCGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTACGCAGACGAGCGCGAAACCGGCCAAGAGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGG‑TTTAGTTTGTGGTGCAGTAACt                                                                                                                                                                                                                               >  1:429966/1‑251 (MQ=255)
             tGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGG‑TTTAGTTTGTGGTGCAGTAACTg                                                                                                                                                                                                                              <  1:157038/250‑1 (MQ=255)
                  tGCTCAGGTTGTTTATACCGTGGTTGCCATCGCTGCCTTGGTTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTTGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCTGCAAAGCGGCGCTGGTGATGGGCTGTGTTGGCGTGACTGAAGGTGCGATTCCTTTCGCCGCGGCCGATCCGCTTCGTTTCATGCCGTCCTTCATGGTCGG‑TTTAGTTTGTGGTGCAGTAACTgcggcg                                                                                                                                                                                                                         <  1:198443/250‑1 (MQ=255)
                     ccggggTGTTTATCCCGGGGTGTCCATCGCTGCGGTGTGTTTTTGTATCCCTCTGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTTTCATTCCTTCCAGCATGGTCGG‑TTTAGTGTGTGGTTCATTAACTTCGGCGCTgg                                                                                                                                                                                                                     <  2:285266/248‑1 (MQ=255)
                              ttagaCCGTGGTGGCCACCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGG‑TTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTgcg                                                                                                                                                                                                             <  1:416034/246‑1 (MQ=255)
                              ttATACCGTGGCGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGCTACGCTGCGTGTCATTCCTTCCATCATGGTCGG‑TTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGACGGTgcgc                                                                                                                                                                                                            >  2:336650/1‑251 (MQ=255)
                                      tggtggCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGAGGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCCAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGATGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTTTCATTCCTTCCATCCTGGGCGG‑TTTAGTTTGTGGTGTCGTAACTGCGGCGCTGGTCGGTGCGCCGTGCTAt                                                                                                                                                                                                    >  2:174304/1‑251 (MQ=255)
                                      gggtggCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGG‑TTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTAt                                                                                                                                                                                                    <  2:94992/250‑1 (MQ=255)
                                        gtggCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGG‑TTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGc                                                                                                                                                                                                  >  1:158482/1‑251 (MQ=255)
                                        gtggCCAGCGCCGCCGTGGGTATTTGTATCCCGCCGGTGGGGGAGGGGCTGGCGACGCTGGTTGGGCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGG‑TTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGc                                                                                                                                                                                                  <  1:105987/251‑1 (MQ=255)
                                               ccgtggCCGTGGGTATTTGTAGCCCGCCGCTGGGGATGGGGCTGGCGACGCTGTGTGGTCTTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCATCTCTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCGGCGTGTCAGTCCTTCCATCATGGTCGG‑TTTAGTTTGGGGTTCAGTAACTGCGTCGCTTGTCGGTGCGCAGTGCCATGCAGGCTg                                                                                                                                                                                            <  2:124423/245‑1 (MQ=255)
                                                    gccttGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGTCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGG‑TTTAGTTTGTGGTGCAGTAACTGCGGCTCTTGTCGGTGCGCAGTGCTATGCAGGCTGGggt                                                                                                                                                                                        <  1:348819/245‑1 (MQ=255)
                                                      cGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGG‑TTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTc                                                                                                                                                                                    >  1:103534/1‑251 (MQ=255)
                                                       gTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGG‑TTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCt                                                                                                                                                                                   <  2:103534/251‑1 (MQ=255)
                                                                    ccccGCCGCTGGGGTTGGGGCGGGCGACGCTGATTTGTCGTAAAAAGTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGG‑TTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTTTGCTGCCg                                                                                                                                                                      <  2:347773/250‑1 (MQ=255)
                                                                                                         tCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGG‑TTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGC‑AGTgg                                                                                                                                >  1:295036/1‑251 (MQ=255)
                                                                                                                  tttctgccccgaagtGCGAGAAACCGACAAACCAGCACTGGTGATCGGCTGAGAAGGCGGTAATGGCGGGGAGAATCCTTTCGCCGCTGCCGATCCGCTGCGTGTCAGTCCTTCCATCATGGTCGGATTC‑GTTTGTGGAGCAGTAACTGCGGCGCTGGTCGGTGCGCGGTGCTATGCAGGCTGGGGGAGTCGGATTGTGCTGCCGGAGGTGGAAGGCAAGCTGGGTTATATCGCAGCCAG‑GGCTGTCggcg                                                                                                                       >  2:394569/11‑249 (MQ=255)
                                                                                                                       ccGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGG‑TTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGC‑AGTGGCTGTCGGAGCAg                                                                                                                    <  1:67571/249‑1 (MQ=255)
                                                                                                                          cAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGG‑TTTAGTTTGTGGTGCAGTAACTTCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGC‑AGTGGCTGTCGGAGCAgtggtg                                                                                                               <  1:198943/251‑1 (MQ=255)
                                                                                                                             aagacggcgAAACCGGCAAAGCGGCACTGGAGAGGGGCTGTGTTGGCTTTTCTGAAGGGGCGAATCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTGCCATCATGGTCGG‑TTTAGTTTGAGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGGTATATCGCAGC‑AGTGGCTGTCGGAGCAGTGGTGACg                                                                                                            <  2:130462/245‑1 (MQ=255)
                                                                                                                              agagCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGG‑TTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGC‑AGTGGCTGTCGGAGCAGTGGTGACgg                                                                                                           <  1:100686/251‑1 (MQ=255)
                                                                                                                               gagCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGGTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGGTGCCGGTCCGCTGCGGGTCATTCCTGCCATCATGGTCGG‑TTTAGGTTGGGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGGGGTTGAAGGCAAGCTGGGTTATATCGCAGC‑AGTGGCTGTCGGAGCAGTGGTGACGGc                                                                                                          <  1:155842/251‑1 (MQ=255)
                                                                                                                                  cgcgAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGG‑TTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGC‑AGTGGCTGTCGGAGCAGTGGTGACGGCTGt                                                                                                       >  1:32438/1‑251 (MQ=255)
                                                                                                                                  cgcgAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGG‑TTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGC‑AGTGGCTGTCGGAGCAGTGGTGACGGCTGt                                                                                                       <  1:362902/251‑1 (MQ=255)
                                                                                                                                       aaCCGGCAAAGCGGCTCTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCGTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGG‑TTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGC‑AGTGGCTGTCGGAGCAGTGTTGACGGCTGTTtgtg                                                                                                  <  2:436333/251‑1 (MQ=255)
                                                                                                                                        cccGGCAAAGCGTCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATGCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGG‑TTCATTTTGTGGTGCATTAACTGCGGCGCTTGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGTTCTGATTGTGCTGCCGGTGGTTTAAGGCAACCTGGGTTATATCGCAGC‑AGTGGCTGTCGGATCAGTTGTGACGGCTGTTtgtgt                                                                                                 <  2:134855/250‑1 (MQ=255)
                                                                                                                                        aCCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGG‑TTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGC‑AGTGGCTGTCGGAGCAGTGGTGACGGCTGTTtgtgt                                                                                                 <  1:113076/251‑1 (MQ=255)
                                                                                                                                         ccGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGG‑TTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGC‑AGTGGCTGTCGGAGCAGTGGTGACGGCTGTTtgtgtg                                                                                                >  1:265228/1‑251 (MQ=255)
                                                                                                                                         ccGGCAAAGCGGAACTGGTGATGGGATGTGCTGGATTTACTGAAGGGGCGATTCCGTTCGCGGCTGCCGATCCGCTGCGTGTCATTTCTTCCATCATGGTCGG‑TGTAGTTTGTGGTGCAGTAACTGCGGCGCTAGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGGCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGC‑AGTGGCTGTCGGAGCAGTGGTGACGGATGTTtgtgtg                                                                                                <  1:313221/251‑1 (MQ=255)
                                                                                                                                          cGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGG‑TTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGC‑AGTGGCTGTCGGATCAGTGGTGACGGCTGTGTGTGTGa                                                                                               <  1:108637/251‑1 (MQ=255)
                                                                                                                                             caaaggtgCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATCCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGG‑TGTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGC‑AGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACg                                                                                            <  2:191221/244‑1 (MQ=255)
                                                                                                                                               aaGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGG‑TTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGC‑AGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTg                                                                                          <  1:76471/251‑1 (MQ=255)
                                                                                                                                                  cGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGG‑TTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGC‑AGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTg                                                                                       <  1:283804/251‑1 (MQ=255)
                                                                                                                                                      aCTGGTGCTGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGG‑TTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGC‑AGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGaaa                                                                                    >  1:119564/1‑250 (MQ=255)
                                                                                                                                                                        tGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGG‑TTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGC‑AGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAAt                                                                  <  1:386441/250‑1 (MQ=255)
                                                                                                                                                                               aCTGAAGGGGCGAGTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGG‑TTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGC‑AGTGGCTGTCGGAGCAGTGGTTACGGCGGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCTTAAAAATGGGTCTTc                                                          <  2:366205/251‑1 (MQ=255)
                                                                                                                                                                                  gAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGG‑TTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGC‑AGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCg                                                         <  1:421984/249‑1 (MQ=255)
                                                                                                                                                                                   aaGGGGCGATTCCTTTCGACGCTGCCGATCCGCTGCGTGTCGTTCCTTCAATCATGGTCGG‑TTTAGTTTGTGGTGCCGTAACTGCGGCGATGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGGGGGTGAAGGCAAGCTGGGTTATATCGCCGC‑CGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAAAGTGCTGAAACGTCTGGGTGGTCAAAATGGGTCTTCGACt                                                      >  2:200317/1‑251 (MQ=255)
                                                                                                                                                                                          gATTCCTTTCTCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGG‑TTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGC‑AGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCGACGGATGaaa                                               <  1:322978/251‑1 (MQ=255)
                                                                                                                                                                                                  tCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGG‑TTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGC‑AGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCGACTGATGAAAAagaaga                                        <  2:361230/250‑1 (MQ=255)
                                                                                                                                                                                                    gccgcTGCCGATCCGCTGCGTGTCATTCCTTCCAGCATGGTCGG‑TTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGGGGTTGAAGGCAAGCTTGGTTATATCGCAGC‑AGTGGCTGTCGGAGCAGTGGTGACGGCTTTTTGTGTGTACGGGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTGCGACTGATGAAAAAGAAgacga                                     <  1:394569/251‑1 (MQ=255)
                                                                                                                                                                                                      cgcTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGG‑TTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGC‑AGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCGACTGATGAAAAAGAAGACGAcc                                   <  1:336650/251‑1 (MQ=255)
                                                                                                                                                                                                          gCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGG‑TTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGC‑AGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCGACTGATGAAAAAGAAGACGACCtg                                 <  1:68650/249‑1 (MQ=255)
                                                                                                                                                                                                            cGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGG‑TTTAGTTTGTGGTGCAGTAACGGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGGTGGGGTGGTCTGAGTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGG‑AGTGGCTTTCGGAGCAGTGGTGACGGCTTTTTGTGTGAACGTGCTGAAAAGGCTGGCGCGTAAAAATGGGTCTTCTGCTGATGAAAAAGAAGACGTCCtggatt                             <  2:263184/251‑1 (MQ=255)
                                                                                                                                                                                                             gATCCGCTGCGTGTCATTCCTTCCATCATGGTCGG‑TTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGC‑AGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCGACTGATGAAAAAGAAGACGACCtggattt                            >  1:370295/1‑251 (MQ=255)
                                                                                                                                                                                                              aTCCGCTGCGTGTCATTCCTTCCATCATGGTCGG‑TTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGC‑AGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCGACTGATGAAAAAGAAGACGACCtggatttg                           <  2:398493/251‑1 (MQ=255)
                                                                                                                                                                                                                    tGCGTGTCATTCCTTCCATCATGGTCGG‑TTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGC‑AGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCGACTGATGAAAAAGAAGACGACCTGGATTTGGAtttt                     >  1:351701/1‑251 (MQ=255)
                                                                                                                                                                                                                    tGCGTGTCATTCATTCCATCATGGTAGG‑TTTAGTTTGTGGTGCAGTAACTGCGACGCTGGTCGGTGCGAAGTGCTATGCAGGCTGGGGCGGTCTGATTGTGCTGCCGGTGGTTGAAGGCACGCTGGGTTATATCGCAGC‑AGTGGCTGTCGGAGGCGTGGTGGCGGACGTTTGTGTGAACGAACTGACAAGACGGGCGCGTAAAAATGGGTCTTCGACTGATGAAAAAGAAGACGACCaggatttgactttt                     >  2:411207/1‑245 (MQ=255)
                                                                                                                                                                                                                         gtCATTCCTTCCATCATGGTCGG‑TTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGC‑AGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCGACTGATGAAAAAGAAGACGACCTGGATTTGGATTTTGaaa                 <  2:143485/250‑1 (MQ=255)
                                                                                                                                                                                                                                       catGGTCGG‑TTTAGGTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGC‑AGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCGACTGATGAAAAAGAAGACGACCTGGATTTGGATTTTGAAATTAATTAATTGAGGa  <  2:104488/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                
CTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGG‑TTCAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGC‑AGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCGACTGATGAAAAAGAAGACGACCTGGATTTGGATTTTGAAATTAATTAATTGAGGA  >  NC_000913/4143135‑4143616

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