Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,799,370 Δ1 bp 9.3% intergenic (‑100/+23) 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,3700A. 9.3% 505.3 / 42.2 151intergenic (‑100/+23)pheM/rplTphenylalanyl‑tRNA synthetase operon leader peptide/50S ribosomal subunit protein L20
Reads supporting (aligned to +/- strand):  ref base A (66/71);  new base . (8/6);  total (74/77)
Fisher's exact test for biased strand distribution p-value = 5.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.25e-01

Only 100 of 152 total aligned reads displayed.

GGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATG  >  NC_000913/1799222‑1799513
                                                                                                                                                    |                                                                                                                                                
 gATTCAGGTGCTAATTTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTGGCGTTACCTTTTGGTACTCTTCTAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGAc                                                                                                                                                 <  4:618029/149‑1 (MQ=255)
   ttCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTa                                                                                                                                               <  2:105945/149‑1 (MQ=255)
   ttCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTa                                                                                                                                               >  3:247707/1‑149 (MQ=255)
      aGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCt                                                                                                                                           <  1:247176/149‑1 (MQ=255)
      aGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCt                                                                                                                                           >  3:729669/1‑149 (MQ=255)
             aaaGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCtt                                                                                                                                    <  4:420493/149‑1 (MQ=255)
                 tAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGGGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCaa                                                                                                                                >  4:954463/1‑149 (MQ=255)
                   aaaaaaGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACt                                                                                                                              <  3:725071/149‑1 (MQ=255)
                   aaaaaaGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACt                                                                                                                              >  4:397384/1‑149 (MQ=255)
                     aaaaGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTgg                                                                                                                            >  4:369741/1‑149 (MQ=255)
                        agaagCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCtt                                                                                                                         <  2:314781/149‑1 (MQ=255)
                           agCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATg                                                                                                                      <  4:142514/149‑1 (MQ=255)
                            gCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGc                                                                                                                     >  2:433080/1‑149 (MQ=255)
                                aaaaTAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCaga                                                                                                                 >  3:1029226/1‑149 (MQ=255)
                                aaaaTAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCaga                                                                                                                 <  4:451317/149‑1 (MQ=255)
                                    tAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGCTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAgctg                                                                                                             >  4:702527/1‑149 (MQ=255)
                                     agcagcATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCATAAGCCagttgaaa                                                                                                            <  1:752676/149‑7 (MQ=255)
                                     agcagcATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCATAAGCCagttgaaa                                                                                                            <  4:853119/149‑7 (MQ=255)
                                         gcaTTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTc                                                                                                        >  1:275860/1‑149 (MQ=255)
                                           attcattGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGAc                                                                                                                                                 >  4:198357/1‑107 (MQ=255)
                                           attcattGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGAc                                                                                                                                                 <  3:198357/107‑1 (MQ=255)
                                           attcattGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGc                                                                                                      <  4:1029226/149‑1 (MQ=255)
                                            ttcattGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCt                                                                                                     >  2:1069090/1‑149 (MQ=255)
                                              cattGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCttt                                                                                                   <  3:810596/149‑1 (MQ=255)
                                               attGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCtttt                                                                                                  <  2:678981/149‑1 (MQ=255)
                                               attGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCGACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCtttt                                                                                                  >  2:944060/1‑149 (MQ=255)
                                                 tgcttgcGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTc                                                                                                <  3:649983/149‑1 (MQ=255)
                                                   cttgcGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCttt                                                                                                   <  1:179724/144‑1 (MQ=255)
                                                      gcGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTc                                                                                                                                  >  1:127198/1‑110 (MQ=255)
                                                        gTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCaga                                                                                         <  2:229270/149‑1 (MQ=255)
                                                        gTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCaga                                                                                         >  1:487288/1‑149 (MQ=255)
                                                        gTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGCTAGTCTCCCTCtt                                                                                                                                    >  4:125940/1‑106 (MQ=255)
                                                             cTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGt                                                                                    <  2:661379/149‑1 (MQ=255)
                                                              ttttGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTg                                                                                   >  2:847942/1‑149 (MQ=255)
                                                              ttttGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTg                                                                                   >  2:635509/1‑149 (MQ=255)
                                                                  ggTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACg                                                                               >  2:355051/1‑149 (MQ=255)
                                                                   gTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGc                                                                              <  3:173481/149‑1 (MQ=255)
                                                                       tctTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACt                                                                          >  4:734290/1‑149 (MQ=255)
                                                                         tttaaaaGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCATTGCTTTGATGTCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTGTCAACTCGCTGATACCAGAGCTGCTATCGCTTTTTCAACCAGAGCTTTGAACGCTACttt                                                                        >  2:168876/4‑149 (MQ=255)
                                                                         ttCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACttt                                                                        <  1:1050489/149‑1 (MQ=255)
                                                                           cAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGt                                                                      <  2:536545/149‑1 (MQ=255)
                                                                           cAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGt                                                                      <  3:235955/149‑1 (MQ=255)
                                                                           cAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGt                                                                      >  3:235159/1‑149 (MQ=255)
                                                                            aaaaGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTc                                                                     >  1:365343/1‑149 (MQ=255)
                                                                            aaaaGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTc                                                                     <  1:6209/149‑1 (MQ=255)
                                                                            aaaaGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTc                                                                     <  3:1053410/149‑1 (MQ=255)
                                                                              aaGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAGAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGa                                                                   >  4:1031548/1‑149 (MQ=255)
                                                                                 aCCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATa                                                                >  3:59228/1‑149 (MQ=255)
                                                                                   cTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACt                                                              >  3:864537/1‑149 (MQ=255)
                                                                                    tttATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTc                                                                     <  2:438122/141‑1 (MQ=255)
                                                                                    tttATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTg                                                             >  2:215636/1‑149 (MQ=255)
                                                                                    tttATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTg                                                             <  1:355051/149‑1 (MQ=255)
                                                                                     ttATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGc                                                            <  3:1031548/149‑1 (MQ=255)
                                                                                     ttATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGc                                                            <  1:433080/149‑1 (MQ=255)
                                                                                         tGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGCTAGTCTCCctc                                                                                                                                      <  1:942037/71‑1 (MQ=38)
                                                                                          gAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGatat                                                       >  1:83141/1‑149 (MQ=255)
                                                                                              aGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGc                                                   <  4:864537/149‑1 (MQ=255)
                                                                                                gCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCa                                                 >  4:109234/1‑149 (MQ=255)
                                                                                                    cGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGAt                                             <  2:365343/149‑1 (MQ=255)
                                                                                                    cGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGAt                                             <  3:397384/149‑1 (MQ=255)
                                                                                                       cGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCtt                                          <  3:994122/149‑1 (MQ=255)
                                                                                                        gATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTa                                         <  2:122031/149‑1 (MQ=255)
                                                                                                        gATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTa                                         <  3:369741/149‑1 (MQ=255)
                                                                                                           aaaaGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGa                                                                                  >  4:871241/1‑105 (MQ=255)
                                                                                                           aaaaGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGa                                                                                  <  3:871241/105‑1 (MQ=255)
                                                                                                               gTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGa                                  <  3:830994/149‑1 (MQ=255)
                                                                                                                  aaTGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGAttt                               <  4:59228/149‑1 (MQ=255)
                                                                                                                    tGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAAcc                                                                                            <  3:1002172/86‑1 (MQ=255)
                                                                                                                      ttttGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAAc                           >  3:398435/1‑149 (MQ=255)
                                                                                                                      ttttGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAAc                           >  3:477896/1‑149 (MQ=255)
                                                                                                                         tGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACaga                        <  1:97773/149‑1 (MQ=255)
                                                                                                                           aTGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCCACAGAgg                      >  2:375672/1‑149 (MQ=255)
                                                                                                                           aTGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAgg                      <  4:536405/149‑1 (MQ=255)
                                                                                                                             ggCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGa                                                          >  4:49448/1‑111 (MQ=255)
                                                                                                                             ggCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGa                                                          <  3:49448/111‑1 (MQ=255)
                                                                                                                              gCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGcc                                                  >  1:957482/1‑118 (MQ=255)
                                                                                                                              gCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCCTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCtt                   >  1:837831/1‑149 (MQ=255)
                                                                                                                                 ttGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCttttt                >  1:328611/1‑149 (MQ=255)
                                                                                                                                 ttGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCttttt                >  1:176893/1‑149 (MQ=255)
                                                                                                                                 ttGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCttttt                >  4:32791/1‑149 (MQ=255)
                                                                                                                                   gAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTc              >  4:819231/1‑149 (MQ=255)
                                                                                                                                    aaaCGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCa             >  3:675407/1‑149 (MQ=255)
                                                                                                                                     aaCGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAg            <  1:749237/149‑1 (MQ=255)
                                                                                                                                     aaCGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAg            <  2:429943/149‑1 (MQ=255)
                                                                                                                                     aaCGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAg            <  3:317480/149‑1 (MQ=255)
                                                                                                                                     aaCGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAACACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAg            <  3:322934/149‑1 (MQ=255)
                                                                                                                                      aCGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAgg           >  4:573210/1‑149 (MQ=255)
                                                                                                                                      aCGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTCTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAgg           >  2:608961/1‑149 (MQ=255)
                                                                                                                                        gAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCa                                                 >  4:28041/1‑109 (MQ=255)
                                                                                                                                        gAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCa                                                 <  3:28041/109‑1 (MQ=255)
                                                                                                                                        gAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGcc         <  3:573210/149‑1 (MQ=255)
                                                                                                                                          aaaGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCAt       >  1:465172/1‑149 (MQ=255)
                                                                                                                                          aaaGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTCCAACAGAGGCTTTTTTCAGGCCAt       >  3:328807/1‑149 (MQ=255)
                                                                                                                                           aaGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCAtt      <  3:141976/149‑1 (MQ=255)
                                                                                                                                           aaGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCAtt      <  1:931357/149‑1 (MQ=255)
                                                                                                                                           aaGAGGGAGCTAGTCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGa                                  <  2:948881/121‑1 (MQ=255)
                                                                                                                                            agagGGAGCTAGTCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTc                                                                     <  2:168692/85‑1 (MQ=39)
                                                                                                                                            agagGGAGCTAGTCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATtg     >  1:534206/1‑149 (MQ=255)
                                                                                                                                            agagGGAGCTAGTCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATtg     >  4:754246/1‑149 (MQ=255)
                                                                                                                                              agGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATtgat   <  2:811269/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATG  >  NC_000913/1799222‑1799513

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