Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 814,955 A→T 6.2% H477L (CAC→CTC)  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_000913814,9550AT 6.2% 195.7 / 3.7 65H477L (CAC→CTC) uvrBexision nuclease of nucleotide excision repair, DNA damage recognition component
Reads supporting (aligned to +/- strand):  ref base A (31/30);  new base T (2/2);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

TGGCGACACAGGTTGATGATCTTCTTTCGGAGATTCGTCAGCGAGCGGCAATTAACGAACGCGTACTGGTCACCACACTGACCAAGCGGATGGCGGAAGATCTTACCGAATATCTCGAAGAACATGGCGAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGTCGAACGTATGGAGATTATCCGCGACTTGCGTCTGGGTGAGTTCGACGTGCTGGTAGGGATCAACTTACTGCGCGAAGGTCTGGATATGCCGGAAGTGTCGCTGGTGGCGATCCTCGACGCTGACAAAGAA  >  NC_000913/814808‑815103
                                                                                                                                                   |                                                                                                                                                    
tGGCGACACAGGTTGATGATCTTCTTTCGGAGATTCGTCAGCGAGCGGCAATTAACGAACGCGTACTGGTCACCACACTGACCAAGCGGATGGCGGAAGATCTTACCGAATATCTCGAAGAACATGGCGAGCGCGTGCGTTATCTTCAc                                                                                                                                                     <  8:61630/149‑1 (MQ=255)
  gCGACACAGGTTGATGATCTTCTTTCGGAGATTCGTCAGCGAGCGGCAATTAACGAACGCGTACTGGTCACCACACTGACCAAGCGGATGGCGGAAGATCTTACCGAATATCTCGAAGAACATGGCGAGCGCGTGCGTTATCTTCACt                                                                                                                                                    >  6:16629/1‑148 (MQ=255)
       acaGGTTGATGATCTTCTTTCGGAGATTCGTCAGCGAGCGGCAATTAACGAACGCGTACTGGTCACCACACTGACCAAGCGGATGGCGGAAGATCTTACCGAATATCTCGAAGAACATGGCGAGCGCGTGCGTTATCTTCACTCAGata                                                                                                                                              >  6:137033/1‑149 (MQ=255)
         aGGTTGATGATCTTCTTTCGGAGATTCGTCAGCGAGCGGCAATTAACGAACGCGTACTGGGCACCACACTGACCAAGCGGATGGCGGAAGATCTTACCGAATATCTCGAAGAACATGGCGAGCGCGTGCGTTATCTTCACTCAGATATc                                                                                                                                            <  3:53895/149‑1 (MQ=255)
               atgatCTTCTTTCGGAGATTCGTCAGCGAGCGGCAATTAACGAACGCGTACTGGTCACCACACTGACCAAGCGGATGGCGGAAGATCTTACCGAATATCTCGAAGAACATGGCGAGCGCGTGCGTTATCTTCACTCAGATATCGacac                                                                                                                                       <  5:16629/148‑1 (MQ=255)
                 gatCTTCTTTCGGAGATTCGTCAGCGAGCGGCAATTAACGAACGCGTACTGGTCACCACACTGACCAAGCGGATGGCGGAAGATCTTACCGAATATCTCGAAGAACATGGCGAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGt                                                                                                                                    >  4:171428/1‑149 (MQ=255)
                  atCTTCTTTCGGAGATTCGTCAGCGAGCGGCAATTAACGAACGCGTACTGGTCACCACACTGACCAAGCGGATGGCGGAAGATCTTACCGAATATCTCGAAGAACATGGCGAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGTc                                                                                                                                   >  7:146032/1‑149 (MQ=255)
                     ttcttTCGGAGATTAGTCAGCGAGCGGCAATTAACGAACGCGTACTGGTAACCACACTGACCAAGCGGATGGCGGAAGATCTTACCGAATATCTCGAAGAACATGGCGAGCGCGTGCGTTATCTTCACTCAGATAACGACACCGTCGaa                                                                                                                                >  3:119268/1‑149 (MQ=255)
                      tcttTCGGAGATTCGTCAGCGAGCGGCAATTAACGAACGCGTACTGGTCACCACACTGACCAAGCGGATGGCGGAAGATCTTACCGAATATCTCGAAGAACATGGCGAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGTCGAAc                                                                                                                               >  3:208917/1‑149 (MQ=255)
                           cGGAGATTCGTCAGCGAGCGGCAATTAACGAACGCGTACTGGTCACCACACTGACCAAGCGGATGGCGGAAGATCTTACCGAATATCTCGAAGAACATGGCGAGCGCGGGCGTTATCTTCACTCAGATATCGACACCGTCGAACGTAt                                                                                                                           <  1:281242/148‑1 (MQ=255)
                            ggAGATTCGTCAGCGAGCGGCAATTAACGAACGCGTACTGGTCACCACACTGACCAAGCGGATGGCGGAAGATCTTACCGAATATCTCGAAGAACATGGCGAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGTCGAACGTATg                                                                                                                          >  2:259124/1‑148 (MQ=255)
                             gagaTTCGTCAGCGAGCGGCAATTAACGAACGCGTACTGGTCACCACACTGACCAAGCGGATGGCGGAAGATCTTACCGAATATCTCGAAGAACATGGCGAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGTCGAACGTATGga                                                                                                                        <  7:372857/149‑1 (MQ=255)
                             gagaTTCGGCAGCGAGCTGCAATAAACGAACTCGTACTGGTAACCACACTGACCAAGCGGATGGCGGAAGATCTTACCGAATATCTCGAAGAACATGGAGAGCGCTTGCGTTATCTTCACTCAGATATCGACACCGTCGATCGTATGga                                                                                                                        <  4:251823/149‑1 (MQ=255)
                                   cGTCAGCGAGCGGCAATTACCGAACGCGTACTGGTCACCACACTGACCAAGCGGATGGCGGAAGATCTTACCGAATATCTCGAAGAACATGGCGAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGTCGAACGTATGGAGATTa                                                                                                                   <  4:276566/148‑1 (MQ=255)
                                                  aTTAACGAACGCGTACTGGTCACCACACTGAACAAGCGGATGGCGGAAGATCTTACCGAATATCTCGAAGACAATGGCGAGCGCGTGCGTTATCGTAAATCAGATATCGACACCGTCGAACGTATGGAGATTATCCGCGACTGGAGTCt                                                                                                   >  2:119120/1‑149 (MQ=255)
                                                   ttAACGAACGCGTACTGGTCACCACACTGACCAAGCGGATGGCGGAAGATCTTACCGAATATCTCGAAGAACATGGCGAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGTCGAACGTATGGAGATTATCCGCGACTTGCGTCTg                                                                                                  >  7:2797/1‑149 (MQ=255)
                                                     aacgaacgCGTACTGGTCACCACACTGACCAAGCGGATGGCGGAAGATCTTACCGAATATCTCGAAGAACATGGCGAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGTCGAACGTATGGAGATTATCCGCGACTTGCGTCTggg                                                                                                >  4:193311/1‑149 (MQ=255)
                                                     aacgaacgCGTACTGGTCACCACACTGACCAAGCGGATGGCGGAAGATCTTACCGAATATCTCGAAGAACATGGCGAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGTCGAACGTATGGAGATTATCCGCGACTTGCGTCTgg                                                                                                 <  4:148760/148‑1 (MQ=255)
                                                      acgaacgCGTACTGGTCACCACACTGACCAAGCGGATGGCGGAAGATCTTACCGAATATCTCGAAGAACATGGCGAGCGCGTGCGTTATCTTCTCTCAGATATCGACACCGTCGAACGTATGGAGAt                                                                                                                     >  4:128358/1‑127 (MQ=255)
                                                      acgaacgCGTACTGGTCACCACACTGACCAAGCGGATGGCGGAAGATCTTACCGAATATCTCGAAGAACATGGCGAGCGCGTGCGTTATCTTCTCTCAGATATCGACACCGTCGAACGTATGGAGAt                                                                                                                     >  4:123022/1‑127 (MQ=255)
                                                      acgaacgCGTACTGGTCACCACACTGACCAAGCGGATGGCGGAAGATCTTACCGAATATCTCGAAGAACATGGCGAGCGCGTGCGTTATCTTCTCTCAGATATCGACACCGTCGAACGTATGGAGAt                                                                                                                     <  3:128358/127‑1 (MQ=255)
                                                      acgaacgCGTACTGGTCACCACACTGACCAAGCGGATGGCGGAAGATCTTACCGAATATCTCGAAGAACATGGCGAGCGCGTGCGTTATCTTCTCTCAGATATCGACACCGTCGAACGTATGGAGAt                                                                                                                     <  3:123022/127‑1 (MQ=255)
                                                       cgaacgCGTACTGGTCACCACACTGACCAAGCGGATGGCGGAAGATCTTACCGAATATCTCGAAGAACATGGCGAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGTCGAACGTATGGAGATTATCCGCGACTTGCGTCTGGGt                                                                                               >  2:340642/1‑148 (MQ=255)
                                                                  tGGTCACCACACTGACCAAGCGGATGGCGGAAGATCTTACCGAATATCTCGAAGAACATGGCGAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGTCGAACGTATGGAGATTATCCGCGACTTGCGTCTGGGTGAGTTCGACGTg                                                                                   >  7:38502/1‑149 (MQ=255)
                                                                    gTCACCACACTGACCAAGCGGATGGCGGAAGATCTTACCGAATATCTCGAAGAACATGGCGAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGTCGAACGTATGGAGATTATCCGCGACTTGCGTCTGGGTGAGTTCGACGTGc                                                                                  <  4:217159/148‑1 (MQ=255)
                                                                      caccacACTGACCAAGCGGATGGCGGAAGATCTTACCGAATATCTCGAAGAACATGGCGAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGTCGAACGTATGGAGATTATCCGCGACTTGCGTCTGGGTGAGTTCGACGTGCTgg                                                                               >  3:263488/1‑149 (MQ=255)
                                                                       accacACTGACCAAGCGGATGGCGGAAGATCTTACCGAATATCTCGAAGAACATGGCGAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGTCGAACGTATGGAGATTATCCGCGACTTGCGTCTGGGTGAGTTCGACGTGCTgg                                                                               >  3:43799/1‑148 (MQ=255)
                                                                       accacACTGACCAAGCGGATGGCGGAAGATCTTACCGAATATCTCGAAGAACATGGCGAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGTCGAACGTATGGAGATTATCCGCGACTTGCGTCTGGGTGAGTTCGACGTGCTgg                                                                               <  8:146032/148‑1 (MQ=255)
                                                                         cacaCTGACCAAGCGGATGGCGGAAGATCTTACCGAATATCTCGAAGAACATGGCGAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGTCGAACGTATGGAGATTATCCGCGACTTGCGTCTGGGTGAGTTCGACGTGCTGGTAg                                                                            <  6:182048/149‑1 (MQ=255)
                                                                           cacTGACCAAGCGGATGGCGGAAGATCTTACCGAATATCTCGAAGAACATGGCGAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGTCGAACGTATGGAGATTATCCGCGACTTGCGTCTGGGTGAGTTCGACGTGCTGGTAgg                                                                           >  6:385341/1‑148 (MQ=255)
                                                                           cacGGACCAAGCGGATGGCGGAAGATCTTACCGAATATCTCGAAGAAGATGGCGAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGTCGAACGTATGGAGATTATCCGCGACTTGCGGCTGGGTGAGTTCGACGTGCTGGTAggg                                                                          <  3:40320/149‑1 (MQ=255)
                                                                            acTGACCAAGGGGATGTCGGATGATCTTAACGAATATGTCGAAGAACATGGCGAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGTCGAACGTATGGAGATTATCCGCGACTTGCGTCTGGGTGAGTTCGACGTGCTGGTAGGGa                                                                         <  4:119268/149‑1 (MQ=255)
                                                                                aCCAAGCGGATGGCGGAAGATCTTACCGAATATCTCGAAGAACATGGCGAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGTCGAACGTATGGAGATTATCCGCGACTTGCGTCTGGGTGAGTTCGACGTGCTGGTAGGGATCaa                                                                     >  3:82101/1‑149 (MQ=255)
                                                                                aCCAAGCGGATGGCGGAAGATCTTACCGAATATCTCGAAGAACATGGCGAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGTCGAACGTATGGAGATTATCCGCGACTTGCGTCTGGGTGAGTTCGACGTGCTGGTAGGGATCa                                                                      <  2:358887/148‑1 (MQ=255)
                                                                                 ccAAGCGGATGGCGGAAGATCTTACCGAATATCTCGAAGAACATGGCGAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGTCGAACGTATGGAGATTATCCGCGACTTGCGTCTGGGTGAGTTCGACGTGCTGGTAGGGATCAAc                                                                    <  5:385341/149‑1 (MQ=255)
                                                                                  cAAGCGGATGGCGGAAGATCTTACCGAATATCTCGAAGAACATGGCTAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGTCGAACGTATGGAGATTATCCGCGACTTGCGTCTGGGTGAGTTCGACGTGCTGGTAGGGATCAACt                                                                   <  4:43799/149‑1 (MQ=255)
                                                                                    aGCGGATGGCGGAAGATCTTACCGAATATCTCGAAGAACATGGCGAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGTCGAACGTATGGAGATTATCCGCGACTTGCGTCTGGGTGAGTTc                                                                                         >  5:253778/1‑125 (MQ=255)
                                                                                    aGCGGATGGCGGAAGATCTTACCGAATATCTCGAAGAACATGGCGAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGTCGAACGTATGGAGATTATCCGCGACTTGCGTCTGGGTGAGTTc                                                                                         <  6:253778/125‑1 (MQ=255)
                                                                                    aGCGGATGGCGGAAGATCTTACCGAATATCTCGAAGAACATGGCGAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGTCGAACGTATGGAGATTATCCGCGACTTGCGTCTGGGTGAGTTCGACGTGCTGGTAGGGATCAACTTa                                                                 <  4:263488/149‑1 (MQ=255)
                                                                                             cGGAAGATCTTACCGAATATCTCGAAGAACATGGCGAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGTCGAACGTATGGAGATTATCCGCGACTTGCGTCTGGGTGAGTTCGACGTGCTGGt                                                                              <  5:242410/127‑1 (MQ=255)
                                                                                             cGGAAGATCTTACCGAATATCTCGAAGAACATGGCGAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGTCGAACGTATGGAGATTATCCGCGACTTGCGTCTGGGTGAGTTCGACGTGCTGGt                                                                              >  6:242410/1‑127 (MQ=255)
                                                                                                       tACCGAATATCTCGAAGAACATGGCGAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGTCGAACGTATGGAGATTATCCGCGACTTGCGTCTGGGTGAGTTCGACGTGCTGGTAGGGATCAACTTACTGCGCGAAGGTCTGGata                                              >  7:161284/1‑149 (MQ=255)
                                                                                                        aCCGAATATCTCGAAGAACATGGCGAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGTCGAACGTATGGAGATTATCCGCGACTTGCGTCTGGGTGAGTTCGACGTGCTGGTAGGGATCAACTTACTGCGCGAAGGTCTg                                                  <  1:259124/144‑1 (MQ=255)
                                                                                                        aCCGAATATCTCGAAGAACATGGCGAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGTCGAACGTATGGAGATTATCCGCGACTTGCGTCTGGGTGAGTTCGACGTGCTGGTAGGGATCAACTTACTGCGCGAAGGTCTGGa                                                <  8:38502/146‑1 (MQ=255)
                                                                                                          cGAATATCTCGAAGAACATGGCGAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGTCGAACGTATGGAGATTATCCGCGACTTGCGTCTGGGTGAGTTCGACGTGCTGGTAGGGATCAACTTACTGCGCGaa                                                        <  7:203928/136‑1 (MQ=255)
                                                                                                          cGAATATCTCGAAGAACATGGCGAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGTCGAACGTATGGAGATTATCCGCGACTTGCGTCTGGGTGAGTTCGACGTGCTGGTAGGGATCAACTTACTGCGCGaa                                                        >  8:203928/1‑136 (MQ=255)
                                                                                                           gAATATCTCGAAGAACATGGCGAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGTCGAACGTATGGAGATTATCCGCGACTTGCGTCTGGGTGAGTTCGACGTGCTGGTAGGGATCAACTTACTGCGCGAAGGTCt                                                   <  7:331341/140‑1 (MQ=255)
                                                                                                           gAATATCTCGAAGAACATGGCGAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGTCGAACGTATGGAGATTATCCGCGACTTGCGTCTGGGTGAGTTCGACGTGCTGGTAGGGATCAACTTACTGCGCGAAGGTCt                                                   >  8:331341/1‑140 (MQ=255)
                                                                                                            aaTATCTCGAAGAACATGGCGAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGTCGAACGTATGGAGATTATCCGCGACTTGCGTCTGGGTGAGTTCGACGTGCTGGTAGGGATCAACTTACTGCGCGAAGGTCTGGATATGCCg                                         >  1:204655/1‑149 (MQ=255)
                                                                                                             atatCTCGAAGAACATGGCGAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGTCGAACGTATGGAGATTATCCGCGACTTGCGTCTGGGTGAGTTCGACGTGCTGGTAGGGATCAACTTACTGCGCGAAGGTCTGGATATGCCgg                                        <  4:208917/149‑1 (MQ=255)
                                                                                                              tatGTCGAAGAACATGGCGAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGTCGAACGTATGGAGTTTATCCGCGACTTGCGTCTGGGTGAGTTCGACGTGCTGGTAGGGATCAACTTACTGCGCGAAGGTCTGGATATTCCGGa                                       <  5:9365/149‑1 (MQ=255)
                                                                                                                tctcGAAGAACATGGCGAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGTCGAACGTATGGAGATTATCCGCGACTTGCGTCTGGGTGAGTTCGACGTGCTGGTAGGGATCAACTTACTGCGCGAAGGTCTGGATATGCCGGAAg                                     >  8:377767/1‑149 (MQ=255)
                                                                                                                tctcGAAGAACATGGCGAGCGCGGGCGTTATCTTCACTCAGATATCGACACCGTCGAACGTATGGAGATTATCCGCGACTTGCGTCTGGGTGAGTTCGACGTGCTGGTAGGGATCAACTTACTGCGCGAAGGTCTGGATATGCCGGaa                                      <  5:137033/148‑1 (MQ=255)
                                                                                                                 ctcGAAGAACATGGCGAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGTCGAACGTATGGAGATTATCCGCGACTTGCGTCTGGGTGAGTTCGACGTGCTGGTAGGGATCAACTTACTGCGCGAAGGTCTGGATATGCCGGAAgt                                    >  4:66689/1‑149 (MQ=255)
                                                                                                                  tcGAAGAACATGGCGAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGTCGAACGTATGGAGATTATCCGCGACTTGCGTCTGGGTGAGTTCGACGTGCTGGTAGGGATCAACTTACTGCGCGAAGGTCTGGATATGCCGGAAgtg                                   >  8:181204/1‑149 (MQ=255)
                                                                                                                     aagaaCATGGCGAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGTCGAACGTATGGAGATTATCCGCGACTTGCGTCTGGGTGAGTTCGACGTGCTGGTAGGGATCAACTTACTGCGCGAAGGTCTGGATATGCCGGAAGTGTCg                                >  5:70322/1‑149 (MQ=255)
                                                                                                                      agaaCATGGCGAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGTCGAACGTATGGAGATTATCCGCGACTTGCGTCTGGGTGAGTTCGACGTGCTGGTAGGGATCAACTTACTGCGCGAAGGTCTGGATATGCCGGAAGTGTCGc                               <  7:181204/149‑1 (MQ=255)
                                                                                                                         aCATGGCGAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGTCGAACGTATGGAGATTATCCGCGACTTGCGTCTGGGTGAGTTCGACGTGCTGGTAGGGATCAACTTACTGCGCGAAGGTCTGGATATGCCGGAAGTGTCGCtgg                            >  5:317031/1‑149 (MQ=255)
                                                                                                                         aCATGGCGAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGTCGAACGTATGGAGATTATCCGCGACTTGCGTCTGGGTGAGTTCGACGTGCTGGTAGGGATCAACTTACTGCGCGAAGGTCTGGATATGCCGGAAGTGTCGCtgg                            <  3:193311/149‑1 (MQ=255)
                                                                                                                              gCGAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGTCGAACGTATGGAGATTATCCGCGACTTGCGTCTGGGTGAGTTCGACGTGCTGGTAGGGATCAACTTACTGCGCGAAGGTCTGGat                                               >  8:63699/1‑125 (MQ=255)
                                                                                                                              gCGAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGTCGAACGTATGGAGATTATCCGCGACTTGCGTCTGGGTGAGTTCGACGTGCTGGTAGGGATCAACTTACTGCGCGAAGGTCTGGat                                               <  7:63699/125‑1 (MQ=255)
                                                                                                                              gCGAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGTCGAACGTATGGAGATTATCCGCGACTTGCGTCTGGGTGAGTTCGACGTGCTGGTAGGGATCAACTTACTGCGCGAAGGTCTGGATATGCCGGAAGTGTCGCTGGTGGc                        >  2:115856/1‑148 (MQ=255)
                                                                                                                                                ttCACTCAGATATCGACACCGTCGAACGTATGGAGATTATCCGCGACTTGCGTCTGGGTGAGTTCGACGTGCTGGTAGGGATCAACTTACTGCGCGAAGGTCTGGATATGCCGGAAGTGTCGCTGGTGGCGATCCTCGACGCTGACa       >  5:126402/1‑147 (MQ=255)
                                                                                                                                                ttCACTCAGATATCGACACCGTCGAACGTATGGAGATTATCCGCGACTTGCGTCTGGGTGAGTTCGACGTGCTGGTAGGGATCAACTTACTGCGCGAAGGTCTGGATATGCCGGAAGTGTCGCTGGTGGCGATCCTCGACGCTGACa       >  4:283079/1‑147 (MQ=255)
                                                                                                                                                 tCACTCAGATATCGACACCGTCGAACGTATGGAGATTATCCGCGACTTGCGTCTGGGTGAGTTCGACGTGCTGGTAGGGATCAACTTACTGCGCGAAGGTCTGGATATGCCGGAAGTGTCGCtggtgg                         <  4:182216/128‑1 (MQ=255)
                                                                                                                                                   aCTCAGATATCGACACCGTCGAACGTATGGAGATTATCCGCGACTTGCGTCTGGGTGAGTTCGACGTGCTGGTAGGGATCAACTTACTGCGCGAAGGTCTGGATATGCCGGAAGTGTCACTGGTGGCGATCCTCGACGCTGACAaagaa  <  8:188615/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TGGCGACACAGGTTGATGATCTTCTTTCGGAGATTCGTCAGCGAGCGGCAATTAACGAACGCGTACTGGTCACCACACTGACCAAGCGGATGGCGGAAGATCTTACCGAATATCTCGAAGAACATGGCGAGCGCGTGCGTTATCTTCACTCAGATATCGACACCGTCGAACGTATGGAGATTATCCGCGACTTGCGTCTGGGTGAGTTCGACGTGCTGGTAGGGATCAACTTACTGCGCGAAGGTCTGGATATGCCGGAAGTGTCGCTGGTGGCGATCCTCGACGCTGACAAAGAA  >  NC_000913/814808‑815103

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑