Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 761,664 A→C 5.1% E48A (GAA→GCA)  sucB → dihydrolipoyltranssuccinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913761,6640AC 5.1% 421.6 / 11.9 117E48A (GAA→GCA) sucBdihydrolipoyltranssuccinase
Reads supporting (aligned to +/- strand):  ref base A (57/54);  new base C (3/3);  total (60/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-01

Only 100 of 125 total aligned reads displayed.

AGATATTCTGGTCCCTGACCTGCCTGAATCCGTAGCCGATGCCACCGTCGCAACCTGGCATAAAAAACCCGGCGACGCAGTCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACAACGGTAACGTCTCGTCAGATCCTTGGTCGCCTGCGTGAAGGCAACAGCGCCGGTAAAGAAACCAGCGCCAAATCTGAAGAGAA  >  NC_000913/761533‑761802
                                                                                                                                   |                                                                                                                                          
    atTCTGGTCCCTGACCTGCCTGAATCCGTAGCCGATGCCACCGTCGCAACCTGGCATAAAAAACCCGGCGACGCAGTCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAa                                                                                                                                           <  1:653067‑M2/129‑2 (MQ=255)
    atTCTGGTCCCTGACCTGCCTGAATCCGTAGCCGATGCCACCGTCGCAACCTGGCATAAAAAACCCGGCGACGCAGTCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAa                                                                                                                                           >  2:653067‑M2/1‑128 (MQ=255)
      tCTGGTCCCTGACCTGCCTGAATCCGTAGCCGATGCCACCGTCGCAACCTGGCATAAAAAACCCGGCGACGCAGTCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGc                                                                                                                      >  4:408567/1‑148 (MQ=255)
             ccTGACCTGCCTGAATCCGTAGCCGATGCCACCGTCGCAACCTGGCATAAAAAACCCGGCGACGCAGTCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTgg                                                                                                               <  3:408567/148‑1 (MQ=255)
                 aCCTGCCTGAATCCGTAGCCGATGCCACCGTCGCAACCTGGCATAAAAAACCCGGCGACGCAGTCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGa                                                                                                                                            >  1:201750/1‑115 (MQ=255)
                         gAATCCGTAGCCGATGCCACCGTCGCAACCTGGCATAAAAAACCCGGCGACGCAGTCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTgg                                                                                                               >  2:25055/1‑136 (MQ=255)
                         gAATCCGTAGCCGATGCCACCGTCGCAACCTGGCATAAAAAACCCGGCGACGCAGTCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTgg                                                                                                   >  3:926105/1‑148 (MQ=255)
                          aaTCCGTAGCCGATGCCACCGTCGCAACCTGGCATAAAAAACCCGGCGACGCAGTCGTACGTGATGAAGTGCTGGTAGAAATCGAAATTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGa                                                                                                  <  4:14646/148‑1 (MQ=255)
                                tAGCCGATGCCACCGTCGCAACCTGGCATAAAAAACCCGGCGACGCAGTCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTTCCagtaccact                                                                                                                                   >  1:279496‑M2/1‑100 (MQ=255)
                                tAGCCGATGCCACCGTCGCAACCTGGCATAAAAAACCCGGCGACGCAGTCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTTCCagtaccact                                                                                                                                   <  2:279496‑M2/109‑10 (MQ=255)
                                  gCCGATGCCACCGTCGCAACCTGGCATAAAAAACCCGGCGACGCAGTCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAg                                                                                          >  4:788884/1‑148 (MQ=255)
                                     gATGCCACCGTCGCAACCTGGCATAAAAAACCCGGCGACGCAGTCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTa                                                                                       >  1:608456/1‑148 (MQ=255)
                                        gCCACCGTCGCAACCTGGCATAAAAAACCCGGCGACGCAGTCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGaa                                                                                           <  3:114317/141‑1 (MQ=255)
                                        gCCACCGTCGCAACCTGGCATAAAAAACCCGGCGACGCAGTCGTACGTGATGAAGTGCTGATAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACaa                                                                                    <  3:788884/148‑1 (MQ=255)
                                         ccACCGTCGCAACCTGGCATAAAAAACCCGGCGACGCAGTCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGt                                                                                                        >  3:818546/1‑127 (MQ=255)
                                         ccACCGTCGCAACCTGGCATAAAAAACCCGGCGACGCAGTCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGt                                                                                                        <  4:818546/127‑1 (MQ=255)
                                            ccGTCGCAACCTGGCATAAAAAACCCGGCGACGCAGTCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACa                                                                                     <  2:380713/143‑1 (MQ=255)
                                            ccGTCGCAACCTGGCATAAAAAACCCGGCGACGCAGTCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACa                                                                                     >  1:380713/1‑143 (MQ=255)
                                            ccGTCGCAACCTGGCATAAAAAACCCGGCGACGCAGTCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACAACGGt                                                                                >  2:893302/1‑148 (MQ=255)
                                                    aCCTGGCATAAAAAACCCGGCGACGCAGTCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAAACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACAACGGTAACGtct                                                                         <  3:1022155/147‑1 (MQ=255)
                                                     ccTGGCATAAAAAACCCGGCGACGCAGTCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACAACGGTaa                                                                              >  4:142837/1‑141 (MQ=255)
                                                     ccTGGCATAAAAAACCCGGCGACGCAGTCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACAACGGTaa                                                                              <  3:142837/141‑1 (MQ=255)
                                                             aaaaaaCCCGGCGACGCAGTCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATca                                                                                                                               <  2:744010/84‑1 (MQ=255)
                                                             aaaaaaCCCGGCGACGCAGTCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATca                                                                                                                               >  1:744010/1‑84 (MQ=255)
                                                             aaaaaaCCCGGCGACGCAGTCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGGGGTTCTg                                                                                                    >  4:747397/1‑111 (MQ=255)
                                                             aaaaaaCCCGGCGACGCAGTCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTc                                                                                                      >  2:873773/1‑109 (MQ=255)
                                                             aaaaaaCCCGGCGACGCAGTCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTc                                                                                                      <  1:873773/109‑1 (MQ=255)
                                                             aaaaaaCCCGGCGACGCAGTCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGaa                                                                                           >  4:406980/1‑120 (MQ=255)
                                                             aaaaaaCCCGGCGACGCAGTCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACAACGGTAACGTCTCGTCAGATcc                                                               >  2:405329/1‑148 (MQ=255)
                                                             aaaaaaCCCGGCGACGCAGTCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACAACGGTAACGTCTCGTCAGATcc                                                               >  3:417762/1‑148 (MQ=255)
                                                             aaaaaaCCCGGCGACGCAGTCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACAACGGTAACGTCTCGTCAGATcc                                                               <  4:417762/148‑1 (MQ=255)
                                                              aaaaaCCCGGCGACGCAGTCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACAACGGTAACGTCTCGTCAGATCCt                                                              <  4:926105/148‑1 (MQ=255)
                                                               aaaaCCCGGCGACGCAGTCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAagta                                                                                                                                        >  2:113709‑M2/1‑69 (MQ=255)
                                                               aaaaCCCGGCGACGCAGTCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAagta                                                                                                                                        <  1:113709‑M2/73‑5 (MQ=255)
                                                                     cGGCGACGCAGTCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAg                                                                                                <  2:457831/107‑1 (MQ=255)
                                                                     cGGCGACGCAGTCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAg                                                                                                >  1:457831/1‑107 (MQ=255)
                                                                     cGGCGACGCAGTCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACAACGGTAACGTCTCGTCAGATcc                                                               <  1:947730/140‑1 (MQ=255)
                                                                     cGGCGACGCAGTCGTACGTGATGAAGGGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACAACGGTAACGTCTCGTCAGATcc                                                               >  2:947730/1‑140 (MQ=255)
                                                                        cGACGCAGTCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACAACGGTAACGTCTCGTCAGATCCTTGGTCGCCTg                                                    <  2:608456/148‑1 (MQ=255)
                                                                            gCAGTCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTTCCagtaccactttgtca                                                                                                                             <  3:965138‑M2/71‑16 (MQ=255)
                                                                            gCAGTCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTTCCagtaccactttgtca                                                                                                                             >  4:965138‑M2/1‑56 (MQ=255)
                                                                             cAGTCGTATGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAagta                                                                                                                                        >  3:874744‑M2/1‑55 (MQ=255)
                                                                             cAGTCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATc                                                                                                                                >  4:628151/1‑67 (MQ=255)
                                                                             cAGTCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATc                                                                                                                                <  3:628151/67‑1 (MQ=255)
                                                                             cAGTCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCt                                                                                                                 <  4:1069254/82‑1 (MQ=255)
                                                                             cAGTCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCt                                                                                                                 >  3:1069254/1‑82 (MQ=255)
                                                                             cAGTCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACAACGGTAACGTCTCGTCAGATCCTTGGTCGCCTGCGTGa                                               >  1:465334/1‑148 (MQ=255)
                                                                             cAGTCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACAACGGTAACGTCTCGTCAGATCCTTGGTCGCCTGCGTGa                                               <  4:1065453/148‑1 (MQ=255)
                                                                             cAGTCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACAACGGTAACGTCTCGTCAGATCCTTGGTCGCCTGCGTGa                                               >  1:870678/1‑148 (MQ=255)
                                                                                tCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTTCCagtaccactt                                                                                                                                  >  4:653312‑M2/1‑52 (MQ=255)
                                                                                tCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTTCCagtaccactt                                                                                                                                  <  3:653312‑M2/62‑11 (MQ=255)
                                                                                tCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAagta                                                                                                                                        >  4:3593‑M2/1‑52 (MQ=255)
                                                                                tCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAagta                                                                                                                                        <  3:3593‑M2/56‑5 (MQ=255)
                                                                                  gTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACAACGGTAACGTCTCGTCAGATCCTTGGTCGCCTGCGTGAAGGCa                                          <  2:870678/148‑1 (MQ=255)
                                                                                   tACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGa                                                                                                  >  1:168396/1‑91 (MQ=255)
                                                                                   tACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGa                                                                                                  <  2:168396/91‑1 (MQ=255)
                                                                                     cGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGt                                                                                        >  3:605608/1‑99 (MQ=255)
                                                                                     cGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGt                                                                                        <  4:605608/99‑1 (MQ=255)
                                                                                        gatgaAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACAACGGTAACGTCTCGTCAGATCCTTGGt                                                          <  4:1065277/126‑1 (MQ=255)
                                                                                        gatgaAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACAACGGTAACGTCTCGTCAGATCCTTGGt                                                          >  3:1065277/1‑126 (MQ=255)
                                                                                           gaAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACAACGGTAACGtct                                                                         >  3:209051/1‑108 (MQ=255)
                                                                                           gaAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACAACGGTAACGTCTCGTCAGATCCTTGGTCGCCTGCGTGAAGGCAACAGCGCCg                                 >  1:713901/1‑148 (MQ=255)
                                                                                            aaGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATg                                                                                                            >  3:699702/1‑72 (MQ=255)
                                                                                            aaGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAACTAACGGCATCAGCAGACGGCATTCTGGATg                                                                                                            <  4:699702/72‑1 (MQ=255)
                                                                                              gTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACAACGGTAACGTCTCGTCAGATCCtt                                                             <  4:190600/117‑1 (MQ=255)
                                                                                              gTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACAACGGTAACGTCTCGTCAGATCCtt                                                             >  3:190600/1‑117 (MQ=255)
                                                                                               tGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACAACGGTAACGTCTCGTCAGATCCTTGGTCGCCTGCGTGAAGGCAACAGCGCCGGTaa                             <  2:713901/148‑1 (MQ=255)
                                                                                                  tGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACAACGGTAACGTCTCGTCAGATCCTTGGTCGCCTGCGTGAAGGCAACAGCGCCGGTAAAGa                          >  2:996037/1‑148 (MQ=255)
                                                                                                    gTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACAACGGTAACGTCTCGTCAGATCCTTGGTCGCCTGCGTGAAGGCAACAGCGCCGGTAAAGaaa                        <  1:996037/148‑1 (MQ=255)
                                                                                                     tAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGaa                                                                                           <  4:923262/80‑1 (MQ=255)
                                                                                                     tAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACAACGGTAACGTCTCGTCAGATCCTTGGTCGCCTGCGTGAAGGCAACAGCGCCGGTAAAGAAAc                       >  3:79049/1‑148 (MQ=255)
                                                                                                      aGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTc                                                                                                      >  4:258235/1‑68 (MQ=255)
                                                                                                      aGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACAACGGTAACGTCTCGTCAGATCCTTGGTCGCCTGCGTGAAGGCAACAGCGCCGGTAAAGAAAcc                      <  2:660889/148‑1 (MQ=255)
                                                                                                      aGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACAACGGTAACGTCTCGTCAGATCCTTGGTCGCCTGCGTGAAGGCAACAGCGCCGGTAAAGAAAcc                      >  1:660889/1‑148 (MQ=255)
                                                                                                         aaTCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACAACGGTAACGTCTCGTCAGATCCTTGGTCGCCTGCGTGaa                                              <  3:1102375/121‑1 (MQ=255)
                                                                                                         aaTCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACAACGGTAACGTCTCGTCAGATCCTTGGTCGCCTGCGTGaa                                              <  3:567031/121‑1 (MQ=255)
                                                                                                         aaTCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGAGGAAGGTACAACGGTAACGTCTCGTCAGATCCTTGGTCGCCTGCGTGaa                                              >  4:1102375/1‑121 (MQ=255)
                                                                                                          aTCGCAACTGACAGAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACAACGGTAACGTCTCGTCAGATCCTTGGTCGCCTGCGTGaa                                              <  4:934820/120‑1 (MQ=255)
                                                                                                          aTCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACAACGGTAACGTCTCGTCAGATCCTTGGTCGCCTACGTGaa                                              >  3:934820/1‑120 (MQ=255)
                                                                                                            cGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATg                                                                                             >  2:308577/1‑71 (MQ=255)
                                                                                                            cGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATg                                                                                             <  1:308577/71‑1 (MQ=255)
                                                                                                            cGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACAACGGTAACGTCTCg                                                                       <  1:386824/93‑1 (MQ=255)
                                                                                                               aaCTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGa                                                                                                  >  2:540161/1‑63 (MQ=255)
                                                                                                               aaCTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGa                                                                                                  <  1:540161/63‑1 (MQ=255)
                                                                                                                aCTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGaa                                                                                           >  2:967870/1‑69 (MQ=255)
                                                                                                                 cTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACAACGGTAACGTCTCGTCAGATCCTTGGTCGCCTGCGTGAAGGCAACAGCGCCGGTaa                             >  4:174190/1‑130 (MQ=255)
                                                                                                                 cTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACAACGGTAACGTCTCGTCAGATCCTTGGTCGCCTGCGTGAAGGCAACAGCGCCGGTaa                             <  3:174190/130‑1 (MQ=255)
                                                                                                                     cAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACAACGGTAACGTCTCGTCAGATCCTTGGTCGCCTGCg                                                  >  4:397145/1‑105 (MQ=255)
                                                                                                                     cAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACAACGGTAACGTCTCGTCAGATCCTTGGTCGCCTGCg                                                  <  3:397145/105‑1 (MQ=255)
                                                                                                                       aaGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACAACGGTAACGTCTCGTCAGATCCTTGGTCGCCTGCGTGAAGGCAACAGCGCCGGTAAAGAAAc                       <  4:607156/130‑1 (MQ=255)
                                                                                                                          tGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACAACGGTAACGTCTCGTCAGATCCTTGGTCGCCTGCGTGAAGGCAACAGCGCCGGTAAAGAAACCAGCGCCAAATCTGAAGAGaa  <  1:953917/148‑1 (MQ=255)
                                                                                                                            gTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACAACGGTAACGTCTCGTCAGATCCTTGGTCGCCTGCGTGAAGGCAACAGCGCCGGTAAAGAAACCAgcg                  <  3:997584/130‑1 (MQ=255)
                                                                                                                             tACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACAACGGTAACGtct                                                                         >  3:677236/1‑74 (MQ=255)
                                                                                                                              aCTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACAACGGTAACGTCTCGTCAGATCCTTGGTCGCCTGCGTGAAGGCAACAGCGCCGGTAAAGAAACCAGCGCCAAATCTGa        <  4:261750/138‑1 (MQ=255)
                                                                                                                              aCTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACAACGGTAACGTCTCGTCAGATCCTTGGTCGCCTGCGTGAAGGCAACAGCGCCGGTAAAGAAACCAGCGCCAAATCTGa        >  3:261750/1‑138 (MQ=255)
                                                                                                                                tGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACAACGGTAACGTCTCGTCAGATCCTTGGTCGCCTGCg                                                  >  3:360573/1‑94 (MQ=255)
                                                                                                                                   aaGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACAACGGTAACGTCTCGTCAGATCCTTGGTCGCCTGCGTGAAg                                             >  2:18744/1‑96 (MQ=255)
                                                                                                                                   aaGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACAACGGTAACGTCTCGTCAGATCCTTGGTCGCCTGCGTGAAg                                             <  1:18744/96‑1 (MQ=255)
                                                                                                                                   aaGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACAACGGTAACGTCTCGTCAGATCCTTGGTCGCCTGCGTGAAGGCAACAGCGCCGGTAAAGAAACCAGCGCCAAATCTGa        >  4:249431/1‑133 (MQ=255)
                                                                                                                                   aaGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACAACGGTAACGTCTCGTCAGATCCTTGGTCGCCTGCGTGAAGGCAACAGCGCCGGTAAAGAAACCAGCGCCAAATCTGa        <  3:249431/133‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                          
AGATATTCTGGTCCCTGACCTGCCTGAATCCGTAGCCGATGCCACCGTCGCAACCTGGCATAAAAAACCCGGCGACGCAGTCGTACGTGATGAAGTGCTGGTAGAAATCGAAACTGACAAAGTGGTACTGGAAGTACCGGCATCAGCAGACGGCATTCTGGATGCGGTTCTGGAAGATGAAGGTACAACGGTAACGTCTCGTCAGATCCTTGGTCGCCTGCGTGAAGGCAACAGCGCCGGTAAAGAAACCAGCGCCAAATCTGAAGAGAA  >  NC_000913/761533‑761802

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