Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,799,374:1 +T 8.8% intergenic (‑104/+19) pheM ← / ← rplT phenylalanyl‑tRNA synthetase operon leader peptide/50S ribosomal subunit protein L20

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,799,3741.T 8.8% 353.1 / 30.2 91intergenic (‑104/+19)pheM/rplTphenylalanyl‑tRNA synthetase operon leader peptide/50S ribosomal subunit protein L20
Reads supporting (aligned to +/- strand):  ref base . (43/40);  new base T (5/3);  total (48/43)
Fisher's exact test for biased strand distribution p-value = 7.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

TAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGCT  >  NC_000913/1799234‑1799521
                                                                                                                                             |                                                                                                                                                   
tAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCctct                                                                                                                                             >  3:214550/1‑149 (MQ=255)
tAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCctct                                                                                                                                             <  4:479281/149‑1 (MQ=255)
tAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAG‑CTAGTCTCCctct                                                                                                                                             <  4:779837/149‑1 (MQ=255)
 aaaGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAG‑CTAGTCTCCCTCtt                                                                                                                                            <  1:122120/149‑1 (MQ=255)
   aGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTc                                                                                                                                          <  1:21761/149‑1 (MQ=255)
    gTAAAAAAAGAAGCGGAAAGTAGCAGCATTCATTGCTTGTGTTAACTTTTGGTACTCTGCAAAACACCTTTATGGAAAAGGCTAGGGCGATAAAATCCAAGGCTTAGATGGTGTTGAAAAGAAAAGAGGGTGAATAG‑CGCTCTCTTTCa                                                                                                                                         >  2:582210/1‑149 (MQ=255)
         aaaaGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGGACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGACGGCGTTGAATCGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTgg                                                                                                                                    >  3:603464/1‑149 (MQ=255)
                  ggAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCa                                                                                                                           <  3:710328/149‑1 (MQ=255)
                    aaaaTAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCaga                                                                                                                         <  2:778660/149‑1 (MQ=255)
                     aaaTAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCagag                                                                                                                        <  3:809514/149‑1 (MQ=255)
                     aaaTAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCagag                                                                                                                        <  1:655753/149‑1 (MQ=255)
                       aTAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAGA‑AGGGGGACCAG‑GTCCCCCTTTCCACTCGCCTTTGCCAAAgccg                                                                                                                     >  3:561398/1‑147 (MQ=255)
                          gcggcaTTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTTCGGCGATAAAAGTCAGTGTTTTGATGGCGTGGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTATCAACTTACTTATGCCAGAgctgct                                                                                                                   >  2:602380/4‑149 (MQ=255)
                             gcaTTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTc                                                                                                                <  2:373168/149‑1 (MQ=255)
                              cattcattGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAATAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCg                                                                                                               >  4:612630/1‑149 (MQ=255)
                                ttcattGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCt                                                                                                             >  4:732695/1‑149 (MQ=255)
                                  cattGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCttt                                                                                                           <  2:719331/149‑1 (MQ=255)
                                  cattGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCttt                                                                                                           <  4:371550/149‑1 (MQ=255)
                                  cattGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCttt                                                                                                           <  2:536519/149‑1 (MQ=255)
                                   attGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCtttt                                                                                                          >  1:673901/1‑149 (MQ=255)
                                   attGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCCTTCGCtttt                                                                                                          <  4:6350/149‑1 (MQ=255)
                                    ttgcttgcGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCttttt                                                                                                         <  2:306916/149‑1 (MQ=255)
                                        ttgcGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAAc                                                                                                     >  4:234753/1‑149 (MQ=255)
                                        ttgcGTTACCTTTT‑GTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAAcc                                                                                                    >  4:696204/1‑149 (MQ=255)
                                            gTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCaga                                                                                                 >  4:408808/1‑149 (MQ=255)
                                               aCCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCg                                                                                              <  1:574195/149‑1 (MQ=255)
                                                 cTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGt                                                                                            >  3:267960/1‑149 (MQ=255)
                                                 cTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGt                                                                                            >  3:147102/1‑149 (MQ=255)
                                                 cTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTATCGCTTTTTCAACCAGAGCGGt                                                                                            <  1:543624/149‑1 (MQ=255)
                                                  ttttGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCa                                                                                                                           >  3:699876/1‑117 (MQ=255)
                                                  ttttGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCa                                                                                                                           <  4:699876/117‑1 (MQ=255)
                                                  ttttGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTg                                                                                           >  2:467236/1‑149 (MQ=255)
                                                    ttGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGaa                                                                                         <  3:647845/149‑1 (MQ=255)
                                                    ttGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTACAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGaa                                                                                         <  4:409434/149‑1 (MQ=255)
                                                     tGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAAc                                                                                        <  3:170729/149‑1 (MQ=255)
                                                       gTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGc                                                                                      >  3:57751/1‑149 (MQ=255)
                                                         aCTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTa                                                                                    <  1:333101/149‑1 (MQ=255)
                                                         aCTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTa                                                                                    >  4:505789/1‑149 (MQ=255)
                                                          ctctTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTAc                                                                                   >  2:550817/1‑149 (MQ=255)
                                                           tctTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑ATCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAAAGCTACt                                                                                  >  3:242811/1‑149 (MQ=255)
                                                             ttCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACttt                                                                                >  2:649062/1‑149 (MQ=255)
                                                                aaaaGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTc                                                                             >  3:352078/1‑149 (MQ=255)
                                                                aaaaGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTc                                                                             <  4:214550/149‑1 (MQ=255)
                                                                aaaaGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTc                                                                             >  1:748917/1‑149 (MQ=255)
                                                                  aaGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGa                                                                           <  3:118889/149‑1 (MQ=255)
                                                                    gACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAAt                                                                         >  4:679394/1‑149 (MQ=255)
                                                                    gACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAAt                                                                         <  1:122975/149‑1 (MQ=255)
                                                                    gACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAAt                                                                         <  2:673901/149‑1 (MQ=255)
                                                                       cTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTCTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACt                                                                      >  3:690391/1‑149 (MQ=255)
                                                                        tttATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTg                                                                     >  1:367518/1‑149 (MQ=255)
                                                                        tttATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTg                                                                     >  2:726543/1‑149 (MQ=255)
                                                                              gAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTCTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGatat                                                               <  4:690391/149‑1 (MQ=255)
                                                                               aaaaGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATc                                                              <  3:832417/149‑1 (MQ=255)
                                                                               aaaaGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATc                                                              >  1:54884/1‑149 (MQ=255)
                                                                               aaaaAGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATc                                                              >  2:324091/1‑149 (MQ=255)
                                                                                  aGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGc                                                           >  1:709703/1‑149 (MQ=255)
                                                                                          gCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGc                                                                                                                             >  2:528568/1‑75 (MQ=255)
                                                                                          gCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGc                                                                                                                             <  1:528568/75‑1 (MQ=255)
                                                                                          gCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCt                                                   <  3:291489/149‑1 (MQ=255)
                                                                                          gCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCt                                                   >  3:55749/1‑149 (MQ=255)
                                                                                              tAAAAGTCAATGTTTTGATGGCGTTGAACCGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGCCTTATGCCAGAGCAGCCTTCGTTTTGTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACg                                               >  2:639011/1‑149 (MQ=255)
                                                                                              tAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACg                                               >  4:850448/1‑149 (MQ=255)
                                                                                                aaaGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGt                                             >  1:616127/1‑149 (MQ=255)
                                                                                                   ttCAATTTTTTGTTGGCTTTGAAACG‑AAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCGACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCTATATCAGCCAGGATCTTACGCTCGAt                                         <  4:561398/148‑1 (MQ=255)
                                                                                                   gTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGATTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGa                                          <  4:242811/149‑1 (MQ=255)
                                                                                                   gTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGa                                          <  3:421068/149‑1 (MQ=255)
                                                                                                    tCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGAt                                         <  4:147102/149‑1 (MQ=255)
                                                                                                    tCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGAt                                         >  1:602886/1‑149 (MQ=255)
                                                                                                      aaTGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGAttt                                       >  3:750664/1‑149 (MQ=255)
                                                                                                      aaTGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGAttt                                       >  4:640688/1‑149 (MQ=255)
                                                                                                      aaTGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGAttt                                       <  3:234753/149‑1 (MQ=255)
                                                                                                         gTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCaa                                    <  1:550817/149‑1 (MQ=255)
                                                                                                             tGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACaga                                <  3:679394/149‑1 (MQ=255)
                                                                                                              gATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACagag                               <  1:726543/149‑1 (MQ=255)
                                                                                                                       gAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTc                      <  1:649062/149‑1 (MQ=255)
                                                                                                                         aaCGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAg                    <  1:868165/149‑1 (MQ=255)
                                                                                                                          gttgaaaGAGGGAG‑CTAGTCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAgg                   >  3:621820/5‑149 (MQ=255)
                                                                                                                            gAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGcc                 >  1:270066/1‑149 (MQ=255)
                                                                                                                               aaGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACa                                  >  2:338673/1‑129 (MQ=255)
                                                                                                                               aaGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACa                                  <  1:338673/129‑1 (MQ=255)
                                                                                                                               aaGAGGGAG‑CTAGTCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCAtt              >  3:265810/1‑149 (MQ=255)
                                                                                                                               aaGAGGGAG‑CTAGTCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCAtt              >  2:753763/1‑149 (MQ=255)
                                                                                                                                agagGGAG‑CTAGTCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTc                                                                             >  1:635899/1‑85 (MQ=39)
                                                                                                                                agagGGAG‑CTAGTCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTc                                                                             <  2:635899/85‑1 (MQ=39)
                                                                                                                                agagGGAG‑CTAGTCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATtg             >  2:379106/1‑149 (MQ=255)
                                                                                                                                 gagGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATtga            <  4:267960/149‑1 (MQ=255)
                                                                                                                                    ggAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATtgatga         >  2:101700/1‑149 (MQ=255)
                                                                                                                                     gagaCTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTGCGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGaa        >  4:536185/1‑149 (MQ=255)
                                                                                                                                         cTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTg    >  4:5593/1‑149 (MQ=255)
                                                                                                                                          tAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGc   >  1:22530/1‑149 (MQ=255)
                                                                                                                                           aG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGCt  <  3:640688/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
TAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGCT  >  NC_000913/1799234‑1799521

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