Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,572,698 G→A R177C (CGC→TGC)  pqqL ← putative peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,572,6980GA99.6% 957.4 / NA 277R177C (CGC→TGC) pqqLputative peptidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (128/148);  minor base T (0/1);  total (128/149)
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.99e-01

Only 100 of 284 total aligned reads displayed.

CCGTGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCGCCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCATTAC  >  NC_000913/1572608‑1572787
                                                                                          |                                                                                         
ccGTGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCAcc                                                                                         >  1:1962138/1‑93 (MQ=255)
  gTGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCAt                                                                                       <  1:5295402/93‑1 (MQ=255)
   tGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCAtt                                                                                      <  2:1276325/93‑1 (MQ=255)
    ggCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCAttt                                                                                     <  2:3554156/93‑1 (MQ=255)
        aCAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCa                                                                                 <  2:825181/93‑1 (MQ=255)
          aGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCCTc                                                                               >  2:622293/1‑93 (MQ=255)
           gTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATcc                                                                              >  1:4138834/1‑93 (MQ=255)
           gTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATcc                                                                              <  2:1520480/93‑1 (MQ=255)
           gTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCATGAGAGGTGCAGCATTTCGCATcc                                                                              >  2:2379684/1‑93 (MQ=255)
            tATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCt                                                                             <  2:5096664/93‑1 (MQ=255)
            tATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCt                                                                             >  1:4091143/1‑93 (MQ=255)
             aTCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTg                                                                            <  2:3373465/93‑1 (MQ=255)
              tCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGa                                                                           <  1:8072220/93‑1 (MQ=255)
              tCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGa                                                                           >  1:6352081/1‑93 (MQ=255)
              tCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGa                                                                           >  2:3661351/1‑93 (MQ=255)
              tCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGa                                                                           <  1:4292231/93‑1 (MQ=255)
               ccATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGAt                                                                          >  2:9679283/1‑93 (MQ=255)
               ccATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGAt                                                                          >  2:8444792/1‑93 (MQ=255)
                cATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATg                                                                         >  2:982412/1‑93 (MQ=255)
                cATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATg                                                                         <  1:3716382/93‑1 (MQ=255)
                  tCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATggg                                                                       >  2:853549/1‑93 (MQ=255)
                   cAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGc                                                                      <  2:8558828/93‑1 (MQ=255)
                   cAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCACTTCGCATCCTGATGGGc                                                                      >  1:8989909/1‑93 (MQ=255)
                    aGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCa                                                                     >  1:7381807/1‑93 (MQ=255)
                     ggCCGAGAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCAc                                                                    <  2:1465680/93‑1 (MQ=255)
                      gCCGATAGGTTCACGGTCTAAATTACGGGAATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGAGCACCATTTCGCATCCTGATGGGCACg                                                                   >  2:3501614/1‑93 (MQ=255)
                         gATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCa                                                                <  2:9305118/93‑1 (MQ=255)
                         gATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCa                                                                >  1:3970818/1‑93 (MQ=255)
                          aTAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCAt                                                               >  2:1633323/1‑93 (MQ=255)
                           tAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGCGGTGCACCATTTCGCATCCTGATGGGCACGCCAtt                                                              >  2:7562260/1‑93 (MQ=255)
                             ggTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACTCCATTcc                                                            <  1:6605823/93‑1 (MQ=255)
                             ggTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTcc                                                            >  1:2047137/1‑93 (MQ=255)
                             ggTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTcc                                                            >  2:6022669/1‑93 (MQ=255)
                                tCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCa                                                         <  1:9129297/93‑1 (MQ=255)
                                tCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGGGCCCCATTTCGCATCCTGATGGGCACGCCATTCCTcg                                                         >  1:2245331/1‑92 (MQ=255)
                                  aCGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGt                                                       <  2:9524127/93‑1 (MQ=255)
                                   cGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGGATCCTCATGGGTACGCCATTCCTCAGTa                                                      >  2:7675092/1‑93 (MQ=255)
                                   cGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTa                                                      <  2:9212161/93‑1 (MQ=255)
                                    cgTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTaa                                                     <  1:2766104/92‑1 (MQ=255)
                                     gTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAAt                                                    <  2:3059581/93‑1 (MQ=255)
                                     gTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAAt                                                    <  1:3659789/93‑1 (MQ=255)
                                      tCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGGACGCCATTCCTCAGTAAtt                                                   <  2:697706/93‑1 (MQ=255)
                                        tAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTAc                                                 <  1:5329600/93‑1 (MQ=255)
                                        tAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTAc                                                 >  1:5820674/1‑93 (MQ=255)
                                          aaTTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGc                                               <  1:8200401/93‑1 (MQ=255)
                                          aaTTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGc                                               <  2:4471931/93‑1 (MQ=255)
                                             tACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCAc                                            <  2:8688527/93‑1 (MQ=255)
                                             tACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCAc                                            >  2:5847366/1‑93 (MQ=255)
                                               cGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGt                                          >  2:7922080/1‑93 (MQ=255)
                                                 ggTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTc                                        >  1:2709320/1‑93 (MQ=255)
                                                 ggTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGAACGCCATTCCTCAGTAATTACGCCACGTTc                                        >  2:1191794/1‑93 (MQ=255)
                                                  gTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCa                                       <  1:2986742/93‑1 (MQ=255)
                                                   tATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCATCATTTCGCATCCTGATGGCCACGCCATTCCTCAGTAATTACGCCACGTTCAg                                      <  2:5666836/93‑1 (MQ=255)
                                                   tATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAg                                      <  2:250032/93‑1 (MQ=255)
                                                    aTTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGc                                     <  1:515020/93‑1 (MQ=255)
                                                    aTTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGc                                     >  2:1916911/1‑93 (MQ=255)
                                                     tttGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCg                                    >  1:7660074/1‑93 (MQ=255)
                                                      ttGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGt                                   >  1:2347887/1‑93 (MQ=255)
                                                       tGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTc                                  <  1:2655479/93‑1 (MQ=255)
                                                       tGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTc                                  >  1:7083018/1‑93 (MQ=255)
                                                       tGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTc                                  >  1:7240307/1‑93 (MQ=255)
                                                          cagcagGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTAc                               >  2:7062931/1‑93 (MQ=255)
                                                           agcagGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACt                              <  1:9670871/93‑1 (MQ=255)
                                                            gcagGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACtt                             <  1:10030741/93‑1 (MQ=255)
                                                               ggAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGa                          <  2:4090948/93‑1 (MQ=255)
                                                               ggAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGa                          >  2:5247945/1‑93 (MQ=255)
                                                               ggAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGa                          >  1:8365921/1‑93 (MQ=255)
                                                                gAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAg                         <  2:4052208/93‑1 (MQ=255)
                                                                  aaGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGtt                       >  1:3168419/1‑93 (MQ=255)
                                                                  aaGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGtt                       <  1:1578769/93‑1 (MQ=255)
                                                                   aGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGttt                      <  2:7412675/93‑1 (MQ=255)
                                                                   aGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGttt                      <  1:177701/93‑1 (MQ=255)
                                                                   aGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGttt                      >  1:5683622/1‑93 (MQ=255)
                                                                    gggCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGtttt                     <  2:4306548/93‑1 (MQ=255)
                                                                       cggcgCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCa                  <  1:1283440/93‑1 (MQ=255)
                                                                       cggcgCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCa                  <  1:9305568/93‑1 (MQ=255)
                                                                         gcgcgcCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCaaa                >  1:1661239/1‑93 (MQ=255)
                                                                         gcgcgcCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACACCACGTTCAGCGTCTACTTCGAGTTTTTCaaa                >  1:1179048/1‑93 (MQ=255)
                                                                          cgcgcCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAg               <  1:4225317/93‑1 (MQ=255)
                                                                          cgcgcCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAg               <  2:6065460/93‑1 (MQ=255)
                                                                           gcgcCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAgg              <  1:3448881/93‑1 (MQ=255)
                                                                             gcCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGtt            <  2:6164735/93‑1 (MQ=255)
                                                                             gcCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGtt            >  1:9743702/1‑93 (MQ=255)
                                                                             gcCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGtt            <  2:3112875/93‑1 (MQ=255)
                                                                             gcCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGtt            <  1:4591171/93‑1 (MQ=255)
                                                                              ccTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTg           >  1:3894791/1‑93 (MQ=255)
                                                                              ccTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTg           >  1:4784564/1‑93 (MQ=255)
                                                                               cTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGc          <  2:2260981/93‑1 (MQ=255)
                                                                                tGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCg         >  1:1872573/1‑93 (MQ=255)
                                                                                 gagaGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCgg        <  2:8117917/93‑1 (MQ=255)
                                                                                 gagaGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCgg        >  1:9385724/1‑93 (MQ=255)
                                                                                  agagGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGc       <  1:9424815/93‑1 (MQ=255)
                                                                                     ggTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCAtt    >  1:2537169/1‑93 (MQ=255)
                                                                                     ggTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCAtt    <  2:4919672/93‑1 (MQ=255)
                                                                                      gTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCATTa   <  1:4996144/93‑1 (MQ=255)
                                                                                      gTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCATTa   <  1:788351/93‑1 (MQ=255)
                                                                                       tGCACCATTTTGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCATTAc  >  2:6356796/1‑93 (MQ=255)
                                                                                       tGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCATTaa  >  2:2122859/1‑92 (MQ=255)
                                                                                       tGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCATTAc  <  1:4391575/93‑1 (MQ=255)
                                                                                       tGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCATTAc  >  2:6474034/1‑93 (MQ=255)
                                                                                          |                                                                                         
CCGTGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCGCCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCATTAC  >  NC_000913/1572608‑1572787

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