Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,442,668 G→T 9.0% P602P (CCG→CCT topA → DNA topoisomerase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,442,6680GT 9.0% 219.5 / 8.3 74P602P (CCG→CCTtopADNA topoisomerase I
Reads supporting (aligned to +/- strand):  ref base G (29/37);  new base T (3/4);  total (33/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGAGAACCTCGACGACGTCGCCCAGGGCGAGCGTGACTGGAAGAACGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTGGCTACAGCCTGCCGCCGAAAGAGCGCTGCAAGGCTACCGTCAACCTGGTGCCGGGCGACGAGATTGCCGC  >  NC_002947/2442520‑2442805
                                                                                                                                                    |                                                                                                                                          
aaGAGAACCTCGACGACGTCGCCCAGGGCGAGCGTGACTGGAAGAACGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCg                                                                                                                                            >  6:121497/1‑149 (MQ=255)
     aaCCTCGACGACGTCGCCCAGGGCGAGCGTGACTGGAAGAACGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGc                                                                                                                                         >  4:133884/1‑147 (MQ=255)
     aaCCTCGACGACGTCGCCCAGGGCGAGCGTGACTGGAAGAACGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCaa                                                                                                                                       <  8:111934/149‑1 (MQ=255)
        cTCGACGACGTCGCCCAGGGCGAGCGTGACTGGAAGAACGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAgg                                                                                                                                     >  4:298403/1‑148 (MQ=255)
          cgacaaCGTCGCCCAGGGCGAGCGTGACTGGAAGAACGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAAt                                                                                                                                  >  8:305607/1‑149 (MQ=255)
             cgacgTCGCCCAGGGCGAGCGTGACTGGAAGAACGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCACATTCCTTGCACGGAAtgttg                                                                                                                             <  7:14521/149‑3 (MQ=255)
              gacgTCGCCCAGGGCGAGCGTGACTGGAAGAACGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCTTGCACGGAAtgttg                                                                                                                             <  7:115327/149‑3 (MQ=255)
              gacgTCGCCCAGGGCGAGCGTGACTGGAAGAACGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTgg                                                                                                                             <  8:299247/149‑1 (MQ=255)
                cgTCGCCCAGGGCGAGCGTGACTGGAAGAACGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCTTGCACGGAAtgttggt                                                                                                                           <  2:176463/149‑5 (MQ=255)
                  tCGCCCAGGGCGAGCGTGACTGGAAGAACGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCgg                                                                                                                         <  4:792/149‑1 (MQ=255)
                      ccAGGGCGAGCGTGACTGGAAGAACGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAAtgtg                                                                                                                              >  2:301002/1‑140 (MQ=255)
                      ccAGGGCGAGCGTGACTGGAAGAACGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAAtgtg                                                                                                                              <  1:301002/140‑1 (MQ=255)
                      ccAGGGCGAGCGTGACTGGAAGAACGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCga                                                                                                                     >  8:268119/1‑149 (MQ=255)
                         gggCGAGCGTGACTGGAAGAACGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCgatga                                                                                                                  <  1:374263/149‑1 (MQ=255)
                           gCGAGCGTGACTGGAAGAACGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCgatgatg                                                                                                                <  2:157374/149‑1 (MQ=255)
                            cGAGCGTGACTGGAAGAACGTGCGCGACGAGTTCTACGGCGTATTCAGCAAAATGCTGCAGTCTGGCGAGTCCAGCGAAGACGGCATGCTTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCgatgatga                                                                                                               <  6:112094/149‑1 (MQ=255)
                              aGCGTGACTGGAAGAACGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCTTGCAAGGAATGTGGTCGGCCGATGATGAtt                                                                                                             <  1:225541/149‑1 (MQ=255)
                              aGCGTGACTGGAAGAACGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCGATGATGAtt                                                                                                             <  3:298403/149‑1 (MQ=255)
                               gCGTGACTGGAAGAACGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCGATGATGATTc                                                                                                            <  3:133884/149‑1 (MQ=255)
                                  tGACTGGAAGAACGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCGATGATGAtt                                                                                                             <  5:58100/145‑1 (MQ=255)
                                  tGACTGGAAGAACGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCGATGATGAtt                                                                                                             >  6:58100/1‑145 (MQ=255)
                                    aCTGGAAGAACGTGCTCGACGAGTTCTAGGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCGATGATGATTCGCAc                                                                                                        <  8:353758/148‑1 (MQ=255)
                                    aCTGGAAGAACGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGGCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACGAACATTCCGTGCAAGGAATGTGGTCGGCCGATGATGATTCGCAcc                                                                                                       <  5:121497/149‑1 (MQ=255)
                                       ggAAGAACGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCgg                                                                                                                         <  8:296278/128‑1 (MQ=255)
                                       ggAAGAACGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCgg                                                                                                                         >  7:296278/1‑128 (MQ=255)
                                        gaagaaCGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCGATGATGATTCGCACCGCt                                                                                                    >  1:357071/1‑148 (MQ=255)
                                        gaagaaCGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCAGATGATGCTTCGCACCGCtt                                                                                                   >  2:186499/1‑149 (MQ=255)
                                           gaaCGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCGATGATGATTCGCACCGCTTCTa                                                                                                >  3:57793/1‑149 (MQ=255)
                                              cGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCGATGATGATTCGCACCGCTTCTACgg                                                                                             <  4:57793/149‑1 (MQ=255)
                                              cGTGCTCGACGAGTTATACCGCGATTTCAGCAAAAACCTGCAGACTACCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCCTGCAAGGAATGTGGTCGCCCGATGATGATTCGCACAGCTTCTACgg                                                                                             >  8:114918/1‑149 (MQ=255)
                                                 gCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCgat                                                                                                                    <  3:282120/123‑1 (MQ=255)
                                                 gCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCgat                                                                                                                    >  4:282120/1‑123 (MQ=255)
                                                 gCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCGATGATGATTCGCACCGCTTCTACGGgtg                                                                                          >  4:87784/1‑149 (MQ=255)
                                                     gacgaGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCGATGATGATTCGCACCGCTTCTACGGgt                                                                                           >  2:353547/1‑144 (MQ=255)
                                                     gacgaGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCGATGATGATTCGCACCGCTTCTACGGgt                                                                                           <  1:353547/144‑1 (MQ=255)
                                                      acgaGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTc                                                                                     >  5:266742/1‑149 (MQ=255)
                                                        gaGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCGATGATGATTCGCACCGc                                                                                                     <  7:76441/131‑1 (MQ=255)
                                                        gaGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCGATGATGATTCGCACCGc                                                                                                     >  8:76441/1‑131 (MQ=255)
                                                              tACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCt                                                                                   <  5:183392/143‑1 (MQ=255)
                                                              tACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCt                                                                                   >  6:183392/1‑143 (MQ=255)
                                                                  gCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGctct                                                                         <  6:86324/149‑1 (MQ=255)
                                                                    gATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTgg                                                                       <  3:87784/149‑1 (MQ=255)
                                                                            cAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTgg                                                                       >  4:38275/1‑141 (MQ=255)
                                                                            cAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTgg                                                                       <  3:38275/141‑1 (MQ=255)
                                                                                      cAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGctct                                                                         >  6:85405/1‑129 (MQ=255)
                                                                                      cAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGctct                                                                         <  5:85405/129‑1 (MQ=255)
                                                                                             ccGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCGATGATGATTCGCACCGCTTCTACGGgt                                                                                           <  4:127246/104‑1 (MQ=255)
                                                                                             ccGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCGATGATGATTCGCACCGCTTCTACGGgt                                                                                           >  3:127246/1‑104 (MQ=255)
                                                                                             ccGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTg                                                                             >  6:159992/1‑118 (MQ=255)
                                                                                             ccGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTg                                                                             <  5:159992/118‑1 (MQ=255)
                                                                                                  tCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTGGCTACAGCCTGCCGCCGAAAGAGCGCTGCaa                                         <  2:357071/149‑1 (MQ=255)
                                                                                                    cAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTGGCTACAGCCTGCCGCCGAAAGAGCGCTGCAAgg                                       <  6:266742/149‑1 (MQ=255)
                                                                                                    cAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTGGCTACAGCCTGCCGCCGAAAGAGCGCTGCAAgg                                       >  2:359118/1‑149 (MQ=255)
                                                                                                      gCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTGGCTACAGCCTGCCGCCGAAAGAGCGCTGCAAGGc                                      >  8:273129/1‑148 (MQ=255)
                                                                                                          aGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTGGCTACAGCCTGCCGCCGAAAGAGCGCTGCAAGGCTACCg                                 <  7:273129/149‑1 (MQ=255)
                                                                                                            aCGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTGGCTACAGCCTGCCGCCGAAAGAGCGCTGCAAGGCt                                     <  2:80588/143‑1 (MQ=255)
                                                                                                            aCGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTGGCTACAGCCTGCCGCCGAAAGAGCGCTGCAAGGCt                                     >  1:80588/1‑143 (MQ=255)
                                                                                                                cATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCgatgat                                                                                                                 >  5:76019/1‑63 (MQ=255)
                                                                                                                cATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCgatgat                                                                                                                 <  6:76019/63‑1 (MQ=255)
                                                                                                                        ccAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTGGCTACAGCCTGCCGCCGAAAGAGCGCTGCAAGGCTACCGTCAACCTGGTGcc                    <  5:333909/148‑1 (MQ=255)
                                                                                                                        ccAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTGGCTACAGCCTGCCGCCGAAAGAGCGCGGCAAGGCTACCGTCAACCTGGTGCCg                   >  3:54905/1‑149 (MQ=255)
                                                                                                                          aaCCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGctct                                                                         >  7:89144/1‑93 (MQ=255)
                                                                                                                          aaCCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGctct                                                                         <  8:89144/93‑1 (MQ=255)
                                                                                                                           aCCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTGGCTACAGCCTGCCGCCGAAAGAGCGCTGCAAGGCTACCGTCAACCTGGTGCCGGGc                >  1:186301/1‑149 (MQ=255)
                                                                                                                              aGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTGGCTCCAGCCTGCCGCCGAAAGAGCGCTGCAAGGCTACCGTCAACCTGGTGCCGGGcgac             <  7:305607/149‑1 (MQ=255)
                                                                                                                               gCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTg                                                                        <  6:45067/89‑1 (MQ=255)
                                                                                                                                ccGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTGGCTACAGCCTGCCGCCGAAAGAGCGCTGCAAGGCTACCGTCAACCTGGTGCCGGGcgacga           <  7:268119/149‑1 (MQ=255)
                                                                                                                                  gagcATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCGTGGGTTGCTCTGGCTACAGCCTGCCGCCGAAAGAGCGCTGCAAGGCTACCGTCAACCTGGTGCCGGGCGTCgaga         <  7:114918/146‑1 (MQ=255)
                                                                                                                                         gaccaCATTCCTTGCACGGAATGTTGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTGGCTACAGCCTGCCGCCGAAAGAGCGCTGCAAGGCTACCGTCAACCTGGTGCCGGGCGACGAGATTGCCg   >  2:157665/5‑149 (MQ=255)
                                                                                                                                          accaCATTCCTTGCACGGAATGTTGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTGGCTACAGCCTGCCGCCGAAAGAGCGCTGCAAGGCTACCGTCAACCTGGTGCCGGGCGACGAGATTGCCGc  >  5:107811/4‑149 (MQ=255)
                                                                                                                                          accaCATTCCTTGCACGGAATGTTGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTGGCTACAGCCTGCCGCCGAAAGAGCGCTGCAAGGCTACCGTCAACCTGGTGCCGGGCGACGAGATTGCCGc  >  4:276545/4‑149 (MQ=255)
                                                                                                                                                ttCCGTGCAAGGAATGTGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTGGCTACAGCCTGCCGCCGAAAGAGCGCTGCAAGGCTAc                                   <  4:88307/109‑1 (MQ=255)
                                                                                                                                                ttCCGTGCAAGGAATGTGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTGGCTACAGCCTGCCGCCGAAAGAGCGCTGCAAGGCTAc                                   >  3:88307/1‑109 (MQ=255)
                                                                                                                                                ttCCGTGCAAGGAATGTGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTGGCTACAGCCTGCCGCCGAAAGAGCGCTGCAAGGCTACCGt                                >  2:218647/1‑112 (MQ=255)
                                                                                                                                                ttCCGTGCAAGGAATGTGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTGGCTACAGCCTGCCGCCGAAAGAGCGCTGCAAGGCTACCGt                                <  1:218647/112‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                          
AAGAGAACCTCGACGACGTCGCCCAGGGCGAGCGTGACTGGAAGAACGTGCTCGACGAGTTCTACGGCGATTTCAGCAAAAAGCTGCAGACTGCCGAGTCCAGCGAAGACGGCATGCGTGCCAACCAGCCGACCATGACCAACATTCCGTGCAAGGAATGTGGTCGGCCGATGATGATTCGCACCGCTTCTACGGGTGTGTTCCTGGGTTGCTCTGGCTACAGCCTGCCGCCGAAAGAGCGCTGCAAGGCTACCGTCAACCTGGTGCCGGGCGACGAGATTGCCGC  >  NC_002947/2442520‑2442805

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: