Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,142,645 G→T 19.8% L39F (TTG→TTT frwC → putative enzyme IIC component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,142,6450GT19.8% 158.5 / 4.4 74L39F (TTG→TTTfrwCputative enzyme IIC component of PTS
Reads supporting (aligned to +/- strand):  ref base G (10/48);  new base T (13/3);  total (23/51)
Fisher's exact test for biased strand distribution p-value = 3.70e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

ATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTAACAGCTTTGGTGCGGGCTTCTTTGGTGCACTGATCGCCGGGATTATCGGCGGCATCGTGGTGCAT  >  NC_000913/4142412‑4142882
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                 
atatTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCCTCCTTGTGTATgg                                                                                                                                                                                                                                              >  1:804074/1‑239 (MQ=255)
  atTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGCCATGTTGTGTGGcc                                                                                                                                                                                                                                            >  1:174339/1‑238 (MQ=255)
  atTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGGTTCCGTCCTGTTGTATGGc                                                                                                                                                                                                                                             >  1:104366/1‑238 (MQ=255)
   ttACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCaa                                                                                                                                                                                                                                           <  1:164088/239‑1 (MQ=255)
        gttgtaaagggacaggCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTCCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGgcg                                                                                                                                                                                                                                      <  1:880028/225‑1 (MQ=255)
        aaggTCCATAAAATGGGCTTTTTATCCACTTTTTTCACCGCCTGCTATTGCTTATTGTTTCCCCCATCAAAAATATGGACAGGGGCTTTCGAAAGCCCGGGAACAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGgcg                                                                                                                                                                                                                                      <  2:843330/236‑1 (MQ=255)
        aaatttctgTAAATGGCTTTTTTTACCACTGTTTTCCCCGCCTGCTATTGCCCATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGAAAGCCCGGGATCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGgcg                                                                                                                                                                                                                                      <  2:859027/231‑1 (MQ=255)
        aaatgccAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGgcg                                                                                                                                                                                                                                      <  1:4308/234‑1 (MQ=255)
        aaaGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGgcg                                                                                                                                                                                                                                      <  1:110290/239‑1 (MQ=255)
        aaaGGCCAGTAAATGGCCTTTTTATCCAATGGTTGCTCCGCCTGCGATTGCCTATTGTTGCCTCCATCAAAAATATGGACATGGGATTTCGCAAGCCCGGGAGCAGCTTATGAATGAGGTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGgcg                                                                                                                                                                                                                                      <  1:749662/239‑1 (MQ=255)
         aattCCCGTAAATTGCCTTTTTATCCACTGTTTGCTCCGCCTGGGATTGCCTATTGTTTCCTCCAACAAAAATATGGACATGGGCTTTCGAAAGCCCGGGAGGAGCGTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGgcgc                                                                                                                                                                                                                                     <  2:896506/236‑1 (MQ=255)
          aGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGAATTTTGCTGGCGGTTTCCCTCCTGTTTTCTGGCAAAAGCGCa                                                                                                                                                                                                                                    >  1:535841/1‑239 (MQ=255)
            tCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGGTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTTCACCCCGTTGTATGGTACAGGGGCAGt                                                                                                                                                                                                                                 >  1:481194/1‑239 (MQ=255)
                  aattGGCCTTTTTACCCACTTTTTGCCCCGCCTGCGATTGCCTATTGTTTCCTCCAACAAAAATATGGACATGGGCTTTCGAAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCgg                                                                                                                                                                                                                            <  2:535841/235‑1 (MQ=255)
                      ggCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGAATTTTGCTGGCGGATTTCCTCTTGTTTTTTGGACAAGgcgcgatgacggggtcc                                                                                                                                                                                                                        >  1:376184/1‑226 (MQ=255)
                           ttttAACCAGTGTGTGCCCCGCCGGCGGTTGCCTAGTGGTTCCTCCATCAAAAATAGGGACAGGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATGTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGc                                                                                                                                                                                                                 <  1:809868/239‑1 (MQ=255)
                                    ctttttGCCCCGCCTGCGTTTTCCTTTTTGTTCCGCCAACAAAAATATGGCCATGGGCTTGCGCAAGCCCGTGGGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATTCCCGGCAGCAGTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCaaa                                                                                                                                                                                                       <  2:376184/236‑1 (MQ=255)
                                     tGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGGATTTTGCTGGCGGGTTCTGCCTAGTTTTGTGGAAAAGGAGCGGTGGCggggggcggaggcggagacaat                                                                                                                                                                                                        >  2:646868/1‑219 (MQ=255)
                                      gTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGGTTTCTTCCTGTTTTATGGCAAAGGGGCAAGGGCGGATGTCCGAGCCCGTCCAAAtt                                                                                                                                                                                                     >  1:427275/1‑238 (MQ=255)
                                         tGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAACATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCATTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTTCGGCTCTCTGTTTGTGAAAGGGGGAATGCCGGAGGGCGTAGCCggcccacaccaga                                                                                                                                                                                                   >  2:859874/1‑227 (MQ=255)
                                              cGCCTGCTATTGCCTTTTTTTTCCTCCATCAAAAATAGGGACATGGGGTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTGCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGaaaaaa                                                                                                                                                                                             <  2:619316/239‑1 (MQ=255)
                                               gCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGTGAGCAGCTTATGCATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGTCGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCCTGCGGACTTTGTCATCCTTTATGGGAACCGCGCGGTGGCGCGTGCCGtttactacctaccactgcacaaa                                                                                                                                                                                              >  2:30362/1‑216 (MQ=255)
                                               gCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTTCTGCATGTTGTTTGGCAAAACGCCAATGGCGGAGGCCGTAGcgcacacacaactgaacaa                                                                                                                                                                                               >  1:101015/1‑219 (MQ=255)
                                               gCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCCTGTTGTATGGCCAAGGCGCAGGGCCGGAGGCCCTAGCCGGTCCCAACCTGcaaaa                                                                                                                                                                                              >  1:643443/1‑237 (MQ=255)
                                               gCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGACGGGATTTTGCTGGCGGGTTTCTTCTCGTTTTATGGGAAAGGCGCAGGGCCGGGAGCGCGTGCCGACGCAAAACTGaaaaaaa                                                                                                                                                                                            >  2:693020/1‑238 (MQ=255)
                                                      aTTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGGTTCCTGCATGTTGTTTTGGAAAGGCGCGGGGCCGGGGGCCGGGACCGCACCCAATCTGaaaaaaaatgttgt                                                                                                                                                                                     >  2:649369/1‑231 (MQ=255)
                                                          ccTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCTTCACTTTTTTTGGGCAAGGGGCAGGTGCCGATGCCGGAGGCCATCCAACTTTGGAAAAAATGTTTGattta                                                                                                                                                                                 >  1:236577/1‑236 (MQ=255)
                                                               tGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCACGTTTTATGGCAAAGGGGGCGTGCCGGATGCCGTACGCGATCCCAAATTTAAAAAAATGTTTTatttttcggt                                                                                                                                                                            >  1:705317/1‑231 (MQ=255)
                                                                  ttCCTCCTACAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATTCCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGc                                                                                                                                                                         <  1:190970/239‑1 (MQ=255)
                                                                  tgtctagaTCAAAAATATGGACAGGGGCTTTCGAAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCGGAAAAATTCCCGTCAGCATTTAATGACGGTCGTCTCACACATGATTCCCTTCGTGTTATCGGGTGGTATTTTGCTGGCGGTTTCCGTCATGTTTTATGGCAAAGGCGCATTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGc                                                                                                                                                                         <  1:369994/232‑1 (MQ=255)
                                                                      tccATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTTCGTCATTTTGTATTGCCAAGGCGCAGGTGCGGGTTGCGGATCCGAGTCCAACTTGGAAAAACATTTTTAGATCGGCGtttggggg                                                                                                                                                                     >  1:567334/1‑233 (MQ=255)
                                                                                atatGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGa                                                                                                                                                          <  1:290018/239‑1 (MQ=255)
                                                                                 tttGGACATGGGGTTTCCAAAGCCCGGGGGCAGTTTATTAATGGGTTGGTTCAGATCCTGAAAAATCCCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGAt                                                                                                                                                         <  2:303030/237‑1 (MQ=255)
                                                                                  atGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTTTATGGCCAAGGGGCGGTGCCCGAAGCCGTAGCCGATCCCAAATTGAAAAAACTGTTTTATATCGGCGGTTCGGGGtttgcgcgcttg                                                                                                                                                        >  1:845608/1‑229 (MQ=255)
                                                                                  atGGACAGGGGCTTTCGCAAGCCCGGGAGCAGGTTATGAATGAGTTGGGGTAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATg                                                                                                                                                        <  1:892857/239‑1 (MQ=255)
                                                                                      ggaTGGGTTTTCGCAAGCCCGGTAGCAGTTTATAAAAGAGTTGGTGAAGATCCTGAAAAATACCCGTCAGCATTTAATGCGGGTCGTTTCACACATGATTCCCGTCGTGGGATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTGGTATGGCGAAGGCGCGGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTg                                                                                                                                                     <  2:831115/236‑1 (MQ=255)
                                                                                           ggCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTTCGGTATGTTTGTTTGCAAAAGGGCGGGGTCCGGTGCCGCATCCCGTCCAAAACTTGAAAAACCGTTTGTTATGCGGGCTGGGGGCGTTGACGcggggggggccct                                                                                                                                                 >  2:734300/1‑226 (MQ=255)
                                                                                           ggCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATTGCAAAGGCGCCGGTCCGGATGCCGCAGCCGGTCCCAATCTGAAAAAAATGCTTGGTATCGGCGTTGCCGGGTTTACGCTGATGGTGCCttt                                                                                                                                                >  1:29350/1‑238 (MQ=255)
                                                                                               ttccccacccccGGGGCAGTTTTTGAATGGGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCg                                                                                                                                           <  2:474374/231‑1 (MQ=255)
                                                                                                 ggaaggCCCGTGGGCAGTTTTAGAATGGGTTGGGTAAAGTCCTGAAAAATACCCGCCAGAATGTAATGGGGGGCGTTTCACACTTGATTCCCTTCGTGGTATGGGGCGGTATTTTGCTGGCGGTTTCCGTCAAGTTTTATGGCAAAGGCGCAGTGCCGGTTGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCa                                                                                                                                         <  2:62450/234‑1 (MQ=255)
                                                                                                        ccgggggatgattATGATTGGGTTGGTGCAGATACGGAAAAACACCCGTCAGCATTTAATGACGGGCGTTTCCAAAATGGTTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCACGCAGCGTAc                                                                                                                                   <  2:174339/227‑1 (MQ=255)
                                                                                                        ccGGGGGCCGCTTATGAATGGGGTGGTGCAGATCCTTGAAAATGCCCTGCAGCATTTAATGACGGGCTGTTCATACCGTATTCCCTTCGTGGTATCGGGCAGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACa                                                                                                                                  <  2:705317/239‑1 (MQ=255)
                                                                                                        ccGGGGGCAGCTTATGAAGGAGTTGGTGCAGATCCGTAAAAATCCCCGTAAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCAGGTTGTATGGCAAAGGCGCAGTCCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACa                                                                                                                                  <  2:603899/239‑1 (MQ=255)
                                                                                                        ccGGGCACAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATGCCCGTCAGCATTTAATGACGGGCGTTTCACGCATGATGCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACa                                                                                                                                  <  1:337994/239‑1 (MQ=255)
                                                                                                        ccGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTGTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACa                                                                                                                                  <  1:548565/239‑1 (MQ=255)
                                                                                                                    tATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCa                                                                                                                      <  1:649369/239‑1 (MQ=255)
                                                                                                                           gttttGGGGCAGATCCTGAAAATTCCCCGTTAGGTTTTAAGTGCGGGCGGTTCACACATGATTCCCTTCGGGGGATGGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCaga                                                                                                               <  2:359892/236‑1 (MQ=255)
                                                                                                                           gaGTTGGTGCAGATCCTGAAAAATCCCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCaga                                                                                                               <  1:485331/239‑1 (MQ=255)
                                                                                                                                      gATCCTGAAAAATCCCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTgcg                                                                                                     <  1:156891/238‑1 (MQ=255)
                                                                                                                                                   aCCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTTCGTCATGTTTTATGGCAAAGGCGCAGTGCCGGATGCCGTAACCGGTCCACATCTGTAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCCTATGGTGCCTTTCCTCTCAGCGGACATCGGGTTTTCCCATTCAGAGCGTTCcgcggcgtcgccgggc                                                                                         >  2:823264/1‑225 (MQ=255)
                                                                                                                                                      cGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCGTCGTGGTATCGGGCGGTACTTTGCTGGCGGGTTCCGTCATGTTGTATGGCAAGGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCATAGCGTTCTGCGCTGGCTCCGTGCGCTAt                                                                                    <  1:65965/239‑1 (MQ=255)
                                                                                                                                                           gCATTTAATGACGGGCGTTTCCCACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTg                                                                               <  1:479580/239‑1 (MQ=255)
                                                                                                                                                             aTTTAATGACGGGCTTTTCCCACATGATTCCCTTCGGGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGc                                                                              <  1:823264/238‑1 (MQ=255)
                                                                                                                                                              ttttATGACTGTCGTTTCCCAATGTATTCCCTTCGTGGTTTCGGGCGGTTTTTTGCTGGCGGTTTCCTTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCCTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCt                                                                            <  2:374787/239‑1 (MQ=255)
                                                                                                                                                              tttAATGACGGGCGTTTCCCACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCt                                                                            <  1:504755/239‑1 (MQ=255)
                                                                                                                                                              tttAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCt                                                                            <  1:646868/239‑1 (MQ=255)
                                                                                                                                                                 aaTGGCGGGCGTTTCCCTCTTGGTTCCCGTCGTGGTATCGGGGGGTATTTTGCTGGCGGTTTCCGTCATGTTGGATGGAAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTggg                                                                         <  2:29350/239‑1 (MQ=255)
                                                                                                                                                                 aaTGACGGGCGTTGCACACATGATTCCCGTCGTGGCATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGCGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTggg                                                                         <  1:859874/239‑1 (MQ=255)
                                                                                                                                                                      cGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGGTTTCCGCATGTTTTATTGCAAAAGCGCAGTGCCGGGTGCCGTTGCCTAACCCAATCTGTAAAAACTGTTTGAGTTCTGGGCTGCGGGCGTGACGCCTGTGGGGCCTTTTCTCCCGCCGCTCATCTGTTATTTCCTTTCCGAGCGGTCTGCGCCGGCCCCGTggccccccgtggccgggggtggg                                                                    >  2:14471/1‑217 (MQ=255)
                                                                                                                                                                      cGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGGCCTTTTTTTTGGCCAAGGCGCAGTGGCGGATTCCCGAGCCGGTCCAAATATGAAAAAACTGTTTGATATCGGCGTTGCGGGCGTTACCCGGATGGGGCCTTTCCCCCCCGGGGACACTGGTTTTTCCATTGCACAGCGTTTTTCCCTCGCTCCGCGCCCTTTCCGGGCCCgggggggg                                                                    >  2:842325/1‑234 (MQ=255)
                                                                                                                                                                             tgcacacaTGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCg                                                                                            <  1:43088/206‑1 (MQ=255)
                                                                                                                                                                                   caTGGTTCCCTTCGTGGTTTCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTAACAGCTTTGGTg                                                       <  1:320195/239‑1 (MQ=255)
                                                                                                                                                                                   caTGGTTCCCTTCGTGGTATCGGGCGGTAGTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTAACAGCTTTGGTg                                                       <  1:634419/239‑1 (MQ=255)
                                                                                                                                                                                     tGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTAACAGCTTTGGTGCg                                                     <  1:321206/239‑1 (MQ=255)
                                                                                                                                                                                     tGATTCCCTTCGGGGTAGCGGGCGGTAGTGTGGGCGCGGTTTCCGGCATGTGGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGAGGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTAACAGCTTTGGTGCg                                                     <  1:30362/239‑1 (MQ=255)
                                                                                                                                                                                        aTTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTc                                                                                                                                            <  1:463325/150‑1 (MQ=255)
                                                                                                                                                                                           cccTTCGTGGTATCTGGCGGTATTTTGCTGGCGGTTTCCGTCAGGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTAACAGCTTTGGTGCGGGctt                                                <  1:494871/239‑1 (MQ=255)
                                                                                                                                                                                             cTTCGGGGTATCGGGGGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTAACAGCTTTGGTGCGGGcttct                                              <  1:253274/239‑1 (MQ=255)
                                                                                                                                                                                              ttCTTGGTTTGGGGGTGTTTTTTGGTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGGGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTAACAGCTTTGGTGCGGGcttctt                                             <  2:481194/239‑1 (MQ=255)
                                                                                                                                                                                                cGTGGTTTCGGGCGGTATTTTGCGGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGt                                                                        <  1:100822/210‑1 (MQ=255)
                                                                                                                                                                                                                          cGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTAACAGCTTTGGTGCGGGCTTCTTTGGTGCACTGATCGCCGGGATTATcgg                  <  1:565967/238‑1 (MQ=255)
                                                                                                                                                                                                                               ttCCGTCATTTTTTATGGCAAAGGCGCCGGGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCTTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTAACAGCTTTGGTGCGGGCTTCTTTGGTGCACTGATCGCCGGGATTATCGGCGGCa             <  2:236577/239‑1 (MQ=255)
                                                                                                                                                                                                                                   gTCATGTTGAATGGCAAAGGCGCAGTGCCGAATGCCGTAGCCGTTCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTAACAGCTTTGGTGCGGGCTTCTTTGGTGCACTGATCGCCGGGATTATCGGCGGCATCgt         <  1:440935/239‑1 (MQ=255)
                                                                                                                                                                                                                                   gTCAGGTTGTATGGCAAAGGGGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTAACAGCTTTGGTGCGGGCTTCTTTGGTGCACTGATCGCCGGGATTATCGGCGGCATCgt         <  1:889540/239‑1 (MQ=255)
                                                                                                                                                                                                                                          tgtATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTAACAGCTTTGGTGCGGGCTTCTTTGGTGCACTGATCGCCGGGATTATCGGCGGCATCGTGGTGCAt  <  1:687644/239‑1 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                 
ATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTAACAGCTTTGGTGCGGGCTTCTTTGGTGCACTGATCGCCGGGATTATCGGCGGCATCGTGGTGCAT  >  NC_000913/4142412‑4142882

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