Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,332,366 C→T 6.2% A440V (GCG→GTG)  topA → DNA topoisomerase I, omega subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,332,3660CT 6.2% 203.9 / 2.5 64A440V (GCG→GTG) topADNA topoisomerase I, omega subunit
Reads supporting (aligned to +/- strand):  ref base C (25/35);  new base T (2/2);  total (27/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.87e-01

TACCAGTTAATCTGGCGTCAGTTCGTTGCCTGCCAGATGACCCCAGCGAAATATGACTCCACGACGCTGACCGTTGGTGCGGGCGATTTCCGCCTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAGCCGCCAGCCCGTTTCAGTGAAGCATCGCTGGTTAAAGAGCTGGAAAAACGC  >  NC_000913/1332218‑1332514
                                                                                                                                                    |                                                                                                                                                    
tACCAGTTAATCTGGCGTCAGTTCGTTGCCTGCCAGATGACACCAGCGAAATATGACTCCACGACGCTGACCGTTGGTGCGGGCGATTTCCGCCTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCtgc                                                                                                                                                      >  6:35825/1‑149 (MQ=255)
tACCAGTTAATCTGGCGTCAGTTCGTGGCCTGCCAGATGACCCCAGCGAAATATGACTCCACGACGCTGACCGTTGGGGCGGGCGATTTCCGCCTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCGGGACAAAAGTGATGCCtgc                                                                                                                                                      <  6:83769/149‑1 (MQ=255)
     gTTAATCTGGCGTCAGTTCGTTGCCTGCCAGATGACCCCAGCGAAATATGACTCCACGACGCTGACCGTTGGTGCGGGCGATTTCCGCCTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCtgcgttgc                                                                                                                                                 <  3:18044/149‑1 (MQ=255)
        aaTCTGGCGTCAGTTCGTTGCCTGCCAGATGACCCCAGCGAAATATGACTCCACGACGCTGACCGTTGGTGCGGGCGATTTCCGCCTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCtgcgttgcgt                                                                                                                                               >  3:324879/1‑148 (MQ=255)
              gCGTCAGTTCGTTGCCTGCCAGATGACCCCAGCGAAATATGACTCCACGACGCTGACCGTTGGTGCGGGCGATTTCCGCCTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAgg                                                                                                                                          <  6:17609/147‑1 (MQ=255)
                gTCAGTTCGTTGCCTGCCAGATGACCCCAGCGAAATATGACTCCACGACGCTGACCGTTGGTGCGGGCGATTTCCGCCTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGAt                                                                                                                                      <  8:26703/149‑1 (MQ=255)
                    gTTCGTTGCCTGCCAGATGACCCCAGCGAAATATGACTCCACGACGCTGACCGTTGGTGCGGGCGATTTCCGCCTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAg                                                                                                                                  <  3:250715/149‑1 (MQ=255)
                        gTTGCCTGCCAGATGACCCCAGCGAAATATGACTCCACGACGCTGACCGTTGGTGCGGGCGATTTCCGCCTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATc                                                                                                                               <  3:190064/148‑1 (MQ=255)
                           gcctgccAGATGACCCCAGCGAAATATGACTCCACGACGCTGACCGTTGGTGCGGGCGATTTCCGCCTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCa                                                                                                                            <  6:65005/148‑1 (MQ=255)
                            cctgccAGATGACCCCAGCGAAATATGACTCCACGACGCTGACCGTTGGTGCGGGCGATTTCCGCCTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATc                                                                                                                          >  5:98787/1‑149 (MQ=255)
                             ctgccAGATGACCCCAGCGAAATATGACTCCACGACGCTGACCGTTGGTGCGGGCGATTTCCGCCTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCt                                                                                                                         >  2:236093/1‑149 (MQ=255)
                                ccAGATGACCCCAGCGAAATATGACTCCACGCCGCTGACCGTTGGTGCGGGCGATTTCCGCCTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCtt                                                                                                                        <  5:337313/147‑1 (MQ=255)
                                ccAGATGACCCCAGCGAAATATGACTCCACGACGCTGACCGTTGGTGCGGGCGATTTCCGCCTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTAc                                                                                                                      <  1:345569/149‑1 (MQ=255)
                                    aTGACCCCAGCGAAATATGACTCCACGACGCTGACCGTTGGTGCGGGCGATTTCCGCCTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACcag                                                                                                                   >  7:337003/1‑148 (MQ=255)
                                         cccAGCGAAATATGACTCCACGACGCTGACCGTTGGTGCGGGCGATTTCCGCCTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGtt                                                                                                              >  5:308256/1‑148 (MQ=255)
                                                aaaTATGACTCCACGACGCTGACCGTTGGTGCGGGCGATTTCCGCCTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAgg                                                                                                      <  6:98787/149‑1 (MQ=255)
                                                   tatGACTCCACGACGCTGACCGTTGGTGCGGGCGATTTCCGCCTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGa                                                                                                   <  6:174913/149‑1 (MQ=255)
                                                    atGACTCCACGACGCTGACCGTTGGTGCGGGCGATTTCCGCCTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGAt                                                                                                  >  1:180256/1‑149 (MQ=255)
                                                         tCCACGACGCTGACCGTTGGTGCGGGCGATTTCCGCCTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGt                                                                                                               >  3:198111/1‑131 (MQ=255)
                                                         tCCACGACGCTGACCGTTGGTGCGGGCGATTTCCGCCTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGt                                                                                                               <  4:198111/131‑1 (MQ=255)
                                                           cACGACGCTGACCGTTGGTGCGGGCGATTTCCGCCTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGATGctct                                                                                             <  5:149253/147‑1 (MQ=255)
                                                           cACGACGCTGACCGTTGGTGCGGGCGATTTCCGCCTGAAAGCACGCGGTCGTATTTTGAGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGATGCTCTGa                                                                                           >  8:325980/1‑149 (MQ=255)
                                                           cACGACGCTGACCGTTGGTGCGGGCGATTTCCGCCGGAAAGGACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGATGCTCTGa                                                                                           <  1:90198/149‑1 (MQ=255)
                                                            acgacgCTGACCGTTGGTGCGGGCGATTTCCGCCTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGATGCTCTGa                                                                                           <  2:11613/148‑1 (MQ=255)
                                                            acgacgCTGACAGTTGGTGCGGGCGATTTCCGCCTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGATGCTCTGa                                                                                           >  1:11613/1‑148 (MQ=255)
                                                                 gCTGACCGTTGGTGCGGGCGATTTCCGCCTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGAt                                                                                                                                >  6:155963/1‑106 (MQ=255)
                                                                 gCTGACCGTTGGTGCGGGCGATTTCCGCCTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGAt                                                                                                                                <  5:155963/106‑1 (MQ=255)
                                                                   tGACCGTTGGTGCGGGCGATTTCCGCCTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGGTGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGtt                                                                                   >  4:301081/1‑149 (MQ=255)
                                                                    gACCGTTGCTGCGGGCGATTTCCGCCTGAAAGCACGCGGTCGTATTTTCCGTTTTGATGGCTGCACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGATGCTCTGCCGCTCGtt                                                                                   >  6:316222/1‑148 (MQ=255)
                                                                     aCCGTTGGTGCGGGCGATTTCCGCCTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGATGCTCTGACTCTCGTTGa                                                                                 <  5:35825/149‑1 (MQ=255)
                                                                     aCCGTTGGTGCGGGCGATTTCCGCCTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTg                                                                                  >  4:89438/1‑148 (MQ=255)
                                                                        gTTGGTGCGGGCGATTTCCGCCTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACt                                                                              <  8:337003/149‑1 (MQ=255)
                                                                               cGGGCGATTTCCGCCTGAAAGCACGCGGTCGTATTTTGAGTTATGATGGCTGGACAAAAGTGAAGCCTGCGTTGCGTAAAGGCGATGAACATCGCATCTTACCAGCAGTTAATAAAGGCGATGCTCTGACTCTCGTTGAACTTACAcc                                                                        <  5:316222/148‑1 (MQ=255)
                                                                                gggCGATTTCCGCCTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAgg                                                                                                      >  3:30043/1‑117 (MQ=255)
                                                                                gggCGATTTCCGCCTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAgg                                                                                                      <  4:30043/117‑1 (MQ=255)
                                                                                gggCGATTTCCGCCTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTacac                                                                         <  4:149932/146‑1 (MQ=255)
                                                                                gggCGATTTCCGCCTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTacac                                                                         >  3:149932/1‑146 (MQ=255)
                                                                                 ggCGATTTCCGCCTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAg                                                                                                       <  4:5066/115‑1 (MQ=255)
                                                                                  gCGATTTCCGCCTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGtt                                                                                   <  2:259261/134‑1 (MQ=255)
                                                                                  gCGATTTCCGCCTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGtt                                                                                   >  1:259261/1‑134 (MQ=255)
                                                                                  gCGATTTCCGCCTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACAccagcc                                                                    >  2:282107/1‑149 (MQ=255)
                                                                                    gATTTCCGCCTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACAccagccca                                                                  <  2:180256/149‑1 (MQ=255)
                                                                                    gATTTCCGCCTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACAccagcc                                                                    >  6:55632/1‑147 (MQ=255)
                                                                                     aTTTCCGCCTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACAccagcccag                                                                 >  7:271190/1‑149 (MQ=255)
                                                                                        tCCGCCTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCa                                                               <  5:258649/148‑1 (MQ=255)
                                                                                         ccGCCTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCAc                                                              <  5:344076/148‑1 (MQ=255)
                                                                                               gAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTAcc                                                        >  4:324103/1‑148 (MQ=255)
                                                                                               gAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTAcc                                                        >  4:239037/1‑148 (MQ=255)
                                                                                                aaaGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTAcc                                                        >  2:191259/1‑147 (MQ=255)
                                                                                                aaaGCACGCGGTCGTATTTTGAGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCa                                                       <  7:325980/148‑1 (MQ=255)
                                                                                                 aaGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAg                                                     <  1:240130/149‑1 (MQ=255)
                                                                                                   gCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAgcc                                                   <  6:308256/149‑1 (MQ=255)
                                                                                                          gTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAGCCGCCAGcc                                            >  2:202011/1‑149 (MQ=255)
                                                                                                          gTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAGCCGCCAGcc                                            <  3:89438/149‑1 (MQ=255)
                                                                                                           tCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGTGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAGCCGCCAGccc                                           >  3:336454/1‑149 (MQ=255)
                                                                                                           tCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGTGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAGCCGCCAGccc                                           >  3:341709/1‑149 (MQ=255)
                                                                                                           tCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAGCCGCCAGccc                                           >  8:132345/1‑149 (MQ=255)
                                                                                                                  ttGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAGCCGCCAGCCCGTTTCa                                     <  3:324103/148‑1 (MQ=255)
                                                                                                                    gCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAGCCGCCAGCCCGTTTCAGt                                   <  8:271190/148‑1 (MQ=255)
                                                                                                                              ggCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAGCCGCCAGCCCGTTTCAGTGAAGCATCGCt                        >  7:194689/1‑149 (MQ=255)
                                                                                                                               gCTGGTCTAAAGTGAAGCCTGCGTTGCGTAAAGGCGATGAAGATCTCATCTTACCCGCAGTTAATAAAGCCGATGCTCTGACGCTCGTTGAACTTTCACCAGCCCAGCACTTTACCCAGCCGCCAGCCCGTTTCAGTGAAGCATCGCTg                       <  6:97960/149‑1 (MQ=255)
                                                                                                                               gCTGGACAAAAGTGATGCCTGTGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAGCCGCCAGCCCGTTTCAGTGAAGCATCGCTg                       <  4:341709/149‑1 (MQ=255)
                                                                                                                               gCTGGACAAAAGTGATGCCTGTGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAGCCGCCAGCCCGTTTCAGTGAAGCATCGCTg                       <  4:336454/149‑1 (MQ=255)
                                                                                                                                   gACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAGCCGCCAGCCCGTTTCAGTGAAGCATCGCTGGtt                    <  8:194689/148‑1 (MQ=255)
                                                                                                                                   gACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAGCCGCCAGCCCGTTTCAGTGAAGCATCGCTGGt                     <  8:109886/147‑1 (MQ=255)
                                                                                                                                     cAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAGCCGCCAGCCCGTTTCAGTGAAGCATCGCTGGTTaaa                 <  5:55632/149‑1 (MQ=255)
                                                                                                                                     cAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAGCCGCCAGCCCGTTTCAGTGAAGCATCGCTGGTTaaa                 <  7:132345/149‑1 (MQ=255)
                                                                                                                                      aaaaGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAGCCGCCAGCCCGTTTCAGTGAAGCATCGCTGGTTa                   >  4:97972/1‑146 (MQ=255)
                                                                                                                                                   gcgttgcgtAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAGCCGCCAGCCCGTTTCAGTGAAGCATCGCTGGTTAAAGAGCTGGAAAAAcg   >  2:73769/1‑149 (MQ=255)
                                                                                                                                                    cgttgcgtAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAGCCGCCAGCCCGTTTCAGTGAAGCATCGCTGGTTAAAGAGCTGGAAAAAcgc  <  7:27581/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TACCAGTTAATCTGGCGTCAGTTCGTTGCCTGCCAGATGACCCCAGCGAAATATGACTCCACGACGCTGACCGTTGGTGCGGGCGATTTCCGCCTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAGCCGCCAGCCCGTTTCAGTGAAGCATCGCTGGTTAAAGAGCTGGAAAAACGC  >  NC_000913/1332218‑1332514

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