Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 350,119 T→G 77.3% G36G (GGT→GGG prpC → 2‑methylcitrate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913350,1190TG77.3% 232.7 / 60.0 117G36G (GGT→GGGprpC2‑methylcitrate synthase
Reads supporting (aligned to +/- strand):  ref base T (12/14);  new base G (44/44);  total (59/58)
Fisher's exact test for biased strand distribution p-value = 8.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.18e-01

Only 100 of 120 total aligned reads displayed.

CAACGATCCTGCAAAACAGTACCCATGTCATTAAACCGAAAAAATCTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGTAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATTCTTGATCTGGCGAAACATTGCGAATTTGAAGAAGTGGCGCATCTGCTGATCCACGGC  >  NC_000913/350022‑350215
                                                                                                 |                                                                                                
caACGATCCTGCAAAACAGTACCCATGTCATTAAACCGAAAAAATCTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGTaa                                                                                                >  2:2222196/1‑100 (MQ=255)
caACGATCCTGCAAAACAGTACCCATGTCATTAAACCGAAAAAATCTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGGaa                                                                                                >  2:1248530/1‑100 (MQ=255)
   cGATCCTGCAAAACAGTACCCATGTCATTAAACCGAAAAAATCTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGGAAAAg                                                                                             >  2:3666479/1‑100 (MQ=255)
   cGATCCTGCAAAACAGTACCCATGTCATTAAACCGAAAAAATCTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGGAAAAg                                                                                             <  1:3457312/100‑1 (MQ=255)
   cGATCCTGCAAAACAGTACCCATGTCATTAAACCGAAAAAATCTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGAAAAGt                                                                                            >  2:2512943/1‑100 (MQ=255)
    gATCCTGCAAAACAGTACCCATGTCATTAAACCGAAAAAATCTGTGGGACTTTCTGGCGTTCCGGCGGGAAATACGGCGCTCTGCACCGTGGGGAAACGt                                                                                            >  2:2464883/1‑100 (MQ=255)
       tcTGCAAAACAGTACCCATGTCATTAAACCGAAAAAATCTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGGAAAAGTGGc                                                                                         <  1:4150480/99‑1 (MQ=255)
        cTGCAAAACATTACATATTTCATTAAACCGAAAAAATCTGTGGCACTTTCGGGCTTTCCGGCGGGCAATACGGCGCTCTGCCCCGTGGGTAAACGAGGCa                                                                                        >  2:4129556/1‑100 (MQ=255)
        cTGCAAAACAGTACCCATGTCATTAAACCGAAAAAATCTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGTAAAAGTGGCa                                                                                        <  2:3389441/100‑1 (MQ=255)
        cTGCAAAACAGTACCCATGTCATTAAACCGAAAAAATCTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGGAAAAGTGGCa                                                                                        >  2:4150446/1‑100 (MQ=255)
         tGCAAAACAGTACCCATGTCATTAAACCGAAAAAATCTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGGAAAAGTGGCaa                                                                                       >  2:871394/1‑100 (MQ=255)
         tGCAAAACAGTACCCATGTCATTAAACCGAAAAAATCTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGGAAAAGTGGCaa                                                                                       >  2:2848273/1‑100 (MQ=255)
         tGCAAAACAGTACCCATGTCATTAAACCGAAAAAATCTGTGGCACTCTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGCGGGGAAAAGTGGCa                                                                                        >  2:1406737/1‑99 (MQ=255)
          gCAAAACAGTACCCATGTCATTAAACCGAAAAAATCTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGGAAAAGTGGCAAt                                                                                      <  1:1211867/100‑1 (MQ=255)
          gCAAAACAGTACCCATGTCATTAAACCGAAAAAATCTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGGAAAAGTGGCAAt                                                                                      >  2:173193/1‑100 (MQ=255)
           cAAAACAGTACCCATGTCATTAAACCGAAAAAATCTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGTAAAAGTGGCAATg                                                                                     <  1:590688/100‑1 (MQ=255)
           cAAAACAGTACCCATGTCATTAAACCGAAAAAATCTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGGAAAAGTGGCAATg                                                                                     <  1:2773519/100‑1 (MQ=255)
             aaaCAGTACCCATGTCATTAAACCGAAAAAATCTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGTAAAAGTGGCAATGAc                                                                                   >  2:1940162/1‑100 (MQ=255)
             aaaCAGTACCCATGTCATTAAACCGAAAAAATCTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGGAAAAGTGGCAATGAc                                                                                   <  1:2904845/100‑1 (MQ=255)
              aaCAGTATCCATGTCATTAAACCGAAAAAATCTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGTAAAAGTGGCAATGAcc                                                                                  <  2:2232691/100‑1 (MQ=255)
                 aGTACCCATGTCATTAAACCGAAAAAATCTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGGAAAAGTGGCAATGACCTGc                                                                               <  2:413060/100‑1 (MQ=255)
                 aGTACCCATGTCATTAAACCGAAAAAATCTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGGAAAAGTGGCAATGACCTGc                                                                               <  2:1691167/100‑1 (MQ=255)
                  gTACCCATGTCATTAAACCGAAAAAATCTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGTAAAAGTGGCAATGACCTGCa                                                                              <  2:2060355/100‑1 (MQ=255)
                    aCCCATGTCATTAAACCGAAAAAATCTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGGAAAAGTGGCAATGACCTGCAtt                                                                            >  2:3654533/1‑100 (MQ=255)
                     cccATGTCATTAAACCGAAAAAATCTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGGAAAAGTGGCAATGACCTGCATTa                                                                           >  2:2962444/1‑100 (MQ=255)
                         tGTCATTAAACCGAAAAAATCTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGTAAAAGTGGCAATGACCTGCATTACcgc                                                                       <  2:1611561/100‑1 (MQ=255)
                          gTCATTAAACCGAAAAAATCTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGGAAAAGTGGCAATGACCTGCATTACcgcg                                                                      <  2:3727048/100‑1 (MQ=255)
                          gTCATTAAACCGAAAAAATCTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGGAAAAGTGGCAATGACCTGCATTACcgcg                                                                      <  2:2683013/100‑1 (MQ=255)
                           tCATTAAACCGCAAAAATCTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGGAAACGTGGCAATGACCTGCATTACcgc                                                                       >  2:2011759/1‑98 (MQ=255)
                           tCATTAAACCGAAAAAATCTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGGAAAAGTGGCAATGACCTGCATTACcgcg                                                                      >  2:2315469/1‑99 (MQ=255)
                            cATTAAACCGAAAAAATCTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGGAAAAGTGGCAATGACCTGCATTACCGCGGc                                                                    >  2:1045890/1‑100 (MQ=255)
                            cATTAAACCGAAAAAATCTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGGAAAAGTGGCAATGACCTGCATTACCGCGGc                                                                    <  1:874109/100‑1 (MQ=255)
                             aTTAAACCGAAAAAATCTGTGGCATTTTCTCGCGTTCCGCCGGGCACTACCGCGCGCTGCACCGTGGCGTAAACTGGCAAGCCCGTGCATTAGcga                                                                       >  2:225353/1‑95 (MQ=255)
                             aTTAAACCGAAAAAATCTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGGAAAAGTGGCAATGACCTGCATTACCGCGGCt                                                                   <  1:540765/100‑1 (MQ=255)
                             aTTAAACCGAAAAAATCTGTGGCACTTTCTGGCGTTCCGGCAGGCAATACGGCGCGCTGCACCGTGGGGAAAAGTGCAAATCACCCGCATTACCGCGGc                                                                    >  2:1451236/1‑99 (MQ=255)
                             aTTAAACCGAAAAAAACTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGGAAAAGTGGCAATGACCTGCATTACCGCGGCt                                                                   >  1:3689118/1‑100 (MQ=255)
                              ttAAACCGAAAAAATCTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGGAAAAGTGGCAATGACCTGCATTACCGCGGCTa                                                                  <  2:3835446/100‑1 (MQ=255)
                              ttAAACCGAAAAAATCTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGACAAAGTGGCAATGACCTGCATTACCGCGGCTa                                                                  >  1:2137417/1‑100 (MQ=255)
                               tAAACCGAAAAAATCTGTGGCACTTTCTGGCGTTCCGGCGGACAATACGGCGCTCTGCACCGAGGGGAAAAGTGGCAATGACCTGCATTACCGCGGCTAc                                                                 >  2:470534/1‑100 (MQ=255)
                                 aaCCGAAAAAATCTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGGAAAAGTGGCAATGACCTGCATTACCGCGGCTACGa                                                               <  1:1370689/100‑1 (MQ=255)
                                  aCCGAACAAATATGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGGAAAAGTGGCAATGACCTGCATTACCGCGGCTACGat                                                              <  1:776473/100‑1 (MQ=255)
                                  aCCGAAAAAATCTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGTAAAAGTGGCAATGACCTGCATTACCGCGGCTACGat                                                              >  2:2955746/1‑100 (MQ=255)
                                  aCCGAAAAAATCTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGGAAAAGTGGCAATGACCTGCATTACCGCGGCTACGat                                                              <  1:3071861/100‑1 (MQ=255)
                                   ccGAAAAAATCTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGGAAAAGTGGCAATGACCTGCATTACCGCGGCTACGata                                                             <  1:2166626/100‑1 (MQ=255)
                                    cGAAAAAATCTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGGAAAAGTGGCAATGACCTGCATTACCGCGGCTACGatat                                                            >  2:491752/1‑100 (MQ=255)
                                    cGAAAAAATCTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGCGGGGCAAAGTGGCAATGACCTGCATTACCGCGGCTACGata                                                             >  2:2415102/1‑99 (MQ=255)
                                       aaaaaTCTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGTAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATTCt                                                         <  1:1716884/100‑1 (MQ=255)
                                       aaaaaTCTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGGAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATTCt                                                         <  1:2495592/100‑1 (MQ=255)
                                        aaaaTCTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGGAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATTCtt                                                        >  2:2934034/1‑100 (MQ=255)
                                        aaaaTCTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGGAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATTCtt                                                        <  2:4062259/100‑1 (MQ=255)
                                         aaaTCTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGTAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATTCTTg                                                       <  2:1346751/100‑1 (MQ=255)
                                         aaaGCTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGGAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATTCTTg                                                       <  1:1110046/100‑1 (MQ=255)
                                            tCTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGGAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATTCTTGATc                                                    <  1:3674459/100‑1 (MQ=255)
                                            tCTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGGAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATTCTTGATc                                                    <  2:3324055/100‑1 (MQ=255)
                                             cTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGTAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATTCTTGATCt                                                   <  2:3360731/100‑1 (MQ=255)
                                              tgtgGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGGAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATTCTTGATCTg                                                  <  1:148436/100‑1 (MQ=255)
                                               gtgGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGGAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATTCTTGATCTgg                                                 <  1:634002/100‑1 (MQ=255)
                                                tgGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGTAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATGCTAGATCTCGc                                                >  2:3828559/1‑100 (MQ=255)
                                                 ggCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGGAAACGTGGCAATGACCTGCCTCACCGCGGCTACGCCATTCTTGATCCGGc                                                >  1:2785641/1‑99 (MQ=255)
                                                 ggCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGGAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATTCTTGATCTGGCg                                               >  1:3825286/1‑100 (MQ=255)
                                                 ggCACTTTCTGGCGTTCCGGCGGGCAATACGGCACTCTGCACCGTGGGTAAAAGTGGCAATGCCCTGCATTACCGCGGCTACGATATTCTTGATCTGGc                                                >  2:2757892/1‑99 (MQ=255)
                                                  gCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGGAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATTCTTGATCTGGCGa                                              <  2:3705105/100‑1 (MQ=255)
                                                   cACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTTCACCGTGGGTAAAAGTGGCTATGACCTGCATTACCGCGGCTACGATATTCTTGATCTGGCGaa                                             >  2:1650427/1‑100 (MQ=255)
                                                        tCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGGAAAAGTGGCAATTACCTGCATTTCCGCGGCTACGATATTCTTTATCCGGGGAAAAAtt                                        >  2:1124324/1‑100 (MQ=255)
                                                         cTGGCGTTCCGGCGGGCAATATGGCGCTCTGCACCGTGGGGAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATTCTTGATCTGGCGAAACATTg                                       <  1:67519/100‑1 (MQ=255)
                                                          tGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGTGAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATTCTTGATCTGGCGAAACATTGc                                      <  2:2265086/100‑1 (MQ=255)
                                                             cGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGTAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATTCTTGATCTGGCGAAACATTGCGa                                    >  2:2683359/1‑99 (MQ=255)
                                                              gTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGGAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATTCTTGATCTGGCGAAACATTGCGAAt                                  <  1:2001302/100‑1 (MQ=255)
                                                              gTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGGAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATTCTTGATCTGGCGAAACATTGCGAAt                                  <  1:2417776/100‑1 (MQ=255)
                                                              gTTCCGGCGGGCAATACCGCGCTCTGCACCGTGGGGAACAGTGGCAATGACCTGCATTACCCCGGCTACGATATTCCGCATCTGGCGAAACATTGCGaa                                   >  2:866486/1‑99 (MQ=255)
                                                               ttCCGGCGGGCAATACGGCGCTCTGCACCGTGGGGAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATTCTTGATCTGGCGAAACATTGCGAAtt                                 <  2:762875/100‑1 (MQ=255)
                                                                tCCGGCGGGCAATACGGCGCTCTGCACCGTAGGCAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATTCTTGATCTGGCGAAACATTGCGAAttt                                >  2:2748899/1‑100 (MQ=255)
                                                                 ccGGCGGGCAATATGGCGCTCTGCACCGTGGGTAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATTCTTGATCTGGCGAAACATTGCGaa                                   >  2:1674487/1‑96 (MQ=255)
                                                                 ccGGCGGGCAATACGGCGCTCTGCACCGTGGGGAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATTCTTGATCTGGCGAAACATTGCGaa                                   >  2:3523723/1‑96 (MQ=255)
                                                                 ccGGCGGGCAATACGGCGCTCTGCACCGTGGGGAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATTCTTGATCTGGCGAAACATTGCGAATTTg                               >  2:2378038/1‑100 (MQ=255)
                                                                 ccGGCGGGCAATACGGCGCTCTGCACCGTGGGGAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATTCTTGATCTGGCGAAACATTGCGAATTTg                               >  2:3459938/1‑100 (MQ=255)
                                                                  nggcggGCAATACGGCGCTCTGCACCGTGGGTAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATTCTTGATCTGGCGAAAAATTGCGAAttt                                >  2:380139/2‑98 (MQ=255)
                                                                   ggcggGCAATACGGCGCTCTGCACCGTGGGTAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATTCTTGATCTGGCGAAACATTGCGAATTTgaa                             <  1:851551/100‑1 (MQ=255)
                                                                   ggcggGCAATACGGCGCTCTGCACCGTGGGTAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATTCTTGATCTGGCGAAACATTGCGAATTTga                              >  2:2667859/1‑99 (MQ=255)
                                                                   ggcggGCAATACGGCGCTCTGCACCGTGGGGAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATTCTTGATCTGGCGAAACATTGCGAATTTgaa                             <  1:3809633/100‑1 (MQ=255)
                                                                    gcggGCAATACGGCGCTCTGCACCGTGGGTAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATTCTTGATCTGGCGAAACATTGCGAATTTgaag                            <  1:3368287/100‑1 (MQ=255)
                                                                    gcggGCAATACGGCGCTCTGCACCGTGGGTAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATTCTTGATCTGGCGAAACATTGCGAATTTgaag                            <  1:3368211/100‑1 (MQ=255)
                                                                    gcggGCAATACGGCGCTCTGCACCGTGGGGAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATTCTTGATCTGGCGAAACATTGCGAATTTaaag                            >  2:2860591/1‑100 (MQ=255)
                                                                     cggGCAATACGGCGCTCTGCACCGTGGGGAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATTCTTGATCTGGCGAAACATTGCGAATTTgaaga                           <  1:142880/100‑1 (MQ=255)
                                                                     cggGCAATACGGCGCTCTGCACCGTGGGGAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATTCTTGACCTGGCGAAACATTGCGAATTTgaaga                           <  1:2538007/100‑1 (MQ=255)
                                                                       ggCAATACGGCGCTCTGCACCGTGGGTAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATTCTTGATCTGGCGAAACATTGCGAATTTGaagaag                         <  1:2632076/100‑1 (MQ=255)
                                                                       ggCAATACGGCGCTCTGCACCGTGGGGAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATTCTTGATCTGGCGAAACATTGCGAATTTga                              <  2:2502920/95‑1 (MQ=255)
                                                                        gCAATACGGCGCTCTGCACCGTGGGGAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATTCTTGATCTGGCGAAACATTGCGAATTTGAAGAAGt                        <  1:2489859/100‑1 (MQ=255)
                                                                        gCAATACGGCGCTCTGCACCGTGGGGAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATTCTTGATCTGGCGAAACATTGCGAATTTGAAGAAGt                        <  1:930479/100‑1 (MQ=255)
                                                                             aCGGCGCTCTGCACCGTGGGGAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATCTTCTTGATCTGGCGCAACATTGCGAATTTGAAGAAGTGgcg                    >  2:1809732/1‑99 (MQ=255)
                                                                                  gcTCTGCACCGTGGGGAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATTCTTGATCTGGCGAAACATTGCGAATTTGAAGAAGTGGCGCATctg              >  2:3629800/1‑100 (MQ=255)
                                                                                   ctctGCACCGTGGGGAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATTCTTGATCTGGCGAAACATTGCGAATTTGAAGAAGTGGCGCATctgc             <  1:2148222/100‑1 (MQ=255)
                                                                                    tctGCACCGTGGGTAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATTCTTGATCTGGCGAAACATTGCGAATTTGAAGAAGTGGCGCATctgct            <  1:3545710/100‑1 (MQ=255)
                                                                                      tGCACCGTGGGGAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATTCTTGATCTGGCGAAACATTGCGCATTTGAAGAAGTGGCGCATctactg           >  2:4020871/1‑99 (MQ=255)
                                                                                          ccGTGGGGAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATTCTTGATCTGGCGAAACATTGCGAATTTGCAGAAGTGACGCATCTGCTGATCCa      >  2:763356/1‑100 (MQ=255)
                                                                                            gTGGGGAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATTCTTGATCTGGCGAAACATTGCGAATTTGAAGAAGTGGCGCATCTGCTGATCCACg    <  1:2737673/100‑1 (MQ=255)
                                                                                             tGGGGAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATTCTTGATCTGGCGAAACATTGCGAATTTGAAGAAGTGGCGCATCTGCTGATCCACgg   <  1:3771277/100‑1 (MQ=255)
                                                                                             tGGGGAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATTCTTGATCTGGCGAAACATTGCGAATTTGAAGAAGTGGCGCATCTGCTGATCCACg    >  2:799960/1‑99 (MQ=255)
                                                                                              ggggAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATTCTTGATCTGGCGAAACATTGCGAGTTTGGAGAAGTGGCGCATATGCTGATCCACGGc  >  2:2082819/1‑100 (MQ=255)
                                                                                              ggggAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATTCTTGATCTGGCGAAACATTGCGAATTTGAAGAAGTGGCGCATCTGCTGATCCACGGc  <  1:2839358/100‑1 (MQ=255)
                                                                                                 |                                                                                                
CAACGATCCTGCAAAACAGTACCCATGTCATTAAACCGAAAAAATCTGTGGCACTTTCTGGCGTTCCGGCGGGCAATACGGCGCTCTGCACCGTGGGTAAAAGTGGCAATGACCTGCATTACCGCGGCTACGATATTCTTGATCTGGCGAAACATTGCGAATTTGAAGAAGTGGCGCATCTGCTGATCCACGGC  >  NC_000913/350022‑350215

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