Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 813,575 Δ1 bp 5.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. 5.7% 317.6 / 12.5 87coding (50/2022 nt)uvrBexision nuclease of nucleotide excision repair, DNA damage recognition component
Reads supporting (aligned to +/- strand):  ref base A (42/40);  new base . (2/3);  total (44/43)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.67e-01

ATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAGCGCCCAACCATGGTACTTGCGCCC  >  NC_000913/813444‑813720
                                                                                                                                    |                                                                                                                                                 
aTCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGa                                                                                                                                   >  2:94894/1‑148 (MQ=255)
 tCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGAc                                                                                                                                  >  2:1106067/1‑148 (MQ=255)
        taATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGa                                                                                                                           <  1:1106067/148‑1 (MQ=255)
                  gCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTg                                                                                                                 <  1:94894/148‑1 (MQ=255)
                       aTTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAggg                                                                                                                     <  4:1126168/139‑1 (MQ=255)
                       aTTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAggg                                                                                                                     >  3:1126168/1‑139 (MQ=255)
                                       aaaaTTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAg                                                                                                             <  2:468518/131‑1 (MQ=255)
                                       aaaaTTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAg                                                                                                             >  1:468518/1‑131 (MQ=255)
                                                  ccTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAag                                                                                                                         >  3:873624/1‑108 (MQ=255)
                                                  ccTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAag                                                                                                                         <  4:873624/108‑1 (MQ=255)
                                                   cTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTgg                                                                                                    >  4:1133759/1‑128 (MQ=255)
                                                   cTGCCCGCCCAACACCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTgg                                                                                                    <  3:1133759/128‑1 (MQ=255)
                                                            cAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCTGATCGCCagaag                                                                                                                                           >  3:724810/1‑78 (MQ=255)
                                                            cAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCTGATCGCCagaag                                                                                                                                           <  4:724810/80‑3 (MQ=255)
                                                            cAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGa                                                                                           <  4:249533/128‑1 (MQ=255)
                                                            cAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGa                                                                                           >  3:249533/1‑128 (MQ=255)
                                                             aaCTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGa                                                                            >  2:1059614/1‑142 (MQ=255)
                                                             aaCTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGa                                                                            <  1:1059614/142‑1 (MQ=255)
                                                              aCTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCa                                                                                                <  3:825786/121‑1 (MQ=255)
                                                              aCTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCa                                                                                                >  4:872878/1‑121 (MQ=255)
                                                              aCTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCa                                                                                                <  3:872878/121‑1 (MQ=255)
                                                              aCTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCa                                                                                                >  4:825786/1‑121 (MQ=255)
                                                              aCTCCTTCAGGTAGCGACTCATGAATAAACCGTTCAAACTGAATTCCGCTCTTACACCTTCTGGCGATCAACCAGAGGAGATTCAACGTCTCGAAGAGGGGCTGTAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTc                                                                     >  4:708579/1‑148 (MQ=255)
                                                                 ccTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCt                                                                                                                  <  1:622121/100‑1 (MQ=255)
                                                                 ccTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCt                                                                                                                  >  2:622121/1‑100 (MQ=255)
                                                                   ttCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAggg                                                                                                                     <  3:468927/95‑1 (MQ=255)
                                                                   ttCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAggg                                                                                                                     >  4:468927/1‑95 (MQ=255)
                                                                   ttCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGa                                                                >  2:748129/1‑148 (MQ=255)
                                                                    tCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCTGATCGCCagaagg                                                                                                                                          <  2:852663/73‑4 (MQ=38)
                                                                    tCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGtc                                                                                                                               >  2:499747/1‑84 (MQ=255)
                                                                    tCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGtc                                                                                                                               <  1:499747/84‑1 (MQ=255)
                                                                      aGGTAGCGGCTCATGAGTAAACCGTTCAAACTGAATTCCGATTTTAAACCATCTGGCGATCAGCCAGAGGCGATACGACGTCACGAAGAGGGGCTGGAAGATGGCCTGGCTCACCAGACGTTTCTTGGCGTGCCTGGCTCAGGGaaaa                                                             <  3:708579/148‑1 (MQ=255)
                                                                                 cATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCAtt                                                  >  1:1114031/1‑148 (MQ=255)
                                                                                   tGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGa                                                                                                                                        <  1:759705/60‑1 (MQ=255)
                                                                                   tGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGa                                                                                                                                        >  2:759705/1‑60 (MQ=255)
                                                                                   tGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAATCtt                                                         >  4:355518/1‑139 (MQ=255)
                                                                                   tGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCtt                                                         <  3:355518/139‑1 (MQ=255)
                                                                                      gTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAg                                                                                                                                                  <  4:359746/47‑1 (MQ=255)
                                                                                      gTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAg                                                                                                                                                  >  3:359746/1‑47 (MQ=255)
                                                                                           ccGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGAt                                                                                                           >  4:1020953/1‑81 (MQ=255)
                                                                                           ccGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGAt                                                                                                           <  3:1020953/81‑1 (MQ=255)
                                                                                           ccGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCa                                                                    <  1:482736/120‑1 (MQ=255)
                                                                                           ccGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCa                                                                    >  2:482736/1‑120 (MQ=255)
                                                                                           ccGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTc                                                        >  1:1052408/1‑132 (MQ=255)
                                                                                           ccGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTc                                                        <  2:1052408/132‑1 (MQ=255)
                                                                                            cGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTgg                                                                                 <  4:1133553/106‑1 (MQ=255)
                                                                                            cGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTgg                                                                                 >  3:1133553/1‑106 (MQ=255)
                                                                                            cGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCAt                                       <  2:1114031/148‑1 (MQ=255)
                                                                                              ttCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGaaaa                                                             <  2:953156/124‑1 (MQ=255)
                                                                                              ttCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGaaaa                                                             >  1:953156/1‑124 (MQ=255)
                                                                                                   aCTGAATTCCGCTTTTAAACCTTCTGGCTGATCGCCagaa                                                                                                                                            <  1:1003941/40‑2 (MQ=18)
                                                                                                   aCTGAATTCCGCTTTTAAACCTTCTGGCTGATCGCCagaa                                                                                                                                            >  2:1003941/1‑39 (MQ=18)
                                                                                                   aCTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTgg                                                                                                    >  2:260777/1‑80 (MQ=255)
                                                                                                   aCTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTgg                                                                                                    >  2:702388/1‑80 (MQ=255)
                                                                                                   aCTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTgg                                                                                                    <  1:702388/80‑1 (MQ=255)
                                                                                                   aCTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTgg                                                                                                    <  1:260777/80‑1 (MQ=255)
                                                                                                   aCTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGAc                                >  1:687192/1‑148 (MQ=255)
                                                                                                         ttCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCa                                                       >  2:566141/1‑119 (MQ=255)
                                                                                                         ttCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCa                                                       <  1:566141/119‑1 (MQ=255)
                                                                                                          tCCGCTTTTAACCCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGaata                                                             <  4:783694/112‑3 (MQ=255)
                                                                                                          tCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGaaaa                                                             >  3:783694/1‑112 (MQ=255)
                                                                                                          tCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGaaaa                                                             >  3:590934/1‑112 (MQ=255)
                                                                                                          tCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGaaaa                                                             <  4:590934/112‑1 (MQ=255)
                                                                                                          tCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATAGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTACTGACCTTCAgc                         >  3:923837/1‑148 (MQ=255)
                                                                                                             gCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAgcgc                       <  3:143037/147‑1 (MQ=255)
                                                                                                             gCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAgcgc                       >  4:143037/1‑147 (MQ=255)
                                                                                                              cTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAGCGccc                     >  1:727593/1‑148 (MQ=255)
                                                                                                                   aaaCCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAGCGCCCAACCa                <  2:687192/148‑1 (MQ=255)
                                                                                                                    aaCCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGaaa                                                              >  4:241515/1‑101 (MQ=255)
                                                                                                                    aaCCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGaaa                                                              <  3:241515/101‑1 (MQ=255)
                                                                                                                     aCCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCtt                                                         >  2:793678/1‑105 (MQ=255)
                                                                                                                     aCCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCtt                                                         <  1:793678/105‑1 (MQ=255)
                                                                                                                     aCCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAGCGCCCAACCATg              <  2:727593/148‑1 (MQ=255)
                                                                                                                      ccTTCTGGCGATCAGCCAGAGGCGGTTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTAGCTCTGGGAAAACCTTCAGCATTGCCAATGTCATTGCTGACCTTCAGCGCCCAACCATgg             >  3:330184/1‑148 (MQ=255)
                                                                                                                       cTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAGCGCCCAACCATGGt            <  4:923837/148‑1 (MQ=255)
                                                                                                                         tCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTa                                                                                      >  2:1007249/1‑72 (MQ=255)
                                                                                                                         tCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTa                                                                                      <  1:1007249/72‑1 (MQ=255)
                                                                                                                         tCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTg                                                 >  3:483476/1‑109 (MQ=255)
                                                                                                                         tCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTg                                                 <  4:483476/109‑1 (MQ=255)
                                                                                                                         tCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCAtt                                      >  4:136274/1‑120 (MQ=255)
                                                                                                                         tCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCAtt                                      <  3:136274/120‑1 (MQ=255)
                                                                                                                              cGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCAcc                                                     >  4:177667/1‑100 (MQ=255)
                                                                                                                              cGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCAcc                                                     <  3:177667/100‑1 (MQ=255)
                                                                                                                              cGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCAcc                                                     >  2:562074/1‑100 (MQ=255)
                                                                                                                              cGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCAcc                                                     <  1:562074/100‑1 (MQ=255)
                                                                                                                                 aTCAGCCAGAGGCGATTCAGCGTCTCGAAGTGGCGCTAGAAGATGGACTGGGGCACGAGACGTTACTTGGCGTGACTGGCTCAGGGAGAACCTTCCCCACTGCCAATGTCATTGCTGACAGACCGCGCCCAACCATGGTACTTGCGcc   <  4:330184/148‑1 (MQ=255)
                                                                                                                                  tCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAGCGCCCAACCATGGTACTTGCGccc  >  2:671919/1‑148 (MQ=255)
                                                                                                                                    aGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTc                                         >  4:478960/1‑107 (MQ=255)
                                                                                                                                    aGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTc                                         <  3:478960/107‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                                 
ATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATCAGCCAGAGGCGATTCGACGTCTCGAAGAGGGGCTGGAAGATGGCCTGGCGCACCAGACGTTACTTGGCGTGACTGGCTCAGGGAAAACCTTCACCATTGCCAATGTCATTGCTGACCTTCAGCGCCCAACCATGGTACTTGCGCCC  >  NC_000913/813444‑813720

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