Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,149,963 G→T G341G (GGC→GGA yegI ← protein kinase‑related putative non‑specific DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,149,9630GT100.0% 230.8 / NA 72G341G (GGC→GGAyegIprotein kinase‑related putative non‑specific DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (31/41);  total (31/41)

AAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTGTTGGCGTAGAGAGTCCAGTGCC  >  NC_000913/2149830‑2150095
                                                                                                                                     |                                                                                                                                    
aaGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCgag                                                                                                                                 >  1:560927/1‑139 (MQ=255)
aaGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCgag                                                                                                                                 <  1:904353/139‑1 (MQ=255)
  ggAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGAt                                                                                                                               <  1:605856/139‑1 (MQ=255)
  ggAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGAt                                                                                                                               <  2:632903/139‑1 (MQ=255)
    aaGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCa                                                                                                                             <  2:862565/139‑1 (MQ=255)
       cGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATa                                                                                                                          <  2:884358/139‑1 (MQ=255)
       cGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATa                                                                                                                          <  2:1270332/139‑1 (MQ=255)
          ccTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATaaa                                                                                                                        <  1:1425906/138‑1 (MQ=255)
          ccTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATaaa                                                                                                                        >  2:1425906/1‑138 (MQ=255)
              ccAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATAAAATAGa                                                                                                                   >  1:48830/1‑139 (MQ=255)
               cAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATAAAATAGAt                                                                                                                  <  2:1520754/139‑1 (MQ=255)
               cAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATAAAATAGAt                                                                                                                  <  2:339328/139‑1 (MQ=255)
                 gTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATAAAATAGATaa                                                                                                                <  2:66762/139‑1 (MQ=255)
                    ggTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATAAAATAGATAacac                                                                                                             <  2:48830/139‑1 (MQ=255)
                        ggAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGTCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATAAAATAGATAACACCTTg                                                                                                         >  1:1308164/1‑139 (MQ=255)
                        ggAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATAAAATAGATAACACCTTg                                                                                                         >  2:27032/1‑139 (MQ=255)
                             tGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACTCCATCACCATCGCCCAGACTTTCGCCAGTACAAAATTCCCGCCGGTGGTAATGGCCTCTTCTCCGAGATCAATAAAATAGATAACACCTTGATTAt                                                                                                    <  2:560927/139‑1 (MQ=255)
                              gATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTATCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATAAAATAGATAACACCTTGATTATc                                                                                                   <  1:1067527/139‑1 (MQ=255)
                                tCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACCGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATAAAATAGATAACACCTTGATTATCCa                                                                                                 >  1:1588353/1‑139 (MQ=255)
                                    gTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATAAAATAGATAACACCTTGATTATCCAgcgc                                                                                             >  1:1055294/1‑139 (MQ=255)
                                         ggCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCAcacc                                                                                        >  2:820682/1‑139 (MQ=255)
                                              ttGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATAAAATAGAt                                                                                                                  <  2:583767/108‑1 (MQ=255)
                                              ttGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATAAAATAGAt                                                                                                                  >  1:583767/1‑108 (MQ=255)
                                              ttGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATAAAATAGATaa                                                                                                                <  1:1574763/110‑1 (MQ=255)
                                              ttGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATAAAATAGATaa                                                                                                                >  2:1574763/1‑110 (MQ=255)
                                              ttGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATAAAATAGATAACACCTTGATTATCCAgcg                                                                                              >  1:279219/1‑128 (MQ=255)
                                              ttGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATAAAATAGATAACACCTTGATTATCCAgcg                                                                                              <  2:279219/128‑1 (MQ=255)
                                              ttGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACggg                                                                                   >  2:389157/1‑139 (MQ=255)
                                              ttGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACggg                                                                                   >  1:146617/1‑139 (MQ=255)
                                               tGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGTCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGc                                                                                  <  2:1055294/139‑1 (MQ=255)
                                                 cAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGcag                                                                                <  1:389157/139‑1 (MQ=255)
                                                 cAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGcag                                                                                <  1:331677/139‑1 (MQ=255)
                                                 cAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGcag                                                                                <  1:1335952/139‑1 (MQ=255)
                                                  aaTGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGcagt                                                                               <  1:1483566/139‑1 (MQ=255)
                                                   aTGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGcagtc                                                                              >  1:1308622/1‑139 (MQ=255)
                                                      cTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATAAAATAGATAacac                                                                                                             >  2:285788/1‑105 (MQ=255)
                                                      cTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATAAAATAGATAacac                                                                                                             <  1:285788/105‑1 (MQ=255)
                                                        ggCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAgtcagtca                                                                         <  2:1190129/139‑1 (MQ=255)
                                                         gCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCaa                                                                        >  1:495303/1‑139 (MQ=255)
                                                              gCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAg                                                                   <  1:820682/139‑1 (MQ=255)
                                                              gCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAg                                                                   <  1:27032/139‑1 (MQ=255)
                                                               cTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGc                                                                  <  2:1521961/139‑1 (MQ=255)
                                                                  cTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCggg                                                               <  2:922455/139‑1 (MQ=255)
                                                                       gCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTaaa                                                          >  2:213525/1‑139 (MQ=255)
                                                                          aTCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATaaaa                                                                                                                       <  1:250872/75‑1 (MQ=255)
                                                                          aTCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATaaaa                                                                                                                       >  2:250872/1‑75 (MQ=255)
                                                                           tCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATg                                                      <  1:82724/139‑1 (MQ=255)
                                                                             aCCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCa                                                    <  2:1308622/139‑1 (MQ=255)
                                                                                 tCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGc                                                >  1:959083/1‑139 (MQ=255)
                                                                                    cccAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGa                                             >  2:698707/1‑139 (MQ=255)
                                                                                     ccAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGaa                                            >  2:916580/1‑139 (MQ=255)
                                                                                             tCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACa                                    >  2:155106/1‑139 (MQ=255)
                                                                                                   gCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATAAAATAGATAaca                                                                                                              <  2:1515957/59‑1 (MQ=255)
                                                                                                   gCACAAAATTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATAAAATAGATAaca                                                                                                              >  1:1515957/1‑59 (MQ=255)
                                                                                                         aaTTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTg                        <  2:1588353/139‑1 (MQ=255)
                                                                                                         aaTTCCCGCCGGTGGTAATGACCTCTTCTCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTg                        <  1:698653/139‑1 (MQ=255)
                                                                                                              ccgccgGTGGTAATGACCTCTTCTCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTGTTGGc                   >  2:32627/1‑139 (MQ=255)
                                                                                                                 ccgGTGGTAATGACCTCTTCTCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTg                        <  2:1563714/131‑1 (MQ=255)
                                                                                                                 ccgGTGGTAATGACCTCTTCTCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTg                        >  1:1563714/1‑131 (MQ=255)
                                                                                                                 ccgGTGGTAATGACCTCTTCTCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTGTTGGCGTa                <  2:146617/139‑1 (MQ=255)
                                                                                                                   ggtggtAATGACCTCTTCTCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTGTTGGCGTcga              <  2:959083/139‑1 (MQ=255)
                                                                                                                          aTGACCTCTTCTCCGAGATCAATAAAATAGATAACGCCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTGTTGGCGTAGAGAGTCCa       <  1:1000537/139‑1 (MQ=255)
                                                                                                                           tGACCTCTTCTCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTGTTGGCGTAGAGAGTCCAg      <  1:698707/139‑1 (MQ=255)
                                                                                                                               ctctTCTCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAg                                                                   >  2:1224932/1‑74 (MQ=255)
                                                                                                                               ctctTCTCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAg                                                                   <  1:1224932/74‑1 (MQ=255)
                                                                                                                               ctctTCTCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTg                        >  1:919833/1‑117 (MQ=255)
                                                                                                                               ctctTCTCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTg                        <  2:919833/117‑1 (MQ=255)
                                                                                                                               ctctTCTCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTGTTGGCGTAGAGAGTCCAGTGcc  >  2:1599817/1‑139 (MQ=255)
                                                                                                                               ctctTCTCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTGTTGGCGTAGAGAGTCCAGTGcc  >  2:1225693/1‑139 (MQ=255)
                                                                                                                               ctctTCTCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTGTTGGCGTAGAGAGTCCAGTGcc  >  1:968076/1‑139 (MQ=255)
                                                                                                                                 ctTCTCCGAGATCAATAAAATAGATAACACCTTGATTATc                                                                                                   <  2:573446/40‑1 (MQ=255)
                                                                                                                                 ctTCTCCGAGATCAATAAAATAGATAACACCTTGATTATc                                                                                                   >  1:573446/1‑40 (MQ=38)
                                                                                                                                     |                                                                                                                                    
AAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTGTTGGCGTAGAGAGTCCAGTGCC  >  NC_000913/2149830‑2150095

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