Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 974,875 C→T P2S (CCG→TCG)  mukE → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913974,8750CT99.6% 1073.4 / ‑5.1 275P2S (CCG→TCG) mukEhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (148/126);  minor base G (1/0);  total (149/126)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 278 total aligned reads displayed.

TCACCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAA  >  NC_000913/974785‑974964
                                                                                          |                                                                                         
  aCCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGGGCATGTCATCGACAAATATTGAACAAGTGATGTCGGt                                                                                       >  1:1621792/1‑93 (MQ=255)
  aCCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGt                                                                                       >  2:2731122/1‑93 (MQ=255)
  aCCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGt                                                                                       <  2:2216201/93‑1 (MQ=255)
    cGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTa                                                                                     <  1:3307155/93‑1 (MQ=255)
     ggACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTaa                                                                                    <  2:328219/93‑1 (MQ=255)
      gACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAg                                                                                   >  2:2944572/1‑93 (MQ=255)
      gACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAg                                                                                   >  2:589670/1‑93 (MQ=255)
      gACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAg                                                                                   >  1:2391924/1‑93 (MQ=255)
      gACTGCCAGCGAAATGGCAGCCGATTAATGATTAAGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAGGTTGTGGGGGTTAAg                                                                                   >  1:4956840/1‑93 (MQ=255)
         tGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTg                                                                                >  2:4783305/1‑93 (MQ=255)
          gCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTgg                                                                               >  1:3645088/1‑93 (MQ=255)
             aGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGgcgc                                                                            <  1:5789553/93‑1 (MQ=255)
                 aaaTGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGgc                                                                        >  2:737766/1‑93 (MQ=255)
                  aaTGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGgcg                                                                       >  1:647731/1‑93 (MQ=255)
                   aTGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGgcgc                                                                      <  2:3584257/93‑1 (MQ=255)
                   aTGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGgcgc                                                                      >  1:1332850/1‑93 (MQ=255)
                    tGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCt                                                                     >  1:4736842/1‑93 (MQ=255)
                     ggCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTg                                                                    <  2:499801/93‑1 (MQ=255)
                     ggCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTg                                                                    >  1:774998/1‑93 (MQ=255)
                     ggCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTg                                                                    >  2:1577279/1‑93 (MQ=255)
                      gCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTgg                                                                   >  2:228980/1‑93 (MQ=255)
                       cAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGc                                                                  <  2:4943012/93‑1 (MQ=255)
                        aGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCg                                                                 <  2:5849650/93‑1 (MQ=255)
                        aGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCg                                                                 >  1:4479801/1‑93 (MQ=255)
                            gATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATc                                                             <  1:6285341/93‑1 (MQ=255)
                             aTTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATcc                                                            >  1:2694874/1‑93 (MQ=255)
                              ttAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCg                                                           >  1:1672462/1‑93 (MQ=255)
                               tAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGt                                                          <  1:5684540/93‑1 (MQ=255)
                               tAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGt                                                          >  1:3832269/1‑93 (MQ=255)
                               tAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGt                                                          >  2:4686524/1‑93 (MQ=255)
                               tAATGATTACGGAGCCAAGGCACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTCAGCTGGCGCAGGCGCTAGCGAATCCGt                                                          >  1:1539026/1‑93 (MQ=255)
                                aaTGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGCATCCGtt                                                         >  2:962411/1‑93 (MQ=255)
                                aaTGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGtt                                                         <  1:2910498/93‑1 (MQ=255)
                                  tGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTAt                                                       <  2:1737268/93‑1 (MQ=255)
                                  tGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTAt                                                       <  2:5299279/93‑1 (MQ=255)
                                     ttACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTc                                                    <  2:3078750/93‑1 (MQ=255)
                                     ttACGGAGCCAAGATACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTc                                                    <  1:1847213/93‑1 (MQ=255)
                                       aCGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCg                                                  <  1:5273083/93‑1 (MQ=255)
                                         ggAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTCAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGgc                                                >  2:1337560/1‑93 (MQ=255)
                                          gAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGgcg                                               >  2:6091782/1‑93 (MQ=255)
                                          gAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGgcg                                               >  2:4168464/1‑93 (MQ=255)
                                          gAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGgcg                                               >  1:4865869/1‑93 (MQ=255)
                                           aGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGgcgc                                              <  2:5834365/93‑1 (MQ=255)
                                           aGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGgcgc                                              <  2:2850097/93‑1 (MQ=255)
                                            gCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCt                                             >  1:19157/1‑93 (MQ=255)
                                             ccAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTg                                            <  1:5552740/93‑1 (MQ=255)
                                              cAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTgg                                           >  2:2711613/1‑93 (MQ=255)
                                               aaGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGa                                          <  1:1158926/93‑1 (MQ=255)
                                               aaGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGa                                          <  1:1959336/93‑1 (MQ=255)
                                               aaGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGa                                          >  2:1396915/1‑93 (MQ=255)
                                                aGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGAc                                         >  2:5036197/1‑93 (MQ=255)
                                                aGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGAc                                         >  1:2507002/1‑93 (MQ=255)
                                                  gTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAg                                       <  1:2217065/93‑1 (MQ=255)
                                                  gTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAg                                       >  1:690205/1‑93 (MQ=255)
                                                   tACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAgc                                      >  2:6291761/1‑93 (MQ=255)
                                                   tACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAgc                                      <  2:5128631/93‑1 (MQ=255)
                                                     cAGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAgcgc                                    >  1:520717/1‑93 (MQ=255)
                                                      aGGCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGcc                                   >  2:17985/1‑93 (MQ=255)
                                                       ggCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCGGGCGCTGGCCAGCGCCt                                  >  1:2775633/1‑93 (MQ=255)
                                                       ggCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCt                                  <  1:1061787/93‑1 (MQ=255)
                                                       ggCGCATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCt                                  >  2:45831/1‑93 (MQ=255)
                                                          gcATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTAc                               >  2:1068423/1‑93 (MQ=255)
                                                          gcATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTAc                               >  2:449432/1‑93 (MQ=255)
                                                          gcATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTAc                               >  2:3819845/1‑93 (MQ=255)
                                                           cATGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACg                              >  2:5960013/1‑93 (MQ=255)
                                                            aTGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGt                             >  1:4033936/1‑93 (MQ=255)
                                                            aTGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGt                             >  1:6070657/1‑93 (MQ=255)
                                                             tGTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGtt                            <  1:4546776/93‑1 (MQ=255)
                                                              gTCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTc                           >  1:3966611/1‑93 (MQ=255)
                                                               tCATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCa                          >  1:3573160/1‑93 (MQ=255)
                                                                cATCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAg                         >  1:6093538/1‑93 (MQ=255)
                                                                 aTCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAgg                        >  2:6387408/1‑93 (MQ=255)
                                                                 aTCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAgg                        <  2:4245178/93‑1 (MQ=255)
                                                                 aTCGACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAgg                        <  1:2085251/93‑1 (MQ=255)
                                                                    gACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACg                     >  1:562392/1‑93 (MQ=255)
                                                                    gACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACg                     >  2:5087938/1‑93 (MQ=255)
                                                                    gACAAATATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACg                     >  1:4223017/1‑93 (MQ=255)
                                                                       aaaTATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCa                  <  2:847240/93‑1 (MQ=255)
                                                                       aaaTATTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCa                  <  2:5406322/93‑1 (MQ=255)
                                                                         atatTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCata                <  2:3127338/93‑1 (MQ=255)
                                                                          tatTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCatat               <  1:3045247/93‑1 (MQ=255)
                                                                          tatTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCatat               >  2:3524532/1‑93 (MQ=255)
                                                                           atTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATAtt              >  2:3905875/1‑93 (MQ=255)
                                                                           atTGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATAtt              <  2:6184220/93‑1 (MQ=255)
                                                                             tGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTgg            <  1:2278456/93‑1 (MQ=255)
                                                                             tGAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATACGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTgg            <  2:263010/93‑1 (MQ=255)
                                                                              gAACAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGc           >  2:2001657/1‑93 (MQ=255)
                                                                                aCAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCt         <  1:2025554/93‑1 (MQ=255)
                                                                                aCAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCt         <  1:928382/93‑1 (MQ=255)
                                                                                aCAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCt         >  1:5866182/1‑93 (MQ=255)
                                                                                aCAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCt         >  2:247981/1‑93 (MQ=255)
                                                                                 cAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTc        <  2:5771875/93‑1 (MQ=255)
                                                                                 cAAGTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTc        <  2:897971/93‑1 (MQ=255)
                                                                                    gTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTcgac     >  2:5217266/1‑93 (MQ=255)
                                                                                    gTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTcgac     >  1:5509801/1‑93 (MQ=255)
                                                                                    gTGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTcgac     <  2:5991304/93‑1 (MQ=255)
                                                                                     tGATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTcgacg    <  1:4974344/93‑1 (MQ=255)
                                                                                      gATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTcgacga   <  1:5392474/93‑1 (MQ=255)
                                                                                      gATGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTcgacga   <  1:1764513/93‑1 (MQ=255)
                                                                                       aTGTCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGaa  <  2:2571598/93‑1 (MQ=255)
                                                                                          |                                                                                         
TCACCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAA  >  NC_000913/974785‑974964

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