Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,330,421 A→G L450L (TTA→TTG topA → DNA topoisomerase type I omega protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,330,4210AG98.8% 531.7 / ‑3.0 165L450L (TTA→TTGtopADNA topoisomerase type I omega protein
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (85/78);  total (87/78)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.30e-01

Only 100 of 169 total aligned reads displayed.

CCGCCTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAGCCGCCAGCCCGTTTCAGTGAAGC  >  NC_000913/1330331‑1330510
                                                                                          |                                                                                         
 cGCCTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCca                                                                                        <  1:1680355/93‑1 (MQ=255)
 cGCCTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCca                                                                                        <  2:3477277/93‑1 (MQ=255)
  gCCTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCcag                                                                                       >  2:2154523/1‑93 (MQ=255)
    cgggaaGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCcagca                                                                                     <  2:2306477/89‑1 (MQ=255)
    cTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCcagca                                                                                     >  1:3332239/1‑93 (MQ=255)
    cTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCcagca                                                                                     >  2:2691934/1‑93 (MQ=255)
    cTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCcagca                                                                                     >  1:2762634/1‑93 (MQ=255)
    cTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCcagca                                                                                     <  1:4355507/93‑1 (MQ=255)
    cTGAAAGCACGCGGTAGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCcagca                                                                                     <  1:2219192/93‑1 (MQ=255)
     tGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCcagcag                                                                                    >  2:391165/1‑93 (MQ=255)
     tGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCcagcag                                                                                    >  1:2404647/1‑93 (MQ=255)
      gAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGt                                                                                   <  2:4067595/93‑1 (MQ=255)
      gAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGt                                                                                   >  1:1259806/1‑93 (MQ=255)
       aaaGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGtt                                                                                  <  2:4121691/93‑1 (MQ=255)
        aaGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGTta                                                                                 >  1:2107757/1‑93 (MQ=255)
           cACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGTtaata                                                                              <  2:2089914/93‑1 (MQ=255)
            aCGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGTtaataa                                                                             >  2:654359/1‑93 (MQ=255)
             cgcgGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATaaa                                                                            <  2:1980749/93‑1 (MQ=255)
              gcgGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAg                                                                           >  1:2520915/1‑93 (MQ=255)
              gcgGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAg                                                                           <  2:2986802/93‑1 (MQ=255)
               cgGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAgg                                                                          >  2:1234562/1‑93 (MQ=255)
                ggTCGTATTGTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGc                                                                         <  2:176873/93‑1 (MQ=255)
                 gTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCg                                                                        >  1:3778965/1‑93 (MQ=255)
                 gTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCg                                                                        >  1:3561570/1‑93 (MQ=255)
                  tCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGa                                                                       <  1:3358063/93‑1 (MQ=255)
                   cGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGAt                                                                      <  1:2022802/93‑1 (MQ=255)
                   cGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGAt                                                                      >  2:3033170/1‑93 (MQ=255)
                     tATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGc                                                                    >  2:3179461/1‑93 (MQ=255)
                     tATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGc                                                                    <  1:3172693/93‑1 (MQ=255)
                     tATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGc                                                                    >  1:1471208/1‑93 (MQ=255)
                     tATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGc                                                                    >  2:797178/1‑93 (MQ=255)
                      aTTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGct                                                                   >  2:2932622/1‑93 (MQ=255)
                       ttttGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGctc                                                                  <  1:3239276/93‑1 (MQ=255)
                       ttttGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGctc                                                                  <  1:568348/93‑1 (MQ=255)
                          tGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGa                                                               <  1:3335105/93‑1 (MQ=255)
                           gCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGAc                                                              >  2:4180642/1‑93 (MQ=255)
                           gCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGAc                                                              <  2:247426/93‑1 (MQ=255)
                             gTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGc                                                            <  1:4611445/93‑1 (MQ=255)
                              ttttGATGGCTGGACAAAAGTGATGCCTGCTTTGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCt                                                           <  1:4336921/93‑1 (MQ=255)
                              ttttGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCt                                                           >  1:891105/1‑93 (MQ=255)
                               tttGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTc                                                          <  1:401668/93‑1 (MQ=255)
                               tttGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTc                                                          <  2:3302253/93‑1 (MQ=255)
                                  gATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGtt                                                       <  1:1003957/93‑1 (MQ=255)
                                     ggCTGGACAAAAGCGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGaa                                                    >  1:1761033/1‑93 (MQ=255)
                                      gCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAAc                                                   <  1:1542220/93‑1 (MQ=255)
                                        tGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACtt                                                 >  1:123312/1‑93 (MQ=255)
                                         ggACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTa                                                >  1:4271418/1‑93 (MQ=255)
                                          gACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTac                                               <  2:547533/93‑1 (MQ=255)
                                            cAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTacac                                             >  1:1439856/1‑93 (MQ=255)
                                             aaaaGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGGACTTACAcc                                            <  2:3539028/93‑1 (MQ=255)
                                             aaaaGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACAcc                                            >  1:135666/1‑93 (MQ=255)
                                              aaaGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACAcca                                           >  2:4345185/1‑93 (MQ=255)
                                              aaaGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACAcca                                           <  2:1852774/93‑1 (MQ=255)
                                              aaaGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACAcca                                           <  1:1373914/93‑1 (MQ=255)
                                              aaaGTGATGCCTGCGTGGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACAcca                                           >  2:3129715/1‑93 (MQ=255)
                                               aaGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACAccag                                          >  1:3316512/1‑93 (MQ=255)
                                               aaGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACAccag                                          >  2:2877598/1‑93 (MQ=255)
                                               aaGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACAccag                                          <  1:4227073/93‑1 (MQ=255)
                                               aaGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACAccag                                          >  2:1092211/1‑93 (MQ=255)
                                               aaGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACAccag                                          >  2:1092210/1‑93 (MQ=255)
                                                      gCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCa                                   >  2:2927065/1‑93 (MQ=255)
                                                         tgcgttgcgtAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACtt                                >  2:4476327/1‑93 (MQ=255)
                                                         tgcgttgcgtAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACtt                                <  1:2257096/93‑1 (MQ=255)
                                                            gttgcgtAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTAc                             >  1:4608206/1‑93 (MQ=255)
                                                            gttgcgtAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTAc                             <  1:3926460/93‑1 (MQ=255)
                                                              tgcgtAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCa                           <  1:4542986/93‑1 (MQ=255)
                                                               gcgtAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCaa                          <  1:4633690/93‑1 (MQ=255)
                                                                 gtAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAgc                        <  1:491423/93‑1 (MQ=255)
                                                                 gtAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAgc                        <  2:2733684/93‑1 (MQ=255)
                                                                 gtAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAgc                        >  1:501441/1‑93 (MQ=255)
                                                                 gtAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAgc                        >  2:1365980/1‑93 (MQ=255)
                                                                  tAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAgcc                       >  1:1535032/1‑93 (MQ=255)
                                                                  tAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAgcc                       <  2:323073/93‑1 (MQ=255)
                                                                  tAAAGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAgcc                       <  1:783747/93‑1 (MQ=255)
                                                                    aaGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAgccgc                     >  2:585734/1‑93 (MQ=255)
                                                                    aaGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAgccgc                     >  2:1439331/1‑93 (MQ=255)
                                                                    aaGGCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAgccgc                     >  2:3056589/1‑93 (MQ=255)
                                                                      ggCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAGCCGCCa                   >  2:2762597/1‑93 (MQ=255)
                                                                       gCGATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAGCCGCCAg                  <  1:4262269/93‑1 (MQ=255)
                                                                         gATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAGCCGCCAGcc                >  1:4441248/1‑93 (MQ=255)
                                                                         gATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAGCCGCCAGcc                >  1:1869638/1‑93 (MQ=255)
                                                                         gATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAGCCGCCAGcc                >  1:1939954/1‑93 (MQ=255)
                                                                         gATGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAGCCGCCAGcc                >  1:2433280/1‑93 (MQ=255)
                                                                          aTGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAGCCGACAGccc               >  1:4091778/1‑93 (MQ=255)
                                                                           tGAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAGCCGCCAGCCCg              <  1:2877398/93‑1 (MQ=255)
                                                                            gAAGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAGCCGCCAGCCCGt             <  2:3977352/93‑1 (MQ=255)
                                                                             aaGATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAGCCGCCAGCCCGtt            <  1:642584/93‑1 (MQ=255)
                                                                              aGATCGCATCTTGCCAGCAGTTATTAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAGCCGCCAGCCCGttt           <  2:2586081/93‑1 (MQ=255)
                                                                               gATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAGCCGCCAGCCCGTTTc          <  2:3693942/93‑1 (MQ=255)
                                                                               gATCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAGCCGCCAGCCCGTTTc          >  1:429145/1‑93 (MQ=255)
                                                                                aTCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAGCCGCCAGCCCGTTTCa         >  2:3356166/1‑93 (MQ=255)
                                                                                 tCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAGCCGCCAGCCCGTTTCAg        >  1:12099/1‑93 (MQ=255)
                                                                                 tCGCATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAGCCGCCAGCCCGTTTCAg        >  2:4407232/1‑93 (MQ=255)
                                                                                  cgcgTCTTGCCAGCAGTTAATAAAGGCGATCCTCTGACGCTCGTTGAACTTGCACCAGCCCAGCACTTTACCAAGCCGCCAGCCCGTTTCAGt       <  2:4230218/93‑1 (MQ=255)
                                                                                    cATCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAGCCGCCAGCCCGTTTCAGTGa     >  1:2825847/1‑93 (MQ=255)
                                                                                     aTCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAGCCGCCAGCCCGTTTCAGTGag    >  1:3792608/1‑92 (MQ=255)
                                                                                     aTCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAGCCGCCAGCCCGTTTCAGTGaa    >  1:1458259/1‑93 (MQ=255)
                                                                                      tCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAGCCGCCAGCCCGTTTCAGTGAAg   <  2:1860164/93‑1 (MQ=255)
                                                                                      tCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAGCCGCCAGCCCGTTTCAGTGAAg   >  2:1168358/1‑93 (MQ=255)
                                                                                      tCTTGCCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAGCCGCCAGCCCGTTTCAGTGAAg   <  2:935583/93‑1 (MQ=255)
                                                                                          |                                                                                         
CCGCCTGAAAGCACGCGGTCGTATTTTGCGTTTTGATGGCTGGACAAAAGTGATGCCTGCGTTGCGTAAAGGCGATGAAGATCGCATCTTACCAGCAGTTAATAAAGGCGATGCTCTGACGCTCGTTGAACTTACACCAGCCCAGCACTTTACCAAGCCGCCAGCCCGTTTCAGTGAAGC  >  NC_000913/1330331‑1330510

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