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‑‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTGGCCTTCGGCG  >  NC_000913/3100894‑3101165
                                                                                                                                             |                                                                                                                                       
aaaTAGAGGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑g                                                                                                                                          <  1:368417/139‑1 (MQ=255)
 aaTAGAGGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑G‑‑c                                                                                                                                       <  2:1373798/139‑1 (MQ=255)
     gagGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑G‑‑ccccc                                                                                                                                  >  1:1218695/1‑139 (MQ=255)
      agGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑G‑‑cccccc                                                                                                                                 <  1:504523/139‑1 (MQ=255)
      agGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑G‑‑cccccc                                                                                                                                 >  1:1146111/1‑139 (MQ=255)
      agGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑G‑‑cccccc                                                                                                                                 >  2:1057653/1‑139 (MQ=255)
       ggATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑G‑‑ccccccc                                                                                                                                >  1:1237526/1‑139 (MQ=255)
           tattTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑G‑‑ccccccccagg                                                                                                                            <  1:818732/139‑5 (MQ=255)
           tattTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑G‑‑ccccccccagg                                                                                                                            <  2:1792298/139‑5 (MQ=255)
                    ttCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGc                                                                                                                    <  2:1415207/139‑1 (MQ=255)
                     tCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGca                                                                                                                   <  2:1237526/139‑1 (MQ=255)
                     tCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGca                                                                                                                   <  2:1220113/139‑1 (MQ=255)
                       ccGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTc                                                                                                                         <  2:1454127/131‑1 (MQ=255)
                       ccGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTc                                                                                                                         >  1:1454127/1‑131 (MQ=255)
                       ccGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGcaca                                                                                                                 >  1:1678870/1‑139 (MQ=255)
                       ccGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGc                                                                                                                    >  2:875229/1‑136 (MQ=255)
                       ccGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGc                                                                                                                    <  1:875229/136‑1 (MQ=255)
                        cGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTCGGATTGAGGATTGCACCC‑‑CCCCCCCAGGCTTCGCAGcaca                                                                                                                 <  2:1146111/139‑1 (MQ=255)
                         gACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCCAGGCTTCCCAGCACAt                                                                                                                <  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‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTAccg                                                                                               >  2:269006/1‑139 (MQ=255)
                                          tATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTAccgc                                                                                              <  1:1057653/139‑1 (MQ=255)
                                          tATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTAccgc                                                                                              <  2:1766975/139‑1 (MQ=255)
                                            tCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTAccgccg                                                                                            <  1:354650/139‑1 (MQ=255)
                                                   gagTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCa                                                                                     >  2:596581/1‑139 (MQ=255)
                                                      tCACTACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCa                                                                                  >  1:542554/1‑139 (MQ=255)
                                                      tCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACAt                                                                                                                >  1:48589/1‑109 (MQ=255)
                                                      tCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACAt                                                                                                                <  2:48589/109‑1 (MQ=255)
                                                      tCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑G‑‑cccccc                                                                                                                                 >  1:1265608/1‑91 (MQ=255)
                                                      tCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑G‑‑cccccc                                                                                                                                 <  2:1265608/91‑1 (MQ=255)
                                                        acaACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTAc                                                                                                 >  2:1384313/1‑122 (MQ=255)
                                                        acaACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTAc                                                                                                 <  1:1384313/122‑1 (MQ=255)
                                                        acaACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCAcc                                                                                >  1:1281261/1‑139 (MQ=255)
                                                        acaACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTCTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCAcc                                                                                >  2:936595/1‑139 (MQ=255)
                                                              cgGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTc                                                                          >  1:31658/1‑139 (MQ=255)
                                                                    cAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGccc                                                                    <  2:1382122/139‑1 (MQ=255)
                                                                      aaCATCTTTGCGAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCAt                                                                  <  2:1218695/139‑1 (MQ=255)
                                                                      aaCATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCAt                                                                  <  1:1602124/139‑1 (MQ=255)
                                                                      aaCATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCAt                                                                  >  2:1602124/1‑139 (MQ=255)
                                                                      aaCATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGAACCGCTTTCGCGCCCAt                                                                  <  1:269006/139‑1 (MQ=255)
                                                                              tGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTCCCACTCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGAGCCCATTTCATTaa                                                          >  2:1150285/1‑139 (MQ=255)
                                                                               gCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCCAGGCTTCCCAGCACAt                                                                                                                <  2:1120549/85‑1 (MQ=255)
                                                                               gCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAAt                                                         >  2:1696582/1‑139 (MQ=255)
                                                                               gCCAAACTGCTGGAAATCTTTTTCCTGATTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCCAGGCTTCCCAGCACAt                                                                                                                >  1:1120549/1‑85 (MQ=255)
                                                                                ccAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATg                                                        <  2:31658/139‑1 (MQ=255)
                                                                                ccAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATg                                                        <  2:1281261/139‑1 (MQ=255)
                                                                                   aaCTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTg                                                     >  2:591406/1‑139 (MQ=255)
                                                                                   aaCTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑G‑‑CCCCACCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTg                                                     >  1:1614296/1‑139 (MQ=255)
                                                                                        ctgGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTg                                                >  1:1700113/1‑139 (MQ=255)
                                                                                           gAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCCGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTaa                                             <  1:1614006/139‑1 (MQ=255)
                                                                                                cTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCt                                                      >  1:122358/1‑125 (MQ=255)
                                                                                                cTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCt                                                      <  2:122358/125‑1 (MQ=255)
                                                                                                cTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTc                                        <  2:1700113/139‑1 (MQ=255)
                                                                                                cTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTc                                        >  2:1563194/1‑139 (MQ=255)
                                                                                                cTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTc                                        <  2:1571402/139‑1 (MQ=255)
                                                                                                  ttttCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACAt                                                                                                                <  2:1192380/65‑1 (MQ=255)
                                                                                                  ttttCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACAt                                                                                                                >  1:1192380/1‑65 (MQ=255)
                                                                                                    ttCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTACAcc                                    <  2:1678870/139‑1 (MQ=255)
                                                                                                           tttttCCACTCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGcc                             <  1:596581/139‑1 (MQ=255)
                                                                                                              ttCCACTCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGa                          <  2:542554/139‑1 (MQ=255)
                                                                                                               tCCACTCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGAt                         <  1:591406/139‑1 (MQ=255)
                                                                                                                    tCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTc                                        >  2:502740/1‑119 (MQ=255)
                                                                                                                    tCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTc                                        <  1:1014561/119‑1 (MQ=255)
                                                                                                                    tCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTc                                        <  1:502740/119‑1 (MQ=255)
                                                                                                                    tCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTc                                        >  2:1014561/1‑119 (MQ=255)
                                                                                                                    tCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACt                    <  1:1696582/139‑1 (MQ=255)
                                                                                                                    tCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACt                    >  2:50249/1‑139 (MQ=255)
                                                                                                                    tCGGTTTGGATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACt                    >  2:1164664/1‑139 (MQ=255)
                                                                                                                           ggATTGAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCt             >  1:1412833/1‑139 (MQ=255)
                                                                                                                                gAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACAtt                                                                                                               >  1:336818/1‑36 (MQ=255)
                                                                                                                                gAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACAtt                                                                                                               <  2:336818/36‑1 (MQ=255)
                                                                                                                                gAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTGGCCt        >  2:658519/1‑139 (MQ=255)
                                                                                                                                gAGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTGGCCt        >  2:881124/1‑139 (MQ=255)
                                                                                                                                 aGGATTGCA‑‑GCCCCCCCCAGGCTTCCCCGCACATTCTGTAATGGGTTACCGCCGGTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCAGTGCCTGATTCACTa                   <  2:1142143/127‑1 (MQ=255)
                                                                                                                                 aGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTc                       <  2:360286/123‑1 (MQ=255)
                                                                                                                                 aGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTc                       >  1:360286/1‑123 (MQ=255)
                                                                                                                                 aGGATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTa                   >  1:1142143/1‑127 (MQ=255)
                                                                                                                                  ggATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCAt                                                                  >  2:1746998/1‑79 (MQ=255)
                                                                                                                                  ggATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCAt                                                                  <  1:1746998/79‑1 (MQ=255)
                                                                                                                                  ggATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTGGCCtt       <  1:894243/138‑1 (MQ=255)
                                                                                                                                  ggATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTGGCCtt       >  2:894243/1‑138 (MQ=255)
                                                                                                                                  ggATTGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTGGCCTTc      >  2:248543/1‑139 (MQ=255)
                                                                                                                                      tGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAg                                                                                    <  1:826640/57‑1 (MQ=255)
                                                                                                                                      tGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAg                                                                                    >  2:826640/1‑57 (MQ=255)
                                                                                                                                      tGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTGGCCTTCGgcg  <  1:1563194/139‑1 (MQ=255)
                                                                                                                                      tGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTGGCCTTCGgcg  >  2:1349302/1‑139 (MQ=255)
                                                                                                                                      tGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTGGCCTTCGgcg  <  1:1166681/139‑1 (MQ=255)
                                                                                                                                      tGCA‑‑GCCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTGGCCTTCGgcg  <  1:1164664/139‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                       
AAATAGAGGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCA‑‑G‑‑CCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTGGCCTTCGGCG  >  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: