Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,138,130 C→T P254S (CCG→TCG)  flgJ → putative peptidoglycan hydrolase FlgJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,138,1300CT98.0% 139.2 / ‑5.6 52P254S (CCG→TCG) flgJputative peptidoglycan hydrolase FlgJ
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (28/22);  minor base G (0/1);  total (28/24)
Fisher's exact test for biased strand distribution p-value = 4.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.84e-01

TCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCA  >  NC_000913/1137883‑1138378
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        
ccatTTGATCCCCGCCCAGGCGGCCCGGGAATCTGGTTGGGGGAAACGGCAATTCCGCCGCGAAAACGGCGAGCCGAGCTATTACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTcgc                                                                                                                                                                                                                                                       <  2:581723/250‑1 (MQ=255)
      gacccTCGCCCAGGCGGCACTTGAATCGGGTGGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGGTACTGAAATCACCACGTCTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTGTCGCGTCTACAGCTCGGAGCTGGAAGCCTTGGCGGATGATGTTGGGCTGTTAACGCGTAACTCGCGCTAc                                                                                                                                                                                                                                                  <  2:522731/247‑1 (MQ=255)
        tCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTCTAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTAGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTAcgcc                                                                                                                                                                                                                                               >  1:558701/1‑251 (MQ=255)
        tCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCAGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGAATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTAcgcc                                                                                                                                                                                                                                               >  2:188593/1‑251 (MQ=255)
        cccTCGCCCAGGCGGCACTGGAATCTGGTTGGGGGAAACGGCAAAGCCGCCGCGAAAACGGCGAGCCGAGGTATAACCTGTTTGGTGTCAAAGCCTCTTGCAAATGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTAcgcc                                                                                                                                                                                                                                               <  2:122486/250‑1 (MQ=255)
          cTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTAcgccgc                                                                                                                                                                                                                                             >  2:65060/1‑251 (MQ=255)
                 aGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCCACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTAGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCAGCAGTgacg                                                                                                                                                                                                                                      >  2:184540/1‑251 (MQ=255)
                  ggcggaattggacgatggctGGGGGCACCGGCAACTCCACCGCGAACACGGCGAGCCGCGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCCACTGGAAAGGGCCAGTTACTGAAATCACCACGCCTGAATATGAAACCGGCGACGCGCAGAAAGTAAAAGAGAAGTTTCGCGTCTACAGCACGTATCTGGAAGCCTTGTCGGCTTACGTTGGGCTGTTCACGCGTAACTCGCGCTACGCCCCCGTgacg                                                                                                                                                                                                                                      >  2:271206/20‑250 (MQ=255)
                    ggtcACTGGAATCTGGTTGGGGTAAAGGTCAAATCCGCCGCGAAAATGGCGAGACGAGCTATACCCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGTCCAGTTACTGAAATCCCAACGACTGAATATGAAAACGTCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCTTCTACAGCACGTATCTGGAAGCCTTGCCGGATTACGTTGGGCTGTTAACGCGTACCTCGCGTTACGCCGCCTTGACGAcc                                                                                                                                                                                                                                   <  2:129538/248‑1 (MQ=255)
                    cggcACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGAcc                                                                                                                                                                                                                                   >  1:60900/1‑251 (MQ=255)
                                                               aaaCGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCGGGCAATGGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATAGGAAAACGGCGAAGTGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCTTAACTCGCGCTACGCCGCCGTTTCCCCCTCCCCTAGTTTTGAACAGGGGGCGCAGGCACTAAAGGAcgcgt                                                                                                                                                                                        <  1:395969/251‑2 (MQ=255)
                                                                  cGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCt                                                                                                                                                                                     >  2:572643/1‑251 (MQ=255)
                                                                   ggCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGTAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCtc                                                                                                                                                                                    >  2:425665/1‑250 (MQ=255)
                                                                           gagcTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCg                                                                                                                                                                            >  1:501798/1‑251 (MQ=255)
                                                                            agcTATAAGCTGTTTGGTGTCAAAACCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCATGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACCGTGGGCGCAGGCCCGACAGGACGCGGGCTATGCCCCCgc                                                                                                                                                                           >  2:486621/1‑250 (MQ=255)
                                                                             gcTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACCGGACGCGGGCTATGCCACCGAt                                                                                                                                                                          >  2:408726/1‑251 (MQ=255)
                                                                                   aCCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGCACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCAc                                                                                                                                                                    >  2:131407/1‑251 (MQ=255)
                                                                                    ccTGTTTGGTGTCAAAGACTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACt                                                                                                                                                                   >  2:208937/1‑251 (MQ=255)
                                                                                            gtgtCAAGGCCTCTGGCAACGGAAAAGGGCCAGTTACTGAAATCCCCACATCTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAATTTTCGCGTCTTCAGATCTTATCTGGAAGCTTTGTCGGATTACGTTGGGCTGTTACCGCGTACCTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGc                                                                                                                                                           <  2:437698/251‑1 (MQ=255)
                                                                                             tgtCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGCAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGc                                                                                                                                                           >  1:46424/1‑250 (MQ=255)
                                                                                                   aGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAACTCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGCGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGCTCCTCACTCTGACCGCAAACTc                                                                                                                                                     >  2:372395/1‑250 (MQ=255)
                                                                                                          ggCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACGGAATATGAAAACGGAGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCGGTTACCGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCTACCGATCCTCACTATGCCCGCAAACTCACCAACAt                                                                                                                                             <  1:165079/251‑1 (MQ=255)
                                                                                                           gCAACTGGAAAGGGCCAGATACTGAAAACACCACGACTGAATATGAAACCGCCGAAGCGAAGAAAGAAAAGGCGAAGTTTCGCGTCTACAGCGCGTATCGGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCTCAACTCGCGCTCCGCCGCCGTGACGACCGCCGCGAGGGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATg                                                                                                                                            <  1:501181/251‑1 (MQ=255)
                                                                                                              aCTGGAAAGGGCAAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAGCGCGTAACTCTCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGAtt                                                                                                                                         <  2:90859/251‑1 (MQ=255)
                                                                                                                 ggAAGGGGCAAGTTACGGAAATCACCACGACTGAATAGGAAGACGGCGAAGCGAAGAAAGGAAAGGCGAAGTTGCGGGGGTACAGCGCGTATCTGGAGGCCTTGTCGGATTACGTGGGGCGGTTAACGCGTAACTCGGGCTACGCCGCCGGGACGGCCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCGCACTATGCCCGCAAACTCACCAACAGGATGcag                                                                                                                                      <  2:46870/251‑1 (MQ=255)
                                                                                                                  gAAAGGGCCAGTTACTGCAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTcagc                                                                                                                                     >  1:458578/1‑251 (MQ=255)
                                                                                                                         ccAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGCAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCACACTCACCAACATGATTCAGCAGATGa                                                                                                                               >  1:185262/1‑250 (MQ=255)
                                                                                                                         ccAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGGCGGATTACGTTGGGCTGTTAACGCGTAACGCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGAGCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGaa                                                                                                                              <  1:69329/251‑1 (MQ=255)
                                                                                                                              tACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCg                                                                                                                          <  1:461762/250‑1 (MQ=255)
                                                                                                                                  aaaaTCACCACGACTGAAGATGAAACCGGCGAAGCGAAGAAAGTAAAAGCGAAGTGTCGCTCCTACAGCTCGTATCTGGAAGCCTTTTCGGATGACGTTGGGCTGTTACCGGGTATCCCGCGCTAGGCCGCCTTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGGCCTACAGGACGCGGGCTATGCCACCGAGCCTGACGATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATaa                                                                                                                      <  1:259482/249‑1 (MQ=255)
                                                                                                                                              aCTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTg                                                                                                           >  2:440721/1‑249 (MQ=255)
                                                                                                                                              aCTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAg                                                                                                         <  1:124040/251‑1 (MQ=255)
                                                                                                                                                   atatGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTGGGCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCaaaa                                                                                                    <  1:125482/251‑1 (MQ=255)
                                                                                                                                                        aaaaCGGCGAAGGAAAGAAAGTAAAAGCGAAGTTTCGCGTCTACGGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTCCGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTAc                                                                                               <  2:389745/251‑1 (MQ=255)
                                                                                                                                                         aaaCGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCCCTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTAc                                                                                               >  1:306851/1‑250 (MQ=255)
                                                                                                                                                             ggCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGCGCAGGCCCTACCGGACGCGGGCTATGCCAACGGTCCTAACTATGCCCGCAAACTCACCAACATGGTTTCGCCGATGAAATCGCAAAGCGAAAAGGGGAGGAAAACCTACAGTa                                                                                           >  1:241273/1‑249 (MQ=255)
                                                                                                                                                                      aagaaAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTg                                                                                  >  1:524467/1‑250 (MQ=255)
                                                                                                                                                                        caaaGTGAAAGCGAAGTTTCGCGTCTAGAGCTCGTATCTGGAAGCCTAGTCGGATTACGTTGGGCTGTTAACGCGAAACTCGCGATACGCCGCCGTGACGACCGCCGAGAGTGCGGAACAGGGGGCGCCGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCACACTCACCAACATGATTCAGCAGCTGAAATCGATCAGAGACACGGTGAGCAAAACCTACCGAATGAACATTGAt                                                                                >  2:71237/2‑250 (MQ=255)
                                                                                                                                                                          aaGTAAAAGCGAATTTTCGCGTCTTCAGCGCGTATCTGGAAGCCTTGCCGGATTACGTTGGGCTGTTAACGCGTACCTCGCGCTCCGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAATGTGAGCAAAACCTACAGTATGAACATTGATAAt                                                                             <  1:184540/251‑1 (MQ=255)
                                                                                                                                                                               aaaGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGtt                                                                        <  1:271082/251‑1 (MQ=255)
                                                                                                                                                                                                cTACAGCTCGTATTTAGAAACCTTGCCGCCTTACGTTGGGCTGTTAACGCGTAACTATCGCTACGACGCCGTATCGACCGCCGCGAGTAAGGAACAGGGGGCGCAAGCACTACATGACGCGCGCTATGCCACCCATCCTCCCTATGCCCGCAAACTCACCAACATGATTCAGCCGATGAAATCGATAACCGACAAGGTGCGCAAAACCTCCAGTATGTACATTGATAATCTGTTCTGAATAACTCAAATcc                                                       >  2:41591/1‑251 (MQ=255)
                                                                                                                                                                                                 tACAGCTCGTATCTGGCAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCg                                                                                                                          >  1:182575/1‑183 (MQ=255)
                                                                                                                                                                                                  aCAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAATCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCcgg                                                     >  1:476222/1‑251 (MQ=255)
                                                                                                                                                                                                          tATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAACCTCACCAACATGATTCAGCAGATGAAATCGATAAGCGCCAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGc                                             >  2:238398/1‑251 (MQ=255)
                                                                                                                                                                                                           atcgggAAGTCTGGTCGGATTACGTTGGGGTGTTAACGGGTAACTCGCGGTACGCGGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAAGCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCt                                            <  2:185262/247‑1 (MQ=255)
                                                                                                                                                                                                               ggAAGCCTTGTGGGATGCCGTTGGGCTGTTAACGCGTACCTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGTGTCGCTg                                           <  2:388118/248‑1 (MQ=255)
                                                                                                                                                                                                               ggAAGCCTTGTCGGCTTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCt                                            <  1:435029/247‑1 (MQ=255)
                                                                                                                                                                                                                    ccTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGataat                                   >  1:529247/1‑251 (MQ=255)
                                                                                                                                                                                                                             aTTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTaa                          <  1:133283/251‑1 (MQ=255)
                                                                                                                                                                                                                                 cGTTGGGCTGTGAACGCGAAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGCCGCTGCCGATAATACTCTGTAATTGa                      <  1:152951/251‑1 (MQ=255)
                                                                                                                                                                                                                                   ttGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGa                      >  1:44932/1‑249 (MQ=255)
                                                                                                                                                                                                                                          gTTAGCGCGGAACTCGCGCTCCGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGAAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTtat              <  2:469027/250‑1 (MQ=255)
                                                                                                                                                                                                                                                     ccccGCGCTCCGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCa  <  1:208937/250‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        
TCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCA  >  NC_000913/1137883‑1138378

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