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.C97.4% 262.0 / ‑1.1 76coding (592/720 nt)yggNDUF2884 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (35/39);  total (36/40)
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.94e-01

AAATAGAGGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑G‑‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTGGCCTTCGGCG  >  NC_000913/3100894‑3101165
                                                                                                                                             |                                                                                                                                       
aaaTAGAGGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑g                                                                                                                                          <  1:368417/139‑1 (MQ=255)
 aaTAGAGGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑G‑‑c                                                                                                                                       <  2:1373798/139‑1 (MQ=255)
     gagGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑G‑‑cccc‑c                                                                                                                                  >  1:1218695/1‑139 (MQ=255)
      agGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑G‑‑cccc‑cc                                                                                                                                 <  1:504523/139‑1 (MQ=255)
      agGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑G‑‑cccc‑cc                                                                                                                                 >  1:1146111/1‑139 (MQ=255)
      agGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑G‑‑cccc‑cc                                                                                                                                 >  2:1057653/1‑139 (MQ=255)
       ggATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑G‑‑cccc‑ccc                                                                                                                                >  1:1237526/1‑139 (MQ=255)
           tattTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑G‑‑cccc‑ccccagg                                                                                                                            <  1:818732/139‑5 (MQ=255)
           tattTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑G‑‑cccc‑ccccagg                                                                                                                            <  2:1792298/139‑5 (MQ=255)
                    ttCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGc                                                                                                                    <  2:1415207/139‑1 (MQ=255)
                     tCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGca                                                                                                                   <  2:1237526/139‑1 (MQ=255)
                     tCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGca                                                                                                                   <  2:1220113/139‑1 (MQ=255)
                       ccGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTc                                                                                                                         <  2:1454127/131‑1 (MQ=255)
                       ccGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTc                                                                                                                         >  1:1454127/1‑131 (MQ=255)
                       ccGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGcaca                                                                                                                 >  1:1678870/1‑139 (MQ=255)
                       ccGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGc                                                                                                                    >  2:875229/1‑136 (MQ=255)
                       ccGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGc                                                                                                                    <  1:875229/136‑1 (MQ=255)
                        cGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTCGGATTGAGGATTGCACCC‑‑CCCC‑CCCAGGCTTCGCAGcaca                                                                                                                 <  2:1146111/139‑1 (MQ=255)
                         gACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GCCCCCC‑CCCAGGCTTCCCAGCACAt                                                                                                                <  2:999599/139‑1 (MQ=255)
                                    tGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑G‑‑c                                                                                                                                       <  1:1626067/104‑1 (MQ=255)
                                    tGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑G‑‑c                                                                                                                                       >  2:1626067/1‑104 (MQ=255)
                                    tGCAGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑G‑‑c                                                                                                                                       <  1:953856/104‑1 (MQ=255)
                                    tGCAGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑G‑‑c                                                                                                                                       >  2:953856/1‑104 (MQ=255)
                                         cTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTAccg                                                                                               >  2:269006/1‑139 (MQ=255)
                                          tATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTAccgc                                                                                              <  1:1057653/139‑1 (MQ=255)
                                          tATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTAccgc                                                                                              <  2:1766975/139‑1 (MQ=255)
                                            tCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTAccgccg                                                                                            <  1:354650/139‑1 (MQ=255)
                                                   gagTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCa                                                                                     >  2:596581/1‑139 (MQ=255)
                                                      tCACTACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCa                                                                                  >  1:542554/1‑139 (MQ=255)
                                                      tCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACAt                                                                                                                >  1:48589/1‑109 (MQ=255)
                                                      tCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACAt                                                                                                                <  2:48589/109‑1 (MQ=255)
                                                      tCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑G‑‑cccc‑cc                                                                                                                                 >  1:1265608/1‑91 (MQ=255)
                                                      tCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑G‑‑cccc‑cc                                                                                                                                 <  2:1265608/91‑1 (MQ=255)
                                                        acaACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTAc                                                                                                 >  2:1384313/1‑122 (MQ=255)
                                                        acaACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTAc                                                                                                 <  1:1384313/122‑1 (MQ=255)
                                                        acaACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCAcc                                                                                >  1:1281261/1‑139 (MQ=255)
                                                        acaACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTCTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCAcc                                                                                >  2:936595/1‑139 (MQ=255)
                                                              cgGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTc                                                                          >  1:31658/1‑139 (MQ=255)
                                                                    cAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGccc                                                                    <  2:1382122/139‑1 (MQ=255)
                                                                      aaCATCTTTGCGAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCAt                                                                  <  2:1218695/139‑1 (MQ=255)
                                                                      aaCATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCAt                                                                  <  1:1602124/139‑1 (MQ=255)
                                                                      aaCATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCAt                                                                  >  2:1602124/1‑139 (MQ=255)
                                                                      aaCATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGAACCGCTTTCGCGCCCAt                                                                  <  1:269006/139‑1 (MQ=255)
                                                                              tGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTCCCACTCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGAGCCCATTTCATTaa                                                          >  2:1150285/1‑139 (MQ=255)
                                                                               gCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GCCCCCC‑CCCAGGCTTCCCAGCACAt                                                                                                                <  2:1120549/85‑1 (MQ=255)
                                                                               gCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAAt                                                         >  2:1696582/1‑139 (MQ=255)
                                                                               gCCAAACTGCTGGAAATCTTTTTCCTGATTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GCCCCCC‑CCCAGGCTTCCCAGCACAt                                                                                                                >  1:1120549/1‑85 (MQ=255)
                                                                                ccAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATg                                                        <  2:31658/139‑1 (MQ=255)
                                                                                ccAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATg                                                        <  2:1281261/139‑1 (MQ=255)
                                                                                   aaCTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTg                                                     >  2:591406/1‑139 (MQ=255)
                                                                                   aaCTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑G‑‑CCCCACCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTg                                                     >  1:1614296/1‑139 (MQ=255)
                                                                                        ctgGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTg                                                >  1:1700113/1‑139 (MQ=255)
                                                                                           gAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCCGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTaa                                             <  1:1614006/139‑1 (MQ=255)
                                                                                                cTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCt                                                      >  1:122358/1‑125 (MQ=255)
                                                                                                cTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCt                                                      <  2:122358/125‑1 (MQ=255)
                                                                                                cTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTc                                        <  2:1700113/139‑1 (MQ=255)
                                                                                                cTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTc                                        >  2:1563194/1‑139 (MQ=255)
                                                                                                cTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTc                                        <  2:1571402/139‑1 (MQ=255)
                                                                                                  ttttCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACAt                                                                                                                <  2:1192380/65‑1 (MQ=255)
                                                                                                  ttttCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACAt                                                                                                                >  1:1192380/1‑65 (MQ=255)
                                                                                                    ttCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTACAcc                                    <  2:1678870/139‑1 (MQ=255)
                                                                                                           tttttCCACTCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGcc                             <  1:596581/139‑1 (MQ=255)
                                                                                                              ttCCACTCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGa                          <  2:542554/139‑1 (MQ=255)
                                                                                                               tCCACTCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGAt                         <  1:591406/139‑1 (MQ=255)
                                                                                                                    tCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTc                                        >  2:502740/1‑119 (MQ=255)
                                                                                                                    tCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTc                                        <  1:1014561/119‑1 (MQ=255)
                                                                                                                    tCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTc                                        <  1:502740/119‑1 (MQ=255)
                                                                                                                    tCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTc                                        >  2:1014561/1‑119 (MQ=255)
                                                                                                                    tCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACt                    <  1:1696582/139‑1 (MQ=255)
                                                                                                                    tCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACt                    >  2:50249/1‑139 (MQ=255)
                                                                                                                    tCGGTTTGGATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACt                    >  2:1164664/1‑139 (MQ=255)
                                                                                                                           ggATTGAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCt             >  1:1412833/1‑139 (MQ=255)
                                                                                                                                gAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACAtt                                                                                                               >  1:336818/1‑36 (MQ=255)
                                                                                                                                gAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACAtt                                                                                                               <  2:336818/36‑1 (MQ=255)
                                                                                                                                gAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTGGCCt        >  2:658519/1‑139 (MQ=255)
                                                                                                                                gAGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTGGCCt        >  2:881124/1‑139 (MQ=255)
                                                                                                                                 aGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCCGCACATTCTGTAATGGGTTACCGCCGGTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCAGTGCCTGATTCACTa                   <  2:1142143/127‑1 (MQ=255)
                                                                                                                                 aGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTc                       <  2:360286/123‑1 (MQ=255)
                                                                                                                                 aGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTc                       >  1:360286/1‑123 (MQ=255)
                                                                                                                                 aGGATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTa                   >  1:1142143/1‑127 (MQ=255)
                                                                                                                                  ggATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCAt                                                                  >  2:1746998/1‑79 (MQ=255)
                                                                                                                                  ggATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCAt                                                                  <  1:1746998/79‑1 (MQ=255)
                                                                                                                                  ggATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTGGCCtt       <  1:894243/138‑1 (MQ=255)
                                                                                                                                  ggATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTGGCCtt       >  2:894243/1‑138 (MQ=255)
                                                                                                                                  ggATTGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTGGCCTTc      >  2:248543/1‑139 (MQ=255)
                                                                                                                                      tGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAg                                                                                    <  1:826640/57‑1 (MQ=255)
                                                                                                                                      tGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAg                                                                                    >  2:826640/1‑57 (MQ=255)
                                                                                                                                      tGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTGGCCTTCGgcg  <  1:1563194/139‑1 (MQ=255)
                                                                                                                                      tGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTGGCCTTCGgcg  >  2:1349302/1‑139 (MQ=255)
                                                                                                                                      tGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTGGCCTTCGgcg  <  1:1166681/139‑1 (MQ=255)
                                                                                                                                      tGCA‑‑GC‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTGGCCTTCGgcg  <  1:1164664/139‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                       
AAATAGAGGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑G‑‑CCCC‑CCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTGGCCTTCGGCG  >  NC_000913/3100894‑3101165

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