Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,332,633 G→T R529L (CGC→CTC)  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,6330GT96.7% 318.6 / ‑2.6 90R529L (CGC→CTC) topADNA topoisomerase I, omega subunit
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base T (41/46);  total (42/48)
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

TATGCGTCGATCATTTCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCGCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAA  >  NC_000913/1332536‑1332729
                                                                                                 |                                                                                                
tATGCGTCGATCATTTCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTcc                                                                                                <  2:1389827/100‑1 (MQ=255)
 aTGCGTCGATCATTTCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCt                                                                                               <  2:1211644/100‑1 (MQ=255)
   gCGTCGATCATTTCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTg                                                                                             >  2:163791/1‑100 (MQ=255)
   gCGTCGATCATTTCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTg                                                                                             <  2:1130707/100‑1 (MQ=255)
     gTCGATCATTTCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTgaa                                                                                           >  2:555929/1‑100 (MQ=255)
     gTCGATCATTTCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAAGGGGCGAAATCGTCACCGATCTCCTTga                                                                                            <  2:795499/99‑1 (MQ=255)
      tCGATCATTTCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTgaag                                                                                          >  2:1984698/1‑100 (MQ=255)
       cGATCATTTCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTgaaca                                                                                         <  1:1626729/100‑3 (MQ=255)
        gATCATTTCGACCATTCAGGATCGTGGCTACGTGCGGGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTgaagaa                                                                                        >  1:1696710/1‑100 (MQ=255)
         atgaTTTCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAATCTTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGaaa                                                                                       <  2:98793/97‑1 (MQ=255)
         aTCATTTCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGaaa                                                                                       >  1:2392963/1‑100 (MQ=255)
            aTTTCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAAtt                                                                                    >  2:1801384/1‑100 (MQ=255)
             tttCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAAttt                                                                                   >  2:2036502/1‑100 (MQ=255)
              ttCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTc                                                                                  <  1:1959063/100‑1 (MQ=255)
              ttCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTc                                                                                  <  1:12910/100‑1 (MQ=255)
               tCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTcc                                                                                 <  1:1801383/100‑1 (MQ=255)
                cGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCcg                                                                                >  2:1818694/1‑100 (MQ=255)
                  aCCATTCAGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCcgcg                                                                              <  1:736790/100‑1 (MQ=255)
                   ccATTCAGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGa                                                                             >  2:1651173/1‑100 (MQ=255)
                   ccATTCAGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGa                                                                             <  1:441206/100‑1 (MQ=255)
                      ttCAGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGtt                                                                          <  1:1745302/100‑1 (MQ=255)
                       tCAGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTa                                                                         <  2:1696711/100‑1 (MQ=255)
                         aGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAAGTTCCGCGAGTTaa                                                                        >  2:2148935/1‑99 (MQ=255)
                          ggATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATg                                                                      >  2:1138954/1‑100 (MQ=255)
                          ggATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATg                                                                      <  1:1607861/100‑1 (MQ=255)
                            aTCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGaa                                                                    <  2:1029157/100‑1 (MQ=255)
                             tCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAAc                                                                   <  2:2392960/100‑1 (MQ=255)
                               gTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTa                                                                 >  1:1233224/1‑100 (MQ=255)
                                  gCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGa                                                              >  1:1241857/1‑100 (MQ=255)
                                   cTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGAc                                                             <  1:163792/100‑1 (MQ=255)
                                   cTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAAGCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGAc                                                             <  1:1505247/100‑1 (MQ=255)
                                    tACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACt                                                            >  1:1554068/1‑100 (MQ=255)
                                     aCGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACtt                                                           >  1:1845769/1‑100 (MQ=255)
                                       gTGCGAGTAGAAAATTGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTGCTTGAAGAAAATTTCCGAGAGTTAATGAACTACGACTTTa                                                         >  1:864532/1‑100 (MQ=255)
                                       gTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTa                                                         <  1:773192/100‑1 (MQ=255)
                                        tGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTAc                                                        >  1:444529/1‑100 (MQ=255)
                                         gCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTAcc                                                       <  2:2002556/100‑1 (MQ=255)
                                         gCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTAcc                                                       >  1:1857473/1‑100 (MQ=255)
                                            aGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACcgcg                                                    <  1:960313/100‑1 (MQ=255)
                                             gTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCgcgc                                                   <  2:562049/100‑1 (MQ=255)
                                              tAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCa                                                  <  2:1857473/100‑1 (MQ=255)
                                               aGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAg                                                 >  1:10761/1‑100 (MQ=255)
                                                gAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGa                                                >  1:1998512/1‑100 (MQ=255)
                                                gAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGa                                                >  1:1940509/1‑100 (MQ=255)
                                                  aaaTCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATg                                              <  1:1397716/100‑1 (MQ=255)
                                                  aaaTCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATg                                              >  1:2322442/1‑100 (MQ=255)
                                                  aaaTCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATg                                              >  2:1726147/1‑100 (MQ=255)
                                                  aaaTCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATg                                              >  1:988027/1‑100 (MQ=255)
                                                    atcttcgtTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGa                                            <  2:1195116/96‑1 (MQ=255)
                                                      cgtcgtTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGaaa                                          <  2:1884806/100‑1 (MQ=255)
                                                        tcgtTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGaagat                                        >  1:1195114/1‑97 (MQ=255)
                                                        tcgtTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAAc                                        <  1:1279544/100‑1 (MQ=255)
                                                         cgtTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACa                                       >  2:1192566/1‑100 (MQ=255)
                                                          gtTTCTATGCGGACAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAg                                      >  2:1482169/1‑100 (MQ=255)
                                                          gtTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAg                                      <  1:557367/100‑1 (MQ=255)
                                                            ttCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGCCTTTACCGCGCAGATGGAAAACAGcc                                    >  2:1992934/1‑100 (MQ=255)
                                                               tATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTATGACTTTACCGCGCAGATCGAAAACAGCCttg                                 <  1:1618648/100‑3 (MQ=255)
                                                               tATGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCg                                 <  1:926001/100‑1 (MQ=255)
                                                                aTGCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGa                                >  1:1315209/1‑100 (MQ=255)
                                                                  gCGGAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGAcc                              <  2:154325/100‑1 (MQ=255)
                                                                  gCGGAAAAAATGGGCGAAATCGTCACCGATCGCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGAcc                              <  1:1669903/100‑1 (MQ=255)
                                                                    ggAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAg                            >  1:1085647/1‑100 (MQ=255)
                                                                    ggAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAg                            >  1:2061893/1‑100 (MQ=255)
                                                                    ggAAAAAATGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAg                            >  1:661037/1‑100 (MQ=255)
                                                                      aaaaaaTGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGt                          <  1:2098140/100‑1 (MQ=255)
                                                                      aaaaaaTAGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGt                          >  1:1268827/1‑100 (MQ=255)
                                                                       aaaaaTGGGCGAAATCGTCACCGATCTCCTTGAAGAAAAATTACGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTg                         <  2:1685724/100‑1 (MQ=255)
                                                                        aaaaTGGGCGAATTCGTCATAGATCGCCTTGAAGAAAATGTCCGCGGGTTATTGAACCACGACTTGACCGCGCAGATGGAAAAGAGCCTCGCCCAGGTgg                        <  2:2364721/100‑1 (MQ=255)
                                                                         aaaTGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGc                       >  2:997548/1‑100 (MQ=255)
                                                                          aaTGGGCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCa                      >  2:1942041/1‑100 (MQ=255)
                                                                             gggCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAAt                   <  1:954715/100‑1 (MQ=255)
                                                                              ggCGAAATCGTCACCGATCTCCTTGAATAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATc                  <  1:250675/100‑1 (MQ=255)
                                                                              ggCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATc                  >  2:437174/1‑100 (MQ=255)
                                                                              ggCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATc                  <  1:1983980/100‑1 (MQ=255)
                                                                               gCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCa                 <  1:1192564/100‑1 (MQ=255)
                                                                               gCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCa                 <  1:969509/100‑1 (MQ=255)
                                                                               gCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCa                 >  1:191878/1‑100 (MQ=255)
                                                                               gCGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCa                 >  2:300730/1‑100 (MQ=255)
                                                                                cGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCAc                <  2:1554069/100‑1 (MQ=255)
                                                                                cGAAATCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCAc                <  2:1845769/100‑1 (MQ=255)
                                                                                  aaaTCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGa              <  1:661688/100‑1 (MQ=255)
                                                                                   aaTCGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGaa             <  2:10761/100‑1 (MQ=255)
                                                                                      cGTCACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCa          <  2:1073192/100‑1 (MQ=255)
                                                                                      cGTCACCGATCTCCTTGAAGAAAATGTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCa          >  2:401637/1‑100 (MQ=255)
                                                                                         cACCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCagag       <  2:444529/100‑1 (MQ=255)
                                                                                          aCCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGt      <  2:2061891/100‑1 (MQ=255)
                                                                                          aCCGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGt      >  1:1562926/1‑100 (MQ=255)
                                                                                           ccGATCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTg     <  1:1984699/100‑1 (MQ=255)
                                                                                           ccGATCGCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTg     >  2:1035688/1‑100 (MQ=255)
                                                                                              aTCTCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGaa  >  2:83816/1‑100 (MQ=255)
                                                                                                 |                                                                                                
TATGCGTCGATCATTTCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCGCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAA  >  NC_000913/1332536‑1332729

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