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

CTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCGACTGATGAAAAAGAAGACGACCTGGATTTGGATTTTGAAATTAATTAATTGAGGA  >  NC_000913/4143135‑4143616
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                
cTTCATGCTGCTTTGCGTTGCTCAGGGTGGTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTCGTTtgtg                                                                                                                                                                                                                                          >  2:183083/1‑251 (MQ=255)
 ttCAAGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGATGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGAGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGGAACCGGCAAAGCGGCACTGGTGATGGGCTGAGGTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTgtgg                                                                                                                                                                                                                                         >  2:373252/1‑251 (MQ=255)
   cATGCTGCTTTGCGTTGCTCAGGGTGTTTATCCCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGCAGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGCAACCGGCAAAGCGGCACTGGTGATGGGCTGCGTTGGCGTTACTGAAGCGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTgtggtg                                                                                                                                                                                                                                       >  1:263184/1‑251 (MQ=255)
   cATGCTGATTTGCGTTGCTCAGTGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGTTATTTGTATCCCGCCGCTGGGTATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGTGCGCGAAACCGGCAAAGCGGCACTGGTGATGTGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTTCTTGTCATTCCTTCCATCATGGTCGGTTTAGTTTgtggt                                                                                                                                                                                                                                        <  1:433352/250‑1 (MQ=255)
      gCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAg                                                                                                                                                                                                                                    <  1:109012/251‑1 (MQ=255)
        tGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGt                                                                                                                                                                                                                                   >  1:113470/1‑250 (MQ=255)
         gATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCTATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGTTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTGTCGCCGCGGCCGATCCGCTGCGTGTCATGCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTaa                                                                                                                                                                                                                                 <  1:52736/251‑1 (MQ=255)
           tttGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGACAATCGTCTCCCGCCGCTCGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTACGCAGACGAGCGCGAAACCGGCCAAGAGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACt                                                                                                                                                                                                                               >  1:429966/1‑251 (MQ=255)
             tGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTg                                                                                                                                                                                                                              <  1:157038/250‑1 (MQ=255)
                  tGCTCAGGTTGTTTATACCGTGGTTGCCATCGCTGCCTTGGTTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTTGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCTGCAAAGCGGCGCTGGTGATGGGCTGTGTTGGCGTGACTGAAGGTGCGATTCCTTTCGCCGCGGCCGATCCGCTTCGTTTCATGCCGTCCTTCATGGTCGGTTTAGTTTGTGGTGCAGTAACTgcggcg                                                                                                                                                                                                                         <  1:198443/250‑1 (MQ=255)
                     ccggggTGTTTATCCCGGGGTGTCCATCGCTGCGGTGTGTTTTTGTATCCCTCTGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTTTCATTCCTTCCAGCATGGTCGGTTTAGTGTGTGGTTCATTAACTTCGGCGCTgg                                                                                                                                                                                                                     <  2:285266/248‑1 (MQ=255)
                              ttagaCCGTGGTGGCCACCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTgcg                                                                                                                                                                                                             <  1:416034/246‑1 (MQ=255)
                              ttATACCGTGGCGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGCTACGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGACGGTgcgc                                                                                                                                                                                                            >  2:336650/1‑251 (MQ=255)
                                      tggtggCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGAGGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCCAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGATGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTTTCATTCCTTCCATCCTGGGCGGTTTAGTTTGTGGTGTCGTAACTGCGGCGCTGGTCGGTGCGCCGTGCTAt                                                                                                                                                                                                    >  2:174304/1‑251 (MQ=255)
                                      gggtggCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTAt                                                                                                                                                                                                    <  2:94992/250‑1 (MQ=255)
                                        gtggCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGc                                                                                                                                                                                                  >  1:158482/1‑251 (MQ=255)
                                        gtggCCAGCGCCGCCGTGGGTATTTGTATCCCGCCGGTGGGGGAGGGGCTGGCGACGCTGGTTGGGCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGc                                                                                                                                                                                                  <  1:105987/251‑1 (MQ=255)
                                               ccgtggCCGTGGGTATTTGTAGCCCGCCGCTGGGGATGGGGCTGGCGACGCTGTGTGGTCTTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCATCTCTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCGGCGTGTCAGTCCTTCCATCATGGTCGGTTTAGTTTGGGGTTCAGTAACTGCGTCGCTTGTCGGTGCGCAGTGCCATGCAGGCTg                                                                                                                                                                                            <  2:124423/245‑1 (MQ=255)
                                                    gccttGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGTCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCTCTTGTCGGTGCGCAGTGCTATGCAGGCTGGggt                                                                                                                                                                                        <  1:348819/245‑1 (MQ=255)
                                                      cGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTc                                                                                                                                                                                    >  1:103534/1‑251 (MQ=255)
                                                       gTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCt                                                                                                                                                                                   <  2:103534/251‑1 (MQ=255)
                                                                    ccccGCCGCTGGGGTTGGGGCGGGCGACGCTGATTTGTCGTAAAAAGTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTTTGCTGCCg                                                                                                                                                                      <  2:347773/250‑1 (MQ=255)
                                                                                                         tCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTgg                                                                                                                                >  1:295036/1‑251 (MQ=255)
                                                                                                                  tttctgccccgaagtGCGAGAAACCGACAAACCAGCACTGGTGATCGGCTGAGAAGGCGGTAATGGCGGGGAGAATCCTTTCGCCGCTGCCGATCCGCTGCGTGTCAGTCCTTCCATCATGGTCGGATTCGTTTGTGGAGCAGTAACTGCGGCGCTGGTCGGTGCGCGGTGCTATGCAGGCTGGGGGAGTCGGATTGTGCTGCCGGAGGTGGAAGGCAAGCTGGGTTATATCGCAGCCAGGGCTGTCggcg                                                                                                                       >  2:394569/11‑249 (MQ=255)
                                                                                                                       ccGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAg                                                                                                                    <  1:67571/249‑1 (MQ=255)
                                                                                                                          cAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTTCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAgtggtg                                                                                                               <  1:198943/251‑1 (MQ=255)
                                                                                                                             aagacggcgAAACCGGCAAAGCGGCACTGGAGAGGGGCTGTGTTGGCTTTTCTGAAGGGGCGAATCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTGCCATCATGGTCGGTTTAGTTTGAGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGGTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACg                                                                                                            <  2:130462/245‑1 (MQ=255)
                                                                                                                              agagCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACgg                                                                                                           <  1:100686/251‑1 (MQ=255)
                                                                                                                               gagCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGGTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGGTGCCGGTCCGCTGCGGGTCATTCCTGCCATCATGGTCGGTTTAGGTTGGGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGGGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGc                                                                                                          <  1:155842/251‑1 (MQ=255)
                                                                                                                                  cgcgAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGt                                                                                                       >  1:32438/1‑251 (MQ=255)
                                                                                                                                  cgcgAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGt                                                                                                       <  1:362902/251‑1 (MQ=255)
                                                                                                                                       aaCCGGCAAAGCGGCTCTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCGTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGTTGACGGCTGTTtgtg                                                                                                  <  2:436333/251‑1 (MQ=255)
                                                                                                                                        cccGGCAAAGCGTCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATGCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCATTTTGTGGTGCATTAACTGCGGCGCTTGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGTTCTGATTGTGCTGCCGGTGGTTTAAGGCAACCTGGGTTATATCGCAGCAGTGGCTGTCGGATCAGTTGTGACGGCTGTTtgtgt                                                                                                 <  2:134855/250‑1 (MQ=255)
                                                                                                                                        aCCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTtgtgt                                                                                                 <  1:113076/251‑1 (MQ=255)
                                                                                                                                         ccGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTtgtgtg                                                                                                >  1:265228/1‑251 (MQ=255)
                                                                                                                                         ccGGCAAAGCGGAACTGGTGATGGGATGTGCTGGATTTACTGAAGGGGCGATTCCGTTCGCGGCTGCCGATCCGCTGCGTGTCATTTCTTCCATCATGGTCGGTGTAGTTTGTGGTGCAGTAACTGCGGCGCTAGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGGCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGATGTTtgtgtg                                                                                                <  1:313221/251‑1 (MQ=255)
                                                                                                                                          cGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGATCAGTGGTGACGGCTGTGTGTGTGa                                                                                               <  1:108637/251‑1 (MQ=255)
                                                                                                                                             caaaggtgCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATCCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTGTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACg                                                                                            <  2:191221/244‑1 (MQ=255)
                                                                                                                                               aaGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTg                                                                                          <  1:76471/251‑1 (MQ=255)
                                                                                                                                                  cGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTg                                                                                       <  1:283804/251‑1 (MQ=255)
                                                                                                                                                      aCTGGTGCTGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGaaa                                                                                    >  1:119564/1‑250 (MQ=255)
                                                                                                                                                                        tGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAAt                                                                  <  1:386441/250‑1 (MQ=255)
                                                                                                                                                                               aCTGAAGGGGCGAGTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTTACGGCGGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCTTAAAAATGGGTCTTc                                                          <  2:366205/251‑1 (MQ=255)
                                                                                                                                                                                  gAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCg                                                         <  1:421984/249‑1 (MQ=255)
                                                                                                                                                                                   aaGGGGCGATTCCTTTCGACGCTGCCGATCCGCTGCGTGTCGTTCCTTCAATCATGGTCGGTTTAGTTTGTGGTGCCGTAACTGCGGCGATGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGGGGGTGAAGGCAAGCTGGGTTATATCGCCGCCGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAAAGTGCTGAAACGTCTGGGTGGTCAAAATGGGTCTTCGACt                                                      >  2:200317/1‑251 (MQ=255)
                                                                                                                                                                                          gATTCCTTTCTCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCGACGGATGaaa                                               <  1:322978/251‑1 (MQ=255)
                                                                                                                                                                                                  tCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCGACTGATGAAAAagaaga                                        <  2:361230/250‑1 (MQ=255)
                                                                                                                                                                                                    gccgcTGCCGATCCGCTGCGTGTCATTCCTTCCAGCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGGGGTTGAAGGCAAGCTTGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTTTTTGTGTGTACGGGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTGCGACTGATGAAAAAGAAgacga                                     <  1:394569/251‑1 (MQ=255)
                                                                                                                                                                                                      cgcTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCGACTGATGAAAAAGAAGACGAcc                                   <  1:336650/251‑1 (MQ=255)
                                                                                                                                                                                                          gCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCGACTGATGAAAAAGAAGACGACCtg                                 <  1:68650/249‑1 (MQ=255)
                                                                                                                                                                                                            cGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACGGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGGTGGGGTGGTCTGAGTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGGAGTGGCTTTCGGAGCAGTGGTGACGGCTTTTTGTGTGAACGTGCTGAAAAGGCTGGCGCGTAAAAATGGGTCTTCTGCTGATGAAAAAGAAGACGTCCtggatt                             <  2:263184/251‑1 (MQ=255)
                                                                                                                                                                                                             gATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCGACTGATGAAAAAGAAGACGACCtggattt                            >  1:370295/1‑251 (MQ=255)
                                                                                                                                                                                                              aTCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCGACTGATGAAAAAGAAGACGACCtggatttg                           <  2:398493/251‑1 (MQ=255)
                                                                                                                                                                                                                    tGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCGACTGATGAAAAAGAAGACGACCTGGATTTGGAtttt                     >  1:351701/1‑251 (MQ=255)
                                                                                                                                                                                                                    tGCGTGTCATTCATTCCATCATGGTAGGTTTAGTTTGTGGTGCAGTAACTGCGACGCTGGTCGGTGCGAAGTGCTATGCAGGCTGGGGCGGTCTGATTGTGCTGCCGGTGGTTGAAGGCACGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGGCGTGGTGGCGGACGTTTGTGTGAACGAACTGACAAGACGGGCGCGTAAAAATGGGTCTTCGACTGATGAAAAAGAAGACGACCaggatttgactttt                     >  2:411207/1‑245 (MQ=255)
                                                                                                                                                                                                                         gtCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCGACTGATGAAAAAGAAGACGACCTGGATTTGGATTTTGaaa                 <  2:143485/250‑1 (MQ=255)
                                                                                                                                                                                                                                       catGGTCGGTTTAGGTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCGACTGATGAAAAAGAAGACGACCTGGATTTGGATTTTGAAATTAATTAATTGAGGa  <  2:104488/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                
CTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCGACTGATGAAAAAGAAGACGACCTGGATTTGGATTTTGAAATTAATTAATTGAGGA  >  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: