Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 813,570:1 +T 8.2% coding (45/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,5701.T 8.2% 340.3 / 27.2 85coding (45/2022 nt)uvrBexision nuclease of nucleotide excision repair, DNA damage recognition component
Reads supporting (aligned to +/- strand):  ref base . (41/37);  new base T (3/4);  total (44/41)
Fisher's exact test for biased strand distribution p-value = 7.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.75e-01

TGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAGCGCCCAACCATGGTAC  >  NC_000913/813424‑813712
                                                                                                                                                   |                                                                                                                                              
tGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCg                                                                                                                                               <  2:692787/148‑1 (MQ=255)
tGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGa                                                                                                                                              >  4:791497/1‑149 (MQ=255)
 ggTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGAt                                                                                                                                             <  4:550641/149‑1 (MQ=255)
 ggTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGAt                                                                                                                                             <  3:411203/149‑1 (MQ=255)
   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                                                                                 >  2:602415/1‑142 (MQ=255)
                                                                    tACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACg                                                                                 <  1:602415/142‑1 (MQ=255)
                                                                    tACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACtt                                                                           <  2:682786/148‑1 (MQ=255)
                                                                    tACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACtt                                                                           >  2:778680/1‑148 (MQ=255)
                                                                     aCCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAGACCTTCTGGCGATCAGCCAGAGGCGATTCGACGtct                                                                                                                      >  3:571693/1‑104 (MQ=255)
                                                                     aCCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAGACCTTCTGGCGATCAGCCAGAGGCGATTCGACGtct                                                                                                                      <  4:571693/104‑1 (MQ=255)
                                                                      ccTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAggg                                                                                                             >  2:598746/1‑112 (MQ=255)
                                                                      ccTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAggg                                                                                                             <  1:598746/112‑1 (MQ=255)
                                                                              cccAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATg                                                                                                  >  1:902367/1‑115 (MQ=255)
                                                                              cccAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATg                                                                                                  <  2:902367/115‑1 (MQ=255)
                                                                                cAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGtc                                                                                                                       <  2:366606/92‑1 (MQ=255)
                                                                                cAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGtc                                                                                                                       >  1:366606/1‑92 (MQ=255)
                                                                                  aCTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTc                                                             <  2:82152/148‑1 (MQ=255)
                                                                                   cTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTg                                                                     >  1:21530/1‑139 (MQ=255)
                                                                                   cTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTg                                                                     <  2:21530/139‑1 (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                                                                                                    <  4:746154/100‑1 (MQ=255)
                                                                                           ggTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGa                                                                                                    >  3:746154/1‑100 (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                                          >  2:178436/1‑148 (MQ=255)
                                                                                                     cATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCAtt                                          >  1:369222/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)
                                                                                                                                                   |                                                                                                                                              
TGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAGCGCCCAACCATGGTAC  >  NC_000913/813424‑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: