Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 813,575 Δ1 bp 8.7% coding (50/2022 nt) uvrB → exision nuclease of nucleotide excision repair, DNA damage recognition component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913813,5750A. 8.7% 278.7 / 19.0 80coding (50/2022 nt)uvrBexision nuclease of nucleotide excision repair, DNA damage recognition component
Reads supporting (aligned to +/- strand):  ref base A (40/33);  new base . (3/4);  total (43/37)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

TGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAGCGCCCAACCATGGTAC  >  NC_000913/813427‑813712
                                                                                                                                                     |                                                                                                                                         
tgatgaACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCa                                                                                                                                           <  3:281536/149‑1 (MQ=255)
tgatgaACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCa                                                                                                                                           <  3:287700/149‑1 (MQ=255)
    gaACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCa                                                                                                                                       >  3:632224/1‑149 (MQ=255)
      aCTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCaga                                                                                                                                     <  3:564302/149‑1 (MQ=255)
        tGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAGTTCCGCTTTTAAACCTTCTGGCTGATCGCCagaag                                                                                                                                   >  4:588507/1‑147 (MQ=255)
        tGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAgg                                                                                                                                   <  4:71269/149‑1 (MQ=255)
         gTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCTGATCGCCagaagg                                                                                                                                  <  3:588507/149‑4 (MQ=255)
          tttttttATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGc                                                                                                                                  >  1:82152/1‑148 (MQ=255)
           ttttttATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCTGATCGCCagaaggt                                                                                                                                 <  1:122792/148‑5 (MQ=255)
                 aTCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGa                                                                                                                           <  1:737547/148‑1 (MQ=255)
                 aTCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGa                                                                                                                           <  1:226952/148‑1 (MQ=255)
                 aTCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGAc                                                                                                                          <  4:811450/149‑1 (MQ=255)
                  tCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCaga                                                                                                                                     <  3:614254/137‑1 (MQ=255)
                  tCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCaga                                                                                                                                     >  4:614254/1‑137 (MQ=255)
                       tataATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCTGATCGCCagaagg                                                                                                                                  >  2:764715/1‑132 (MQ=255)
                       tataATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCTGATCGCCagaagg                                                                                                                                  <  1:764715/135‑4 (MQ=255)
                       tataATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCg                                                                                                                    >  3:539132/1‑149 (MQ=255)
                                ttgGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGc                                                                                                           >  4:453547/1‑149 (MQ=255)
                                 tgGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCt                                                                                                          <  4:632224/149‑1 (MQ=255)
                                     aTAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGaa                                                                                                      >  4:6589/1‑149 (MQ=255)
                                         ttAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATg                                                                                                  <  3:6589/149‑1 (MQ=255)
                                           aaGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATgg                                                                                                 <  1:377277/148‑1 (MQ=255)
                                              tACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTg                                                                                             <  3:791497/149‑1 (MQ=255)
                                                 gacgaGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATg                                                                                                  <  3:1395/141‑1 (MQ=255)
                                                 gacgaGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATg                                                                                                  >  4:1395/1‑141 (MQ=255)
                                                  acgaGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGgcg                                                                                          >  1:682786/1‑148 (MQ=255)
                                                    gaGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCAc                                                                                       <  4:539132/149‑1 (MQ=255)
                                                      gTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCa                                                                                     >  4:187265/1‑149 (MQ=255)
                                                              aCATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTAc                                                                             <  3:330241/149‑1 (MQ=255)
                                                                 tACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACg                                                                                 <  1:602415/142‑1 (MQ=255)
                                                                 tACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACg                                                                                 >  2:602415/1‑142 (MQ=255)
                                                                 tACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACtt                                                                           <  2:682786/148‑1 (MQ=255)
                                                                 tACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACtt                                                                           >  2:778680/1‑148 (MQ=255)
                                                                  aCCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAGACCTTCTGGCGATCAGCCAGAGGCGATTCGACGtct                                                                                                                      <  4:571693/104‑1 (MQ=255)
                                                                  aCCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAGACCTTCTGGCGATCAGCCAGAGGCGATTCGACGtct                                                                                                                      >  3:571693/1‑104 (MQ=255)
                                                                   ccTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAggg                                                                                                             <  1:598746/112‑1 (MQ=255)
                                                                   ccTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAggg                                                                                                             >  2:598746/1‑112 (MQ=255)
                                                                           cccAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATg                                                                                                  <  2:902367/115‑1 (MQ=255)
                                                                           cccAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATg                                                                                                  >  1:902367/1‑115 (MQ=255)
                                                                             cAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGtc                                                                                                                       <  2:366606/92‑1 (MQ=255)
                                                                             cAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGtc                                                                                                                       >  1:366606/1‑92 (MQ=255)
                                                                               aCTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTc                                                             <  2:82152/148‑1 (MQ=255)
                                                                                cTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTg                                                                     <  2:21530/139‑1 (MQ=255)
                                                                                cTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTg                                                                     >  1:21530/1‑139 (MQ=255)
                                                                                cTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAg                                                           >  3:179508/1‑149 (MQ=255)
                                                                                 tCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAgg                                                          >  4:41254/1‑149 (MQ=255)
                                                                                 tCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAg                                                           >  2:722860/1‑148 (MQ=255)
                                                                                    ttCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTgg                                                                         <  1:588395/131‑1 (MQ=255)
                                                                                    ttCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTgg                                                                         >  2:588395/1‑131 (MQ=255)
                                                                                        ggTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGa                                                                                                    >  3:746154/1‑100 (MQ=255)
                                                                                        ggTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGa                                                                                                    <  4:746154/100‑1 (MQ=255)
                                                                                          tAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCt                                                  <  1:722860/148‑1 (MQ=255)
                                                                                           aGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTc                                                >  4:200208/1‑149 (MQ=255)
                                                                                            gCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTc                                                <  2:655389/148‑1 (MQ=255)
                                                                                             cGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGa                                                        <  1:492361/139‑1 (MQ=255)
                                                                                             cGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGa                                                        >  2:492361/1‑139 (MQ=255)
                                                                                                  cATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCAtt                                          >  1:369222/1‑148 (MQ=255)
                                                                                                  cATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCAtt                                          >  2:178436/1‑148 (MQ=255)
                                                                                                   aTGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTg                                         <  1:5693/148‑1 (MQ=255)
                                                                                                   aTGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGc                                        <  3:453547/149‑1 (MQ=255)
                                                                                                    tGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGc                                        >  1:523868/1‑148 (MQ=255)
                                                                                                       gTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCTGATCGCCagaag                                                                                                                                   >  2:900642/1‑52 (MQ=25)
                                                                                                       gTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCTGATCGCCagaag                                                                                                                                   <  1:900642/54‑3 (MQ=25)
                                                                                                        tAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATg                                   >  4:903411/1‑149 (MQ=255)
                                                                                                          aaCCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGAGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATGGCCAATGTc                                 >  4:666013/1‑149 (MQ=255)
                                                                                                            ccGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCa                                >  1:637775/1‑148 (MQ=255)
                                                                                                             cGTTCAAACTGAATTCCGCTTTTAAACCTTCTGTCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCAt                               <  2:231521/148‑1 (MQ=255)
                                                                                                               ttCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTg                             <  1:178436/148‑1 (MQ=255)
                                                                                                                    aCTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGAcc                       >  3:90111/1‑149 (MQ=255)
                                                                                                                       gAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCtt                     >  1:727026/1‑148 (MQ=255)
                                                                                                                       gAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCtt                     >  1:268145/1‑148 (MQ=255)
                                                                                                                           tCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAgc                 >  2:693703/1‑148 (MQ=255)
                                                                                                                           tCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAgc                 >  2:687080/1‑148 (MQ=255)
                                                                                                                            ccGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAgcg                >  1:84079/1‑148 (MQ=255)
                                                                                                                                    aaaCCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAGCGCCCAACCAt       <  3:200208/149‑1 (MQ=255)
                                                                                                                                      aCCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGACGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAGCGCCCAACCATg      <  2:727026/148‑1 (MQ=255)
                                                                                                                                      aCCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAGCGCCCAACCATgg     <  3:41254/149‑1 (MQ=255)
                                                                                                                                      aCCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAGCGCCCAACCATg      >  2:22563/1‑148 (MQ=255)
                                                                                                                                       ccTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAGCGCCCAACCATGGt    >  4:185779/1‑149 (MQ=255)
                                                                                                                                        cTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAGCGCCCAACCATGGt    >  1:558508/1‑148 (MQ=255)
                                                                                                                                          tCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAGCGCCCAACCATGGTAc  >  2:603513/1‑148 (MQ=255)
                                                                                                                                          tCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAGCGCCCAACCATGGTAc  >  1:596193/1‑148 (MQ=255)
                                                                                                                                                     |                                                                                                                                         
TGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAGCGCCCAACCATGGTAC  >  NC_000913/813427‑813712

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