Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,331,126 G→T 57.4% V27L (GTG→TTG)  topA → DNA topoisomerase I, omega subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,331,1260GT57.4% 37.3 / 150.9 106V27L (GTG→TTG) topADNA topoisomerase I, omega subunit
Reads supporting (aligned to +/- strand):  ref base G (21/24);  new base T (25/36);  total (46/60)
Fisher's exact test for biased strand distribution p-value = 6.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.41e-01

Only 100 of 109 total aligned reads displayed.

CAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAG  >  NC_000913/1331029‑1331221
                                                                                                 |                                                                                               
cAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCTGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACgtggtg                                                                                               >  2:1159910/1‑100 (MQ=255)
cAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACgtgttg                                                                                               <  2:1841701/100‑1 (MQ=255)
cAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACgtgttg                                                                                               <  2:3246787/100‑1 (MQ=255)
 aaaTTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGTTGa                                                                                              <  2:708258/100‑1 (MQ=255)
    ttAGGTAAGGTGAAGATGTGTAAAGCTCTGTTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGaaan                                                                                           <  1:1369319/100‑2 (MQ=255)
      aGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATcc                                                                                         <  1:906338/100‑1 (MQ=255)
       ggTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGTTGAAATCCa                                                                                        <  1:168095/100‑1 (MQ=255)
        gTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAg                                                                                       <  2:695191/100‑1 (MQ=255)
         tAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGTTGAAATCCAGc                                                                                      >  2:2008808/1‑100 (MQ=255)
         tAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGTTGAAATCCAGc                                                                                      <  1:3478875/100‑1 (MQ=255)
         tAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGTTGAAATCCAGc                                                                                      <  1:3474105/100‑1 (MQ=255)
            ggTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGTTGAAATCCAGCGTc                                                                                   <  2:949245/100‑1 (MQ=255)
            ggTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTc                                                                                   >  2:1549818/1‑100 (MQ=255)
              tGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCgg                                                                                 <  1:2101080/100‑1 (MQ=255)
                aaTATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTc                                                                               >  2:3372738/1‑100 (MQ=255)
                aaTATGGGTAAAGCTCTTGTAATCGTTGAGTCCCTGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGTTGAAATCCAGCGTCGGTc                                                                               <  1:1460680/100‑1 (MQ=255)
                 atctGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTca                                                                              <  1:3014144/97‑1 (MQ=255)
                 atatGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGTTGAAATCCAGCGTCGGTca                                                                              >  2:435520/1‑100 (MQ=255)
                    tGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGTTGAAATCCAGCGTCGGTCACAt                                                                           <  1:1418743/100‑1 (MQ=255)
                     gggTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGTTGAAATCCAGCGTCGGTCACATc                                                                          >  2:189951/1‑100 (MQ=255)
                     gggTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTcaaga                                                                           >  2:3014166/1‑96 (MQ=255)
                      ggTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGTTGAAATCCAGCGTCGGTCACATcc                                                                         >  2:1064425/1‑100 (MQ=255)
                      ggTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGTTGAAATCCAGCGTCGGTCACATc                                                                          >  1:294684/1‑99 (MQ=255)
                       gTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGGCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCcg                                                                        >  2:2615774/1‑100 (MQ=255)
                          aaGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGTTGAAATCCAGCGTCGGTCACATCCGCGa                                                                     <  1:81443/100‑1 (MQ=255)
                           aGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGAt                                                                    <  2:2857458/100‑1 (MQ=255)
                           aGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGAt                                                                    <  1:144164/100‑1 (MQ=255)
                           aGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGAGCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGAt                                                                    <  1:606764/100‑1 (MQ=255)
                            gCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGTTGAAATCCAGCGTCGGTCACATCCGCGAtt                                                                   <  2:3114498/100‑1 (MQ=255)
                            gCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGAtt                                                                   <  2:1621281/100‑1 (MQ=255)
                             ctctTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGTTGAAATCCAGCGTCGGTCACATCCGCGAttt                                                                  <  1:3499958/100‑1 (MQ=255)
                             ctctTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGTTGAAATCCAGCGTCGGTCACATCCGCGAttt                                                                  >  1:2924276/1‑100 (MQ=255)
                              tctTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGTTGAAATCCAGCGTCGGTCACATCCGCGATTTg                                                                 <  2:383870/100‑1 (MQ=255)
                              tctTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGTTGAAATCCAGCGTCGGTCACATCCGCGATTTg                                                                 >  1:824874/1‑100 (MQ=255)
                              tctTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTg                                                                 <  2:625362/100‑1 (MQ=255)
                               ctTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGTTGAAATCCAGCGTCGGTCACATCCGCGATTTGc                                                                <  2:1325395/100‑1 (MQ=255)
                                ttGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGcc                                                               <  2:2270545/100‑1 (MQ=255)
                                 tGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGTTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCg                                                              >  2:1491987/1‑100 (MQ=255)
                                 tGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCg                                                              >  2:2096255/1‑100 (MQ=255)
                                 tGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCg                                                              <  1:828938/100‑1 (MQ=255)
                                  gTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGTTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGa                                                             <  1:2490829/100‑1 (MQ=255)
                                    cATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGTTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGAcc                                                           >  1:1842974/1‑100 (MQ=255)
                                    cATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGAcc                                                           <  2:1427629/100‑1 (MQ=255)
                                    cATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGAcc                                                           <  1:522462/100‑1 (MQ=255)
                                     aTCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGTTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCa                                                          >  1:2036592/1‑100 (MQ=255)
                                     aTCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGTTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCa                                                          >  2:258263/1‑100 (MQ=255)
                                      tCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTTTCTGGGTAGTGACTACGTGTTGAAATCCAGCGTCGGTCACATCCGCGATTTGTCGACCAg                                                         <  1:1237432/100‑1 (MQ=255)
                                      tCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGTTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAg                                                         <  1:2541097/100‑1 (MQ=255)
                                        gTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATTTGGGTAGTGACTACGTGTTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTg                                                       <  1:461827/100‑1 (MQ=255)
                                        gTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGTTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTg                                                       >  2:63361/1‑100 (MQ=255)
                                            tctCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGTTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTc                                                   <  1:3223295/98‑1 (MQ=255)
                                             gTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGTTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCa                                                  >  1:1311710/1‑100 (MQ=255)
                                             gTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGTTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCa                                                  >  1:2485785/1‑100 (MQ=255)
                                             gTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCa                                                  >  2:1203871/1‑100 (MQ=255)
                                              tCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTTGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAg                                                 <  2:2232736/100‑1 (MQ=255)
                                               ccccGGCAAAAGCCAAAACGATCAGCAAGTATCTGGGTAGTGCCTACTTGTCGAAATCCAGCGTCGGTCACATCCGCGATTTGCTGACCAGTGGCTc                                                   <  1:2583573/97‑1 (MQ=255)
                                               ccccGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGTTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAga                                                >  2:3223319/1‑99 (MQ=255)
                                               ccccGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGTTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGc                                                >  2:2751658/1‑100 (MQ=255)
                                               ccccGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGc                                                <  2:3334147/100‑1 (MQ=255)
                                                cccGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGTTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCt                                               >  2:3022318/1‑100 (MQ=255)
                                                cccGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGTTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCt                                               <  1:63361/100‑1 (MQ=255)
                                                cccGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGTTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCt                                               <  2:2436869/100‑1 (MQ=255)
                                                     cAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaa                                          >  2:3389593/1‑100 (MQ=255)
                                                     cAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACATGTTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaa                                          <  1:1983168/100‑1 (MQ=255)
                                                        aaGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGTTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaaaaa                                       <  2:2864246/100‑1 (MQ=255)
                                                         aGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGTTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAag                                      >  2:1429183/1‑100 (MQ=255)
                                                          gCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGTTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAaga                                     >  2:2270594/1‑100 (MQ=255)
                                                          gCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAaga                                     <  1:1539632/100‑1 (MQ=255)
                                                          gCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAaga                                     >  1:1198940/1‑100 (MQ=255)
                                                           ccAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGTTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAagag                                    >  2:899950/1‑100 (MQ=255)
                                                           ccAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAagag                                    <  1:1159904/100‑1 (MQ=255)
                                                            cAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCGGCTGCCAAAAAGAGt                                   >  2:3646679/1‑100 (MQ=255)
                                                            cAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGt                                   <  2:3322374/100‑1 (MQ=255)
                                                            cAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGt                                   >  1:103898/1‑100 (MQ=255)
                                                             aaaaCGATCAACAAGTATCTGGGTAGTGACTACGTGTTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTg                                  <  2:673200/100‑1 (MQ=255)
                                                                aCGATCAACAAGTATCTGGGTAGTGACTACGTGTTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCg                               <  2:2080636/100‑1 (MQ=255)
                                                                  gATCAACAAGTATCTGGGTAGTGACTACGTGTTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCg                               >  2:2297187/1‑98 (MQ=255)
                                                                  gATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGAc                             >  2:2195661/1‑100 (MQ=255)
                                                                      aacaaGTATCTGGGTAGTGACTACGTGTTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGActct                          <  1:3544865/99‑1 (MQ=255)
                                                                      aacaaGTATCTGGGTAGTGACTACGTGTTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTa                         >  2:1213153/1‑100 (MQ=255)
                                                                      aacaaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTa                         >  2:1654238/1‑100 (MQ=255)
                                                                      aacaaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTa                         >  2:1349982/1‑100 (MQ=255)
                                                                      aacaaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTa                         <  2:3243658/100‑1 (MQ=255)
                                                                      aacaaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTa                         >  2:313493/1‑100 (MQ=255)
                                                                        caaGTATCTGGGTAGTGACTACGTGTTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTAcc                       <  1:1314374/100‑1 (MQ=255)
                                                                        caaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTAcc                       <  1:3434417/100‑1 (MQ=255)
                                                                         aaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCt                      <  1:3451204/100‑1 (MQ=255)
                                                                          aGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTc                     >  2:2300707/1‑100 (MQ=255)
                                                                            tATCTGGGGAGTTACTACGTGTTGAAATCCAGCGTCGGTCACAACCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTcc                    <  1:3004116/99‑1 (MQ=255)
                                                                             aTCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTccac                  <  2:1225907/100‑1 (MQ=255)
                                                                                tGGGTAGTGACTACGTGTTGAAATACAGCGTCGATCACATCCGCGATTTGCCGACCAGTGGCTCAGCAACCAAAAAGAGTGCCGACTCTACCTCCACCaa               <  1:2955670/100‑1 (MQ=255)
                                                                                 nggTAGTGACTACGTGTTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAg              >  2:1877409/2‑100 (MQ=255)
                                                                                      gTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAACGAGTGCCGACGCTACCTCCTCCAAGACg           >  1:3142044/1‑98 (MQ=255)
                                                                                      gTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGc         >  2:1162697/1‑100 (MQ=255)
                                                                                       tGACTACGTGTTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCt        <  1:258263/100‑1 (MQ=255)
                                                                                        gACTACGTGTTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTa       >  2:2838254/1‑100 (MQ=255)
                                                                                        gACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTa       >  1:729750/1‑100 (MQ=255)
                                                                                             cGTGTTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAg  <  2:700436/100‑1 (MQ=255)
                                                                                             cGTGTTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAg  <  1:3651010/100‑1 (MQ=255)
                                                                                             cGTGTTGAAATCCAGCGTCGGTCACATCCGCGATTCGCCGACCAGTGGCGCAGCTGCCACAAAGAGGGCCCACTCTACCTCCACTAGCACGGc         >  1:2904642/1‑93 (MQ=255)
                                                                                                 |                                                                                               
CAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAG  >  NC_000913/1331029‑1331221

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