Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,959,995 A→T 6.2% noncoding (74/103 nt) PP_mr63 ← rnk_pseudo ncRNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,959,9950AT 6.2% 239.8 / 3.4 80noncoding (74/103 nt)PP_mr63rnk_pseudo ncRNA
Reads supporting (aligned to +/- strand):  ref base A (39/36);  new base T (3/2);  total (42/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

TCGTCGAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACT‑ACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGCCGGCGCCGTTC  >  NC_002947/5959856‑5960141
                                                                                                                                              |                                                                                                                                                   
tcgtcgAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCTGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCgggg                                                                                                                                            >  6:125068/1‑149 (MQ=255)
tcgtcgAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCgggg                                                                                                                                            >  6:329385/1‑149 (MQ=255)
tcgtcgAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCgggg                                                                                                                                            >  1:96944/1‑149 (MQ=255)
tcgtcgAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCgg                                                                                                                                              <  4:28997/147‑1 (MQ=255)
tcgtcgAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAAACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCgggg                                                                                                                                            >  5:323535/1‑149 (MQ=255)
  gtcgAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCTTGGGCGGGGt                                                                                                                                           <  3:136247/148‑1 (MQ=255)
  gtcgAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGt                                                                                                                                           <  6:221761/148‑1 (MQ=255)
  gtcgAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAAACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGtt                                                                                                                                          <  6:323535/149‑1 (MQ=255)
   tcgAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGttt                                                                                                                                         >  7:281907/1‑149 (MQ=255)
     gAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTga                                                                                                                                       <  1:316466/149‑1 (MQ=255)
      aGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCTTGGGCGGGGttttgc                                                                                                                                      >  4:40565/1‑146 (MQ=255)
       ggCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTgaga                                                                                                                                     <  6:105610/149‑1 (MQ=255)
       ggCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTgag                                                                                                                                      >  1:67191/1‑148 (MQ=255)
         cTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAg                                                                                                                                   <  5:272032/149‑1 (MQ=255)
               aTCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTg                                                                                                                             <  1:65171/149‑1 (MQ=255)
                  aGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGc                                                                                                                           <  7:48802/148‑1 (MQ=255)
                  aGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGc                                                                                                                           <  2:96944/148‑1 (MQ=255)
                  aGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGc                                                                                                                           <  3:221805/148‑1 (MQ=255)
                  aGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGc                                                                                                                           <  8:281907/148‑1 (MQ=255)
                        tCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCAcc                                                                                                                     >  5:42422/1‑148 (MQ=255)
                        tCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCAcc                                                                                                                     >  4:140410/1‑148 (MQ=255)
                        tCCAGACGCTGTAAGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACACCGCCCAACGGCGGGGTTTGAGAAGGCTTTGGCGTCa                                                                                                                       >  5:351550/1‑146 (MQ=255)
                         ccAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCCCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGACGGCTTTGGCGTCACCgc                                                                                                                   >  6:165864/1‑148 (MQ=255)
                           aGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCTTGGGCGGGGtttt                                                                                                                                        <  1:15918/126‑2 (MQ=255)
                           aGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCTTGGGCGGGGtttt                                                                                                                                        >  2:15918/1‑125 (MQ=255)
                           aGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTg                                                                                                                             <  2:266038/137‑1 (MQ=255)
                           aGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTg                                                                                                                             >  1:266038/1‑137 (MQ=255)
                              cGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCa                                                                                                               <  2:67191/148‑1 (MQ=255)
                                   tACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAAccc                                                                                                         <  3:155572/149‑1 (MQ=255)
                                     cGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGa                                                                                                       <  7:145171/149‑1 (MQ=255)
                                     cGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGa                                                                                                       >  5:307300/1‑149 (MQ=255)
                                      gTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGa                                                                                                       >  1:326510/1‑148 (MQ=255)
                                        ccAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGACCCCGaca                                                                                                     >  4:260677/1‑148 (MQ=255)
                                          aaTCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCggg                                                                                                                                             <  5:71304/106‑1 (MQ=255)
                                          aaTCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCggg                                                                                                                                             >  6:71304/1‑106 (MQ=255)
                                          aaTCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACT‑a                                                                                                 <  2:326510/149‑1 (MQ=255)
                                                       tGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTgg                                                                                                                            >  1:174800/1‑110 (MQ=255)
                                                           ggAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACT‑ACCACAGCCGGGCAAATa                                                                                >  2:230840/1‑149 (MQ=255)
                                                             aaGGCTTTGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTTCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAATGCTTTGTCGTCACCGACGCATAACCCGACACT‑ACCAAAGCCGGGCAAATACt                                                                              <  2:190331/149‑1 (MQ=255)
                                                                      tGCTCATGATCCAGTCAGACTCCTTTGAACAGGCTGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACT‑ACCACAGCCGGGCAAATACTCAAGCCCCt                                                                     <  5:125068/149‑1 (MQ=255)
                                                                      tGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACT‑ACCACAGCCGGGCAAATa                                                                                <  3:151625/138‑1 (MQ=255)
                                                                      tGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACT‑ACCACAGCCGGGCAAATa                                                                                >  4:151625/1‑138 (MQ=255)
                                                                        cTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACT‑ACCACAGCCGGGCAAATACTCAAGCCCCTgc                                                                   <  5:329385/149‑1 (MQ=255)
                                                                           aTGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACT‑ACCACAGCCGGGCAAATACTCAAGCCCCTGCGCt                                                                <  3:140410/149‑1 (MQ=255)
                                                                              aTCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACT‑ACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAg                                                             >  4:152457/1‑149 (MQ=255)
                                                                               tCCATTCAGCCTCCTTTTAACAGGCAGTTTTCAATGTGCAAAAAAACCAAA‑‑‑ACCCCGCCCAAGGGGGGGGTTTGATAAGGCTTTGGCGTCACCGACGCAGAACCCGACACT‑ACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAgt                                                            <  3:260677/149‑1 (MQ=255)
                                                                               tCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTCACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCa                                                              <  1:230840/148‑1 (MQ=255)
                                                                                 cAGTCAGACTCCTTTGAACAGTCAGTTTTCAATGTGCAAATAAAGCCTATTTTCCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACT‑ACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCa                                                              <  1:283839/148‑1 (MQ=255)
                                                                                  aGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACT‑ACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGtgtg                                                         >  8:230787/1‑149 (MQ=255)
                                                                                        aCTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGc                                                                                                                                                <  6:272615/57‑1 (MQ=255)
                                                                                        aCTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGc                                                                                                                                                >  5:272615/1‑57 (MQ=255)
                                                                                          tCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAGGGTTTGGCGTCACCGACGCAGAACCCGACACT‑ACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGcc                                                 >  3:282329/1‑149 (MQ=255)
                                                                                            cTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACT‑ACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTg                                               <  8:196960/149‑1 (MQ=255)
                                                                                             tttGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACT‑ACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGc                                              >  4:129624/1‑149 (MQ=255)
                                                                                                 aaCAGGCAGTTTTCAATGTGCAAATAACGCAAA‑‑‑CCCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACT‑ACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTcgc                                          <  5:165864/149‑1 (MQ=255)
                                                                                                 aaCAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACT‑ACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCa                                                              <  3:239550/129‑1 (MQ=255)
                                                                                                 aaCAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACT‑ACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCa                                                              >  4:239550/1‑129 (MQ=255)
                                                                                                     ggCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACT‑ACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCGTAGCCTGCTCgcgc                                        >  4:63533/1‑147 (MQ=255)
                                                                                                      gCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACT‑ACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAg                                     >  6:338526/1‑149 (MQ=255)
                                                                                                         gTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACT‑ACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGAcc                                  >  1:260976/1‑149 (MQ=255)
                                                                                                           tttCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACT‑ACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTg                                >  3:254405/1‑149 (MQ=255)
                                                                                                               aaTGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACAGACGCAGCACCCGACACT‑ACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGAACTGGGc                             >  3:127273/1‑148 (MQ=255)
                                                                                                                  gtgCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACT‑AGCACAGCCGAACAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTcg                         >  1:110719/1‑149 (MQ=255)
                                                                                                                   tgCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACT‑ACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTcgc                        <  3:129624/149‑1 (MQ=255)
                                                                                                                     cAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACT‑ACCACAGCCGGGCAATTACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTcgcg                       <  1:270925/148‑1 (MQ=255)
                                                                                                                     cAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACT‑ACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCgg                      >  1:308479/1‑149 (MQ=255)
                                                                                                                      aaaTAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACT‑ACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCgg                      >  4:333438/1‑148 (MQ=255)
                                                                                                                       aaTAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACT‑ACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGcc                                                 >  3:181451/1‑120 (MQ=255)
                                                                                                                       aaTAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACT‑ACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGcc                                                 <  4:181451/120‑1 (MQ=255)
                                                                                                                       aaTAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACT‑ACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCa                    <  3:152457/149‑1 (MQ=255)
                                                                                                                       aaTAAAGCAAA‑‑‑ACCCCGCCCAAGGGAGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACT‑ACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCgg                      >  6:166173/1‑147 (MQ=255)
                                                                                                                         tAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACT‑ACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGGGTGCCTTAGCCTGCTCGCGCAAGCCCTCGGCCTCGCGGCAg                   >  8:142502/1‑148 (MQ=255)
                                                                                                                            aGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAGTCTTTGTCGTCACCGACGCAGAACCCGACACC‑ACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCt               <  6:351550/149‑1 (MQ=255)
                                                                                                                             gCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACT‑ACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTg              >  6:48641/1‑149 (MQ=255)
                                                                                                                               caa‑‑‑aCCCCGCCCTTGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACT‑ACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGcc            >  5:309275/2‑149 (MQ=255)
                                                                                                                               aaa‑‑‑aCCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACT‑ACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGcc            >  4:240793/1‑149 (MQ=255)
                                                                                                                                     aCCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACT‑ACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGCCGgc         >  7:187858/1‑149 (MQ=255)
                                                                                                                                       cccGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACT‑ACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGCCGgcgc       >  8:127800/1‑149 (MQ=255)
                                                                                                                                          gCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACT‑ACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGCCGGCGCCGt    <  3:240793/149‑1 (MQ=255)
                                                                                                                                            ccAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACT‑ACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGCCGGCGCCGTTc  <  4:254405/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
TCGTCGAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAA‑‑‑ACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACT‑ACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGCCGGCGCCGTTC  >  NC_002947/5959856‑5960141

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: ‑