Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,101,039 (C)7→8 coding (585/720 nt) yggN ← DUF2884 family putative periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,101,0321.C96.7% 207.9 / ‑0.6 60coding (592/720 nt)yggNDUF2884 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base C (26/32);  total (26/34)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑G‑‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTGGCCTTCGGCG  >  NC_000913/3100900‑3101165
                                                                                                                                      |                                                                                                                                      
agGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑G‑‑cccccc                                                                                                                                 >  1:246718/1‑139 (MQ=255)
agGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑G‑‑cccccc                                                                                                                                 >  2:159294/1‑139 (MQ=255)
agGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑G‑‑cccccc                                                                                                                                 >  1:181704/1‑139 (MQ=255)
     tattTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑G‑‑ccccccccagg                                                                                                                            <  1:778465/139‑5 (MQ=255)
     tattTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑G‑‑ccccccccagg                                                                                                                            <  1:159294/139‑5 (MQ=255)
            aaTTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCa                                                                                                                      >  1:623191/1‑139 (MQ=255)
              ttCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGc                                                                                                                    <  2:623191/139‑1 (MQ=255)
              ttCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGc                                                                                                                    <  1:745599/139‑1 (MQ=255)
               tCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGca                                                                                                                   <  2:181704/139‑1 (MQ=255)
               tCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTCCCCAGca                                                                                                                   <  2:1152441/139‑1 (MQ=255)
                 ccGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGcaca                                                                                                                 >  1:1006219/1‑139 (MQ=255)
                   gACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACAtt                                                                                                               <  2:1054833/139‑1 (MQ=255)
                              tGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGt                                                                                                    >  1:611464/1‑139 (MQ=255)
                                      tCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTAccgccg                                                                                            <  2:634024/139‑1 (MQ=255)
                                      tCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACGGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGG‑‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTAccgccg                                                                                            <  2:735959/139‑1 (MQ=255)
                                                  acaACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCAc                                                                                 >  1:464420/1‑139 (MQ=255)
                                                  acaACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCAcc                                                                                >  2:1042330/1‑139 (MQ=255)
                                                  acaACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑G‑‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCg                                                                               <  2:391239/139‑1 (MQ=255)
                                                        cgGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGcac                                                                                                                  <  1:1170345/99‑1 (MQ=255)
                                                        cgGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGcac                                                                                                                  >  2:1170345/1‑99 (MQ=255)
                                                          gCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTc                                                                                                                         <  2:997179/90‑1 (MQ=255)
                                                          gCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTc                                                                                                                         >  1:997179/1‑90 (MQ=255)
                                                             aCAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGcc                                                                     <  1:1002473/139‑1 (MQ=255)
                                                              cAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGccc                                                                    >  1:1034765/1‑139 (MQ=255)
                                                                aaCATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCAt                                                                  <  2:1034765/139‑1 (MQ=255)
                                                                     cTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCAt                                                             >  2:608783/1‑139 (MQ=255)
                                                                        tGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGca                                                                                                                   >  1:1084230/1‑82 (MQ=255)
                                                                        tGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGca                                                                                                                   <  2:1084230/82‑1 (MQ=255)
                                                                        tGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑G‑‑c                                                                                                                                      >  2:347131/1‑62 (MQ=255)
                                                                        tGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑G‑‑c                                                                                                                                      >  2:388103/1‑62 (MQ=255)
                                                                        tGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑G‑‑c                                                                                                                                      <  1:347131/62‑1 (MQ=255)
                                                                        tGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑G‑‑c                                                                                                                                      <  1:388103/62‑1 (MQ=255)
                                                                         gCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAAt                                                         >  1:844288/1‑139 (MQ=255)
                                                                          ccAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGc                                                                                                                    <  2:997280/79‑1 (MQ=255)
                                                                          ccAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGc                                                                                                                    >  1:997280/1‑79 (MQ=255)
                                                                          ccAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATg                                                        <  2:468843/139‑1 (MQ=255)
                                                                          ccAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑G‑‑ccccccccagg                                                                                                                            >  1:1286883/1‑66 (MQ=255)
                                                                          ccAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑G‑‑ccccccccagg                                                                                                                            <  2:1286883/70‑5 (MQ=255)
                                                                             aaCTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTg                                                     >  2:436292/1‑139 (MQ=255)
                                                                              aCTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGt                                                    >  1:1182598/1‑139 (MQ=255)
                                                                                  ctgGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTAccgccg                                                                                            <  1:937863/95‑1 (MQ=255)
                                                                                  ctgGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTAccgccg                                                                                            >  2:937863/1‑95 (MQ=255)
                                                                                          cTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGTCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTc                                        <  2:246718/139‑1 (MQ=255)
                                                                                          cTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTc                                        <  2:611464/139‑1 (MQ=255)
                                                                                          cTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTc                                        <  1:1286921/139‑1 (MQ=255)
                                                                                          cTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTc                                        <  2:844288/139‑1 (MQ=255)
                                                                                            ttttCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTg                                                     >  2:462841/1‑124 (MQ=255)
                                                                                            ttttCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTg                                                     <  1:462841/124‑1 (MQ=255)
                                                                                            ttttCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCa                                      <  1:1042330/139‑1 (MQ=255)
                                                                                                  tGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCAtt                                <  2:1182598/139‑1 (MQ=255)
                                                                                                   gCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTa                                              >  1:19581/1‑125 (MQ=255)
                                                                                                   gCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GCCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTa                                              <  2:19581/125‑1 (MQ=255)
                                                                                                    cTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGca                                                                                                                   <  2:1201013/54‑1 (MQ=255)
                                                                                                    cTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGca                                                                                                                   >  1:1201013/1‑54 (MQ=255)
                                                                                                    cTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGCCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGc                              >  2:240901/1‑139 (MQ=255)
                                                                                                         tCCACTCGGTTTGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGAt                         <  2:1006219/139‑1 (MQ=255)
                                                                                                                  tttGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCa                                                                                                                      >  1:802505/1‑37 (MQ=255)
                                                                                                                  tttGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCa                                                                                                                      <  2:802505/37‑1 (MQ=255)
                                                                                                                  tttGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCa                                  <  2:890321/121‑1 (MQ=255)
                                                                                                                  tttGGATTGAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCa                                  >  1:890321/1‑121 (MQ=255)
                                                                                                                          gAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGAtt                        >  2:346060/1‑123 (MQ=255)
                                                                                                                          gAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGAtt                        <  1:346060/123‑1 (MQ=255)
                                                                                                                          gAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACt                    <  1:741780/127‑1 (MQ=255)
                                                                                                                          gAGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACt                    >  2:741780/1‑127 (MQ=255)
                                                                                                                           aGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCAcc                                                                                <  1:55468/66‑1 (MQ=255)
                                                                                                                           aGGATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCAcc                                                                                >  2:55468/1‑66 (MQ=255)
                                                                                                                            ggATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACt                    >  2:559573/1‑125 (MQ=255)
                                                                                                                            ggATTGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACt                    <  1:559573/125‑1 (MQ=255)
                                                                                                                                tGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCa                                      >  1:371594/1‑103 (MQ=255)
                                                                                                                                tGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCa                                      <  2:371594/103‑1 (MQ=255)
                                                                                                                                tGCA‑GC‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTGGCCTTCGgcg  <  1:509918/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                      
AGGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑G‑‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTGGCCTTCGGCG  >  NC_000913/3100900‑3101165

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