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

TGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAGCGCCCAACCATGGTAC  >  NC_000913/813427‑813712
                                                                                                                                                     |                                                                                                                                         
tgatgaACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCa                                                                                                                                           <  3:281536/149‑1 (MQ=255)
tgatgaACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCa                                                                                                                                           <  3:287700/149‑1 (MQ=255)
    gaACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCa                                                                                                                                       >  3:632224/1‑149 (MQ=255)
      aCTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCaga                                                                                                                                     <  3:564302/149‑1 (MQ=255)
        tGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAGTTCCGCTTTTAAACCTTCTGGCTGATC‑GCCagaag                                                                                                                                   >  4:588507/1‑147 (MQ=255)
        tGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAgg                                                                                                                                   <  4:71269/149‑1 (MQ=255)
         gTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCTGATC‑GCCagaagg                                                                                                                                  <  3:588507/149‑4 (MQ=255)
          tttttttATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGc                                                                                                                                  >  1:82152/1‑148 (MQ=255)
           ttttttATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCTGATC‑GCCagaaggt                                                                                                                                 <  1:122792/148‑5 (MQ=255)
                 aTCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGa                                                                                                                           <  1:737547/148‑1 (MQ=255)
                 aTCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGa                                                                                                                           <  1:226952/148‑1 (MQ=255)
                 aTCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGAc                                                                                                                          <  4:811450/149‑1 (MQ=255)
                  tCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCaga                                                                                                                                     <  3:614254/137‑1 (MQ=255)
                  tCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCaga                                                                                                                                     >  4:614254/1‑137 (MQ=255)
                       tataATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCTGATC‑GCCagaagg                                                                                                                                  >  2:764715/1‑132 (MQ=255)
                       tataATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCTGATC‑GCCagaagg                                                                                                                                  <  1:764715/135‑4 (MQ=255)
                       tataATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCg                                                                                                                    >  3:539132/1‑149 (MQ=255)
                                ttgGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGc                                                                                                           >  4:453547/1‑149 (MQ=255)
                                 tgGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCt                                                                                                          <  4:632224/149‑1 (MQ=255)
                                     aTAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGaa                                                                                                      >  4:6589/1‑149 (MQ=255)
                                         ttAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATg                                                                                                  <  3:6589/149‑1 (MQ=255)
                                           aaGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATgg                                                                                                 <  1:377277/148‑1 (MQ=255)
                                              tACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTg                                                                                             <  3:791497/149‑1 (MQ=255)
                                                 gacgaGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATg                                                                                                  <  3:1395/141‑1 (MQ=255)
                                                 gacgaGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATg                                                                                                  >  4:1395/1‑141 (MQ=255)
                                                  acgaGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGgcg                                                                                          >  1:682786/1‑148 (MQ=255)
                                                    gaGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCAc                                                                                       <  4:539132/149‑1 (MQ=255)
                                                      gTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCa                                                                                     >  4:187265/1‑149 (MQ=255)
                                                              aCATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTAc                                                                             <  3:330241/149‑1 (MQ=255)
                                                                 tACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACg                                                                                 <  1:602415/142‑1 (MQ=255)
                                                                 tACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACg                                                                                 >  2:602415/1‑142 (MQ=255)
                                                                 tACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACtt                                                                           <  2:682786/148‑1 (MQ=255)
                                                                 tACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACtt                                                                           >  2:778680/1‑148 (MQ=255)
                                                                  aCCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAGACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGtct                                                                                                                      <  4:571693/104‑1 (MQ=255)
                                                                  aCCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAGACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGtct                                                                                                                      >  3:571693/1‑104 (MQ=255)
                                                                   ccTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAggg                                                                                                             <  1:598746/112‑1 (MQ=255)
                                                                   ccTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAggg                                                                                                             >  2:598746/1‑112 (MQ=255)
                                                                           cccAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATg                                                                                                  <  2:902367/115‑1 (MQ=255)
                                                                           cccAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATg                                                                                                  >  1:902367/1‑115 (MQ=255)
                                                                             cAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGtc                                                                                                                       <  2:366606/92‑1 (MQ=255)
                                                                             cAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGtc                                                                                                                       >  1:366606/1‑92 (MQ=255)
                                                                               aCTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTc                                                             <  2:82152/148‑1 (MQ=255)
                                                                                cTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTg                                                                     <  2:21530/139‑1 (MQ=255)
                                                                                cTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTg                                                                     >  1:21530/1‑139 (MQ=255)
                                                                                cTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAg                                                           >  3:179508/1‑149 (MQ=255)
                                                                                 tCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAgg                                                          >  4:41254/1‑149 (MQ=255)
                                                                                 tCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAg                                                           >  2:722860/1‑148 (MQ=255)
                                                                                    ttCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTgg                                                                         <  1:588395/131‑1 (MQ=255)
                                                                                    ttCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTgg                                                                         >  2:588395/1‑131 (MQ=255)
                                                                                        ggTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGa                                                                                                    >  3:746154/1‑100 (MQ=255)
                                                                                        ggTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGa                                                                                                    <  4:746154/100‑1 (MQ=255)
                                                                                          tAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCt                                                  <  1:722860/148‑1 (MQ=255)
                                                                                           aGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTc                                                >  4:200208/1‑149 (MQ=255)
                                                                                            gCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTc                                                <  2:655389/148‑1 (MQ=255)
                                                                                             cGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGa                                                        <  1:492361/139‑1 (MQ=255)
                                                                                             cGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGa                                                        >  2:492361/1‑139 (MQ=255)
                                                                                                  cATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCAtt                                          >  1:369222/1‑148 (MQ=255)
                                                                                                  cATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCAtt                                          >  2:178436/1‑148 (MQ=255)
                                                                                                   aTGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTg                                         <  1:5693/148‑1 (MQ=255)
                                                                                                   aTGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGc                                        <  3:453547/149‑1 (MQ=255)
                                                                                                    tGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGc                                        >  1:523868/1‑148 (MQ=255)
                                                                                                       gTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCTGATC‑GCCagaag                                                                                                                                   >  2:900642/1‑52 (MQ=25)
                                                                                                       gTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCTGATC‑GCCagaag                                                                                                                                   <  1:900642/54‑3 (MQ=25)
                                                                                                        tAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATg                                   >  4:903411/1‑149 (MQ=255)
                                                                                                          aaCCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGAGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATGGCCAATGTc                                 >  4:666013/1‑149 (MQ=255)
                                                                                                            ccGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCa                                >  1:637775/1‑148 (MQ=255)
                                                                                                             cGTTCAAACTGAATTCCGCTTTTAAACCTTCTGTC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCAt                               <  2:231521/148‑1 (MQ=255)
                                                                                                               ttCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTg                             <  1:178436/148‑1 (MQ=255)
                                                                                                                    aCTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGAcc                       >  3:90111/1‑149 (MQ=255)
                                                                                                                       gAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCtt                     >  1:727026/1‑148 (MQ=255)
                                                                                                                       gAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCtt                     >  1:268145/1‑148 (MQ=255)
                                                                                                                           tCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAgc                 >  2:693703/1‑148 (MQ=255)
                                                                                                                           tCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAgc                 >  2:687080/1‑148 (MQ=255)
                                                                                                                            ccGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAgcg                >  1:84079/1‑148 (MQ=255)
                                                                                                                                    aaaCCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAGCGCCCAACCAt       <  3:200208/149‑1 (MQ=255)
                                                                                                                                      aCCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGACGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAGCGCCCAACCATg      <  2:727026/148‑1 (MQ=255)
                                                                                                                                      aCCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAGCGCCCAACCATgg     <  3:41254/149‑1 (MQ=255)
                                                                                                                                      aCCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAGCGCCCAACCATg      >  2:22563/1‑148 (MQ=255)
                                                                                                                                       ccTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAGCGCCCAACCATGGt    >  4:185779/1‑149 (MQ=255)
                                                                                                                                        cTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAGCGCCCAACCATGGt    >  1:558508/1‑148 (MQ=255)
                                                                                                                                          tCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAGCGCCCAACCATGGTAc  >  2:603513/1‑148 (MQ=255)
                                                                                                                                          tCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAGCGCCCAACCATGGTAc  >  1:596193/1‑148 (MQ=255)
                                                                                                                                                     |                                                                                                                                         
TGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGC‑GATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAGCGCCCAACCATGGTAC  >  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: ‑