Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,332,872 C→A 20.9% R609S (CGC→AGC)  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,8720CA20.9% 268.3 / 59.7 114R609S (CGC→AGC) topADNA topoisomerase I, omega subunit
Reads supporting (aligned to +/- strand):  ref base C (43/47);  new base A (13/11);  total (56/58)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

Only 100 of 119 total aligned reads displayed.

GATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAG  >  NC_000913/1332773‑1332971
                                                                                                   |                                                                                                   
gATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTc                                                                                                     >  1:2308820/1‑100 (MQ=255)
 aTAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCg                                                                                                    <  2:1064542/100‑1 (MQ=255)
  tAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGc                                                                                                   <  2:1915143/100‑1 (MQ=255)
  tAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGc                                                                                                   >  1:1323753/1‑100 (MQ=255)
      gCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAg                                                                                               >  1:1467144/1‑100 (MQ=255)
      gCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAg                                                                                               <  1:3559669/100‑1 (MQ=255)
      gCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAg                                                                                               <  1:3302729/100‑1 (MQ=255)
       cTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGc                                                                                              >  2:1708428/1‑100 (MQ=255)
        tGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCg                                                                                             <  2:3480259/100‑1 (MQ=255)
        tGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCg                                                                                             <  2:511511/100‑1 (MQ=255)
        tGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCg                                                                                             <  1:2790045/100‑1 (MQ=255)
        tGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCg                                                                                             >  1:1534193/1‑100 (MQ=255)
         gAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGa                                                                                            >  2:1342116/1‑100 (MQ=255)
         gAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGa                                                                                            >  1:2783176/1‑100 (MQ=255)
         cttccgatctcGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGa                                                                                            <  2:2172126/90‑1 (MQ=255)
           aaaaGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGc                                                                                          >  1:670762/1‑100 (MQ=255)
             aaGATCCGGAAGAGGGTTGTATGCGCGCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGc                                                                                          <  1:1041030/98‑1 (MQ=255)
               gATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCg                                                                                      <  2:197075/100‑1 (MQ=255)
               gATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCa                                                                                         >  1:1011640/1‑97 (MQ=255)
                  ccGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGTCCAGCATTGACTGCCCGACTTGTAGTAGCAAAATGGGGATTCGCACAGCTAGCACCgggg                                                                                   <  1:683358/100‑1 (MQ=255)
                  ccGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCgggg                                                                                   <  2:639655/100‑1 (MQ=255)
                  ccGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCgggg                                                                                   <  2:1897115/100‑1 (MQ=255)
                   cGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGaagatcggaag                                                                                  >  1:2172113/1‑90 (MQ=255)
                   cGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCgggg                                                                                   >  2:934613/1‑99 (MQ=255)
                    ggAAGAGGGTGGTATGCGCCAGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTa                                                                                 >  1:2267523/1‑100 (MQ=255)
                      aaGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTAtt                                                                               <  1:1342110/100‑1 (MQ=255)
                       agagGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTc                                                                              >  2:2836700/1‑100 (MQ=255)
                        gagGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTcc                                                                             <  2:2343639/100‑1 (MQ=255)
                        gagGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTcc                                                                             <  2:856107/100‑1 (MQ=255)
                        gagGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTcc                                                                             <  1:3317424/100‑1 (MQ=255)
                         agGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCt                                                                            <  2:3195751/100‑1 (MQ=255)
                          gggTGGTATGCTCCCGAACCAGAGGGTTCTGACCAGCATTGTCTGCGCGACTTGTGGGCGCAAAATGGGTATTCGCACAGCGAGCACCGGGGTATTCCt                                                                            <  1:2817459/99‑1 (MQ=255)
                          gggTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCtt                                                                           <  1:678345/100‑1 (MQ=255)
                          gggTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCtt                                                                           >  2:718904/1‑100 (MQ=255)
                           ggtggtATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCtt                                                                           <  2:2457104/99‑1 (MQ=255)
                            gtggtATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGGCCGAATTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCAACGGGGTATTCCtt                                                                           >  1:266760/1‑98 (MQ=255)
                            gtggtATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTgg                                                                         <  1:291400/100‑1 (MQ=255)
                               gtATGCGCCCGAACCAGATGGTTCTGACCTGCATTCACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGc                                                                        <  1:2492366/98‑1 (MQ=255)
                                tATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGCGGTATTCCTTGGCTGt                                                                     >  1:2160368/1‑100 (MQ=255)
                                 aTGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGtt                                                                    >  1:201014/1‑100 (MQ=255)
                                 aTGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGtt                                                                    >  2:3123101/1‑100 (MQ=255)
                                   gcgcCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCt                                                                  <  1:3145574/100‑1 (MQ=255)
                                      cccGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGc                                                               >  2:3689671/1‑100 (MQ=255)
                                          aaCCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATg                                                           <  1:1687336/100‑1 (MQ=255)
                                            ccAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATgcg                                                         >  2:684979/1‑100 (MQ=255)
                                            ccAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATgcg                                                         >  2:933453/1‑100 (MQ=255)
                                               gATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTg                                                      <  2:3223632/100‑1 (MQ=255)
                                                 tGGTTCTGACCAGCATTGACTGCCCGCCTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGAACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTgcc                                                    >  1:1330827/1‑100 (MQ=255)
                                                 tGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTgcc                                                    <  1:1128621/100‑1 (MQ=255)
                                                   gTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTgccgc                                                  >  1:3117931/1‑100 (MQ=255)
                                                    ttCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTgccgcc                                                 <  2:2868480/100‑1 (MQ=255)
                                                      cTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGa                                               <  2:1644929/100‑1 (MQ=255)
                                                      cTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGa                                               <  1:70778/100‑1 (MQ=255)
                                                          ccAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAaga                                           <  2:253716/100‑1 (MQ=255)
                                                          ccAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAaga                                           <  1:3478648/100‑1 (MQ=255)
                                                          ccAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAaga                                           >  1:1258242/1‑100 (MQ=255)
                                                              cATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGt                                       >  1:368708/1‑100 (MQ=255)
                                                               aTTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGtt                                      <  1:3292630/100‑1 (MQ=255)
                                                                ttGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGCGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTg                                     <  1:3342632/100‑1 (MQ=255)
                                                                ttGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTg                                     <  1:368782/100‑1 (MQ=255)
                                                                ttGACTGCCCGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTg                                     <  2:442665/100‑1 (MQ=255)
                                                                 tGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTCCTGGCTATGCGCTGCCGCCGAAAGAGcct                                       >  1:3445684/1‑95 (MQ=255)
                                                                  gACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCa                                   >  1:2937752/1‑100 (MQ=255)
                                                                  gACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCa                                   >  1:1177129/1‑100 (MQ=255)
                                                                     tGCCCGACTTGTGGTCGCAAAATGGGGATTCGCGCAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGACTATGCGCTGCCGCCGAAAGAGCATTGCaaaa                                >  1:494993/1‑100 (MQ=255)
                                                                     tGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCaaaa                                <  2:1323759/100‑1 (MQ=255)
                                                                     tGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCaaaa                                <  1:1085905/100‑1 (MQ=255)
                                                                      gCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAac                               >  2:2972696/1‑100 (MQ=255)
                                                                      gCCAGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAACCAGCGTTGCgaa                                 >  1:890896/1‑98 (MQ=255)
                                                                       cccGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAacc                              >  1:2750315/1‑100 (MQ=255)
                                                                        ccGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAacca                             >  1:1303703/1‑100 (MQ=255)
                                                                        ccGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAacca                             >  2:1882543/1‑100 (MQ=255)
                                                                         cGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAaccac                            >  2:2221990/1‑100 (MQ=255)
                                                                         cGACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAacca                             <  1:953756/99‑1 (MQ=255)
                                                                          gACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTAGTCCTTGGCTGTTCGGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAaccacc                           <  1:1473080/100‑1 (MQ=255)
                                                                          gACTTGTGGTCGCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAaccacc                           <  1:1504819/100‑1 (MQ=255)
                                                                             ttGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCAtt                        <  1:2972675/100‑1 (MQ=255)
                                                                              tgtgGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTa                       >  2:1808895/1‑100 (MQ=255)
                                                                                 ggTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAAcc                    >  1:750627/1‑100 (MQ=255)
                                                                                  gTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAAc                     >  2:3394997/1‑98 (MQ=255)
                                                                                  gTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCt                   >  1:2097503/1‑100 (MQ=255)
                                                                                    cttccgatctggATTCGCACAGCGAGCACCGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGt                <  2:1629331/90‑1 (MQ=255)
                                                                                    cGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACAACCATTAACCTg                  >  2:2390537/1‑99 (MQ=255)
                                                                                     gCAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGt                >  2:3259284/1‑100 (MQ=255)
                                                                                      cAAAATGGGGATTAGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTg               >  2:1406031/1‑100 (MQ=255)
                                                                                       aaaaTGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGc              <  2:471783/100‑1 (MQ=255)
                                                                                       aaaaTGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGc              <  2:3055440/100‑1 (MQ=255)
                                                                                        aaaTGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGcc             >  1:931335/1‑100 (MQ=255)
                                                                                         aatgtcgATTAGCACACCGAGCACCGTGGTATTCCTTGGCTCTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCg            <  2:3343066/94‑1 (MQ=255)
                                                                                         aaTGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAAACACCATTAACCTGGTGCCg            >  2:1836569/1‑100 (MQ=255)
                                                                                             gggATTAGCACAGCGAGCACCAGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCAACATTAACCTGGTGCCGGaaa        >  2:3525447/1‑100 (MQ=255)
                                                                                              ggATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGaaaa       <  2:3470685/100‑1 (MQ=255)
                                                                                              ggATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGaaaa       >  1:2306654/1‑100 (MQ=255)
                                                                                              ggATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGaaaa       <  2:2267536/100‑1 (MQ=255)
                                                                                               gATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAAc      >  1:1151408/1‑100 (MQ=255)
                                                                                                aTTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACg     >  1:1416857/1‑100 (MQ=255)
                                                                                                aTTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACg     <  1:3694817/100‑1 (MQ=255)
                                                                                                   cGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAGCCACCATTAACCTGGTGCCGGAAAACGAAg  >  1:2085210/1‑100 (MQ=255)
                                                                                                   cGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAg  >  2:3306474/1‑100 (MQ=255)
                                                                                                   cGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGATATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAg  >  1:2272809/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
GATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAG  >  NC_000913/1332773‑1332971

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: