Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 59,301 A→C 5.2% intergenic (+22/+386) yabP → / ← rluA pseudogene, pentapeptide repeats‑containing/dual specificity 23S rRNA pseudouridine(746), tRNA pseudouridine(32) synthase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091359,3010AC 5.2% 419.0 / 8.1 144intergenic (+22/+386)yabP/rluApseudogene, pentapeptide repeats‑containing/dual specificity 23S rRNA pseudouridine(746), tRNA pseudouridine(32) synthase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (65/70);  new base C (6/2);  total (72/72)
Fisher's exact test for biased strand distribution p-value = 1.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 144 total aligned reads displayed.

ATTATTATTAATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATGGGCATCGCCAACCAGACTGCTAATTAGC  >  NC_000913/59154‑59447
                                                                                                                                                   |                                                                                                                                                  
    ttattaATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGccc                                                                                                                                               <  5:56856/149‑1 (MQ=255)
    ttattaATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGccc                                                                                                                                               >  6:54556/1‑149 (MQ=255)
    ttattaATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGccc                                                                                                                                               <  1:57075/149‑1 (MQ=255)
      attaATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTTAAGAAGTTGCTACAAAAAAAGTACTGCTTAACAAATTACTGAAAACGCAATTATTGAATGATCCAGAATAAGCAAAGGTTGAAGGGGCTGGAACGcccgt                                                                                                                                             >  5:230950/1‑147 (MQ=255)
       ttaATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCtt                                                                                                                                            <  7:257033/149‑1 (MQ=255)
           taaAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAAc                                                                                                                                        >  7:351427/1‑149 (MQ=255)
               aTCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCGTTCCAGCCCCTTCAACCTTa                                                                                                                                    <  3:268554/149‑1 (MQ=39)
                tCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAg                                                                                                                                   <  5:90065/149‑1 (MQ=255)
                  ccccccAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCACAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCa                                                                                                                                 <  2:395959/149‑1 (MQ=255)
                     cccAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACAAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTc                                                                                                                                           >  5:361082/1‑136 (MQ=255)
                          aaaaaTTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTa                                                                                                                              >  4:147527/1‑144 (MQ=255)
                          aaaaaTTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTg                                                                                                                         <  5:54556/149‑1 (MQ=255)
                            aaaTTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTggg                                                                                                                       <  8:285199/149‑1 (MQ=255)
                              aTTGATACTAAAGGAGTTGCTACTGAAGCAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACGACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGgat                                                                                                                     >  5:121761/1‑149 (MQ=255)
                              aTTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAAc                                                                                                                                        >  1:333427/1‑130 (MQ=255)
                              aTTGAGACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACGAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAGCGCCCCTTCAACCTGAGCAGTAGCGGGGgat                                                                                                                     >  6:302714/1‑149 (MQ=255)
                                      tAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGAttt                                                                                                                <  2:186610/146‑1 (MQ=255)
                                      tAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGAttt                                                                                                                >  1:186610/1‑146 (MQ=255)
                                      tAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTca                                                                                                              >  8:155347/1‑148 (MQ=255)
                                      tAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTca                                                                                                              >  1:26737/1‑148 (MQ=255)
                                        aaGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGt                                                                                                                          >  3:6290/1‑134 (MQ=255)
                                        aaGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGt                                                                                                                          <  4:6290/134‑1 (MQ=255)
                                        aaGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACaa                                                                                                           >  1:149853/1‑149 (MQ=255)
                                        aaGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTAGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACaa                                                                                                           >  3:394189/1‑149 (MQ=255)
                                          ggAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAAtt                                                                                                         <  3:161010/149‑1 (MQ=255)
                                              ttGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGa                                                                                                      >  3:314255/1‑148 (MQ=255)
                                                  tACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAgcc                                                                                                 >  5:312926/1‑147 (MQ=255)
                                                     tGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTg                                                                                              >  7:108532/1‑149 (MQ=255)
                                                       aagaagTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGaaa                                                                                                    >  8:34607/1‑141 (MQ=255)
                                                         gaagTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGAc                                                                                                 <  7:389561/142‑1 (MQ=255)
                                                            gTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCatgat                                                                                        >  8:363609/1‑148 (MQ=255)
                                                             ttGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCAtgatga                                                                                      <  3:50915/149‑1 (MQ=255)
                                                               gCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGACCGCCCCTTCAACCATAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGc                                                                                    >  8:326121/1‑149 (MQ=255)
                                                               gCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGc                                                                                    >  7:146589/1‑149 (MQ=255)
                                                                 tactaAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAAc                                                                                                                                        <  2:30692/95‑1 (MQ=255)
                                                                 tactaAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGc                                                                                    >  6:437525/1‑147 (MQ=255)
                                                                   ctaAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTaga                                                                                <  3:134788/149‑1 (MQ=255)
                                                                   ctaAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGGTTAAGGGGCTGTAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTa                                                                                  <  5:104666/147‑1 (MQ=255)
                                                                    taAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGaaa                                                                                                    >  1:287214/1‑128 (MQ=255)
                                                                    taAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTagag                                                                               <  7:165257/149‑1 (MQ=255)
                                                                    taAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTa                                                                                  >  7:255897/1‑146 (MQ=255)
                                                                    taAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTGGAAGGGGCTGGAACGACCCTTCCACCTTAGCAGTAGCGTGGGCTGATTTCACAATTAGAAAGACCTGCAGGATGAGCTa                                                                                  >  5:314181/1‑146 (MQ=255)
                                                                     aaaaaaGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTagag                                                                               >  3:251008/1‑148 (MQ=255)
                                                                     aaaaaaGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTagag                                                                               >  1:87401/1‑148 (MQ=255)
                                                                      aaaaaGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAg                                                                                                                                   <  1:140687/95‑1 (MQ=255)
                                                                      aaaaaGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAg                                                                                                                                   >  2:140687/1‑95 (MQ=255)
                                                                      aaaaaGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAgaga                                                                              <  2:271110/148‑1 (MQ=255)
                                                                      aaaaaGTAATGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGagaa                                                                             <  8:282166/149‑1 (MQ=255)
                                                                       aaaaGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGATCCAGAATAAGCTAATGTTTAGAGGGCTGGAACGCCACTTAAACCTTAGATGTAGCGTGGTATGATTTAACAATTAGAAAGACCTGCATGAAGAGCTagag                                                                               >  7:323640/1‑146 (MQ=255)
                                                                           gTACTGCTGAACAAATTACTGTCAACGAAATTATTGAATGATCCAGAATAAGCTTTGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGc                                                                        <  2:158156/149‑1 (MQ=255)
                                                                              ctgctgAACAAATTACTGACAACGCAATTATTGAATGAGCCAGTATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTa                                                                      <  2:26737/148‑1 (MQ=255)
                                                                              ctgctgAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTa                                                                      <  6:325739/148‑1 (MQ=255)
                                                                              ctgctgAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTa                                                                      <  8:146589/148‑1 (MQ=255)
                                                                                 ctgAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGc                                                                        <  1:324732/143‑1 (MQ=255)
                                                                                 ctgAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGc                                                                        >  2:324732/1‑143 (MQ=255)
                                                                                  tgAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGAc                                                                 >  6:367292/1‑149 (MQ=255)
                                                                                    aaCAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAgaga                                                                              <  5:398349/134‑1 (MQ=255)
                                                                                    aaCAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAgaga                                                                              >  6:398349/1‑134 (MQ=255)
                                                                                    aaCAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACg                                                                >  7:237726/1‑148 (MQ=255)
                                                                                      cAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCaa                                                             <  8:108532/149‑1 (MQ=255)
                                                                                      cAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCaa                                                             <  4:314255/149‑1 (MQ=255)
                                                                                       aaaTTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGTATGATTTCACAATTAGAAAGACATGCATGATGAGc                                                                                    <  4:59282/125‑1 (MQ=255)
                                                                                       aaaTTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTa                                                                      <  2:325660/139‑1 (MQ=255)
                                                                                       aaaTTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTa                                                                      >  1:325660/1‑139 (MQ=255)
                                                                                              tGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTcaca                                                                                                            >  8:18313/1‑94 (MQ=255)
                                                                                               gACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTg                                                    >  8:395462/1‑149 (MQ=255)
                                                                                                aCAACGCACTTATTGAATGAGCCAGTATCAGCTACGGTGGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGGGGGATGATTTCACAATTAGAAAGACCTGCATGATGCGCTAGAGAAGAGGCTAGTGACGCAAGGcgtcgt                                                     <  6:314181/147‑1 (MQ=255)
                                                                                                 cAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCa                                                  <  6:111957/149‑1 (MQ=255)
                                                                                                 cAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCa                                                  <  2:149853/149‑1 (MQ=255)
                                                                                                 cAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCa                                                  <  2:94585/149‑1 (MQ=255)
                                                                                                  aaCGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAg                                                 <  2:223185/149‑1 (MQ=255)
                                                                                                  aaCGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAg                                                 >  4:266745/1‑149 (MQ=255)
                                                                                                    cGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGa                                               <  6:121761/149‑1 (MQ=255)
                                                                                                     gCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGaa                                                                                                     <  6:266795/94‑1 (MQ=255)
                                                                                                       aaTTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCaa                                                             <  7:170178/132‑1 (MQ=255)
                                                                                                       aaTTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCaa                                                             >  8:170178/1‑132 (MQ=255)
                                                                                                        attattGAATGAGCCAGAATAATCTAAGGTTGTAGGGGCTGGAACGCCCCTTCAACATTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGATAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACg                                           <  5:237508/149‑1 (MQ=255)
                                                                                                            ttGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACTATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGAt                                        <  5:395171/148‑1 (MQ=255)
                                                                                                                 tGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCga                                  <  7:395462/149‑1 (MQ=255)
                                                                                                                   aGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCAcc                                    >  1:195258/1‑145 (MQ=255)
                                                                                                                   aGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCgag                                 <  7:377242/148‑1 (MQ=255)
                                                                                                                    gCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAATAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCgaga                                <  1:140494/148‑1 (MQ=255)
                                                                                                                     ccAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATg                              <  7:363609/149‑1 (MQ=255)
                                                                                                                          tactgCTAAGGTTGAAGGGGCGTTCCAGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATGGGCAt                         >  7:145867/5‑149 (MQ=38)
                                                                                                                          aTAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGCCCTGCATGATGAGCGAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATGGGCAt                         >  3:42272/1‑149 (MQ=255)
                                                                                                                            aaGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATgg                             >  8:125003/1‑143 (MQ=255)
                                                                                                                            aaGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATgg                             <  7:125003/143‑1 (MQ=255)
                                                                                                                               cTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATGGGCATCGCCa                    >  1:351469/1‑149 (MQ=255)
                                                                                                                                tAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATGGGCATCGCCaa                   >  7:29667/1‑149 (MQ=255)
                                                                                                                                tAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATGGGCATCGCCaa                   >  8:412334/1‑149 (MQ=255)
                                                                                                                                tAAGGTTGAAGGGGCGTTCCAGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGc                                                                                    <  5:123809/84‑1 (MQ=21)
                                                                                                                                tAAGGTTGAAGGGGCGTTCCAGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGc                                                                                    >  6:123809/1‑84 (MQ=21)
                                                                                                                                  aGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATGGGCATCGCCAAcc                 <  8:365413/149‑1 (MQ=255)
                                                                                                                                    gTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATGGGCATCGCCAACCa                <  1:351118/148‑1 (MQ=255)
                                                                                                                                    gTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAAGTAGAAAGACCTGCATGATGAGCGAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACCCGGAGCACCGAGCTGGGCATCGCCAAc                  >  2:351118/1‑146 (MQ=255)
                                                                                                                                     ttGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATGGGCATAGCCAACCAGa              >  2:306766/1‑149 (MQ=255)
                                                                                                                                        aaGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATGGGCATCGCCAACCAGACTg           <  2:351469/149‑1 (MQ=255)
                                                                                                                                           gggCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATGGGc                           <  5:80316/130‑1 (MQ=255)
                                                                                                                                                 gAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCatgat                                                                                        >  8:26114/1‑63 (MQ=255)
                                                                                                                                                 gAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATGGGCATCGCCAACCAGACTGCTAATTAGc  >  6:363178/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ATTATTATTAATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATGGGCATCGCCAACCAGACTGCTAATTAGC  >  NC_000913/59154‑59447

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