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.C98.0% 174.5 / ‑3.4 51coding (592/720 nt)yggNDUF2884 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (22/28);  total (22/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

AATAGAGGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑CC‑CCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTGGCCTT  >  NC_000913/3100895‑3101160
                                                                                                                                          |                                                                                                                                  
aaTAGAGGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑c                                                                                                                                  <  2:480962/139‑1 (MQ=255)
aaTAGAGGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑c                                                                                                                                  <  2:540534/139‑1 (MQ=255)
    gagGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑cc‑ccc                                                                                                                             >  2:147159/1‑139 (MQ=255)
    gagGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑cc‑ccc                                                                                                                             >  1:338512/1‑139 (MQ=255)
     agGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑cc‑cccc                                                                                                                            >  2:1479705/1‑139 (MQ=255)
     agGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑cc‑cccc                                                                                                                            >  1:1175913/1‑139 (MQ=255)
         ttattTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑cc‑ccccccca                                                                                                                        <  2:175842/139‑4 (MQ=255)
          tattTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑cc‑ccccccagg                                                                                                                       <  2:704234/139‑5 (MQ=255)
                 aaTTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGC‑CC‑CCCCCAGGCTTCCCa                                                                                                                 >  1:208387/1‑139 (MQ=255)
                  aTTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGC‑CC‑CCCCCAGGCTTCCCa                                                                                                                 <  1:430317/138‑1 (MQ=255)
                  aTTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGC‑CC‑CCCCCAGGCTTCCCa                                                                                                                 >  2:430317/1‑138 (MQ=255)
                    tCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑CCTCCCCCAGGCTTCCCAGca                                                                                                              <  1:1438387/139‑1 (MQ=255)
                     cccGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCC‑CCCCCAGGCTTCCCAGca                                                                                                              <  1:708463/139‑1 (MQ=255)
                     cccGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGC‑CC‑CCCCCAGGCTTCCCAGcac                                                                                                             <  1:1479705/139‑1 (MQ=255)
                      ccGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGCTTGAGGATTGCAGC‑CC‑CCCCCAGGCTTCCCAAcaca                                                                                                            >  2:1438387/1‑139 (MQ=255)
                      ccGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGC‑CC‑CCCCCAGGCTTCCCAGcaca                                                                                                            <  2:208387/139‑1 (MQ=255)
                       cGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGC‑CC‑CCCCCAGGCTTCCCAGCACAt                                                                                                           <  2:315248/139‑1 (MQ=255)
                              ggCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCTAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGC‑CC‑CCCCCAGGCTTCCCAGCACATTCTGTaa                                                                                                    >  2:1658153/1‑139 (MQ=255)
                              ggCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGC‑CC‑CCCCCAGGCTTCCCAGCACATTCTGTaa                                                                                                    >  1:1477762/1‑139 (MQ=255)
                                   tGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGC‑CC‑CCCCCAGGCTTCCCAGCACATTCTGTAATGGGt                                                                                               >  2:1619252/1‑139 (MQ=255)
                                          aTCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGC‑CC‑CCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTAccgcc                                                                                        >  2:648597/1‑139 (MQ=255)
                                           tCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGC‑CC‑CCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTAccgccg                                                                                       <  1:1093540/139‑1 (MQ=255)
                                             ttCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑cc‑ccccccagg                                                                                                                       <  2:73037/104‑5 (MQ=255)
                                             ttCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑cc‑ccccccagg                                                                                                                       >  1:73037/1‑100 (MQ=255)
                                                     tCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGAATTGCAGC‑CC‑CCCCCAGGCTTCCCAGca                                                                                                              <  1:1630944/106‑1 (MQ=255)
                                                     tCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGAATTGCAGC‑CC‑CCCCCAGGCTTCCCAGca                                                                                                              >  2:1630944/1‑106 (MQ=255)
                                                       acaACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑cc‑ccccccag                                                                                                                        >  1:1416326/1‑90 (MQ=255)
                                                       acaACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑cc‑ccccccag                                                                                                                        <  2:1416326/93‑4 (MQ=255)
                                                             cgGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTTGATTGAGGATTGCAGC‑CC‑CCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTc                                                                     <  2:1489260/139‑1 (MQ=255)
                                                               gCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGC‑CC‑CCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCg                                                                          <  2:726589/132‑1 (MQ=255)
                                                               gCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGC‑CC‑CCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCg                                                                          >  1:726589/1‑132 (MQ=255)
                                                                  aCAAACATCTTTGCTAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGC‑CC‑CCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGcc                                                                <  1:1658153/139‑1 (MQ=255)
                                                                     aaCATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGC‑CC‑CCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCAt                                                             <  2:1512724/139‑1 (MQ=255)
                                                                     aaCATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGC‑CC‑CCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCAt                                                             <  1:1065936/139‑1 (MQ=255)
                                                                     aaCATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGC‑CC‑CCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCAt                                                             <  2:895472/139‑1 (MQ=255)
                                                                     aaCATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGC‑CC‑CCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCAt                                                             <  1:1619252/139‑1 (MQ=255)
                                                                            ttGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGC‑CC‑CCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTcg                                                                    <  2:480462/125‑1 (MQ=255)
                                                                            ttGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGC‑CC‑CCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTcg                                                                    >  1:480462/1‑125 (MQ=255)
                                                                             tGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGC‑CC‑CCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTaa                                                     >  1:839111/1‑139 (MQ=255)
                                                                             tGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑cc‑ccccccagg                                                                                                                       <  2:1240122/72‑5 (MQ=255)
                                                                             tGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑cc‑ccccccagg                                                                                                                       >  1:1240122/1‑68 (MQ=255)
                                                                              gCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGC‑CC‑CCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAAt                                                    >  1:1412900/1‑139 (MQ=255)
                                                                               ccAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGC‑CC‑CCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATg                                                   <  2:461269/139‑1 (MQ=255)
                                                                               ccAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGC‑CC‑CCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATg                                                   <  2:1175913/139‑1 (MQ=255)
                                                                                       ctgGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGC‑CC‑CCCCCAGGCTTCCCAGc                                                                                                               >  2:855213/1‑71 (MQ=255)
                                                                                       ctgGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGC‑CC‑CCCCCAGGCTTCCCAGc                                                                                                               <  1:855213/71‑1 (MQ=255)
                                                                                               cTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGC‑CC‑CCCCCAGGCTTc                                                                                                                    >  2:756425/1‑58 (MQ=255)
                                                                                               cTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGC‑CC‑CCCCCAGGCTTc                                                                                                                    <  1:756425/58‑1 (MQ=255)
                                                                                               cTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGC‑CC‑CCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTcg                                                                    >  2:228599/1‑106 (MQ=255)
                                                                                               cTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGC‑CC‑CCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTcg                                                                    <  1:228599/106‑1 (MQ=255)
                                                                                               cTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGC‑CC‑CCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTc                                   <  2:1477762/139‑1 (MQ=255)
                                                                                               cTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGC‑CC‑CCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTc                                   >  1:715411/1‑139 (MQ=255)
                                                                                                 ttttCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGC‑CC‑CCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTg                                                >  2:655213/1‑124 (MQ=255)
                                                                                                 ttttCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGC‑CC‑CCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTg                                                <  1:655213/124‑1 (MQ=255)
                                                                                                 ttttCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑c                                                                                                                                  >  1:1024943/1‑42 (MQ=255)
                                                                                                 ttttCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑c                                                                                                                                  <  2:1024943/42‑1 (MQ=255)
                                                                                                        gCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGC‑CC‑CCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTg                          <  1:1650976/139‑1 (MQ=255)
                                                                                                             ttCCACTCGGTTTGGATTGAGGATTGCAGC‑CC‑CCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGa                     <  1:648597/139‑1 (MQ=255)
                                                                                                                               gAGGATTGCAGC‑CC‑CCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTGGCCt   >  1:1583001/1‑139 (MQ=255)
                                                                                                                               gAGGATTGCAGC‑CC‑CCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTGGCCt   >  1:602693/1‑139 (MQ=255)
                                                                                                                                aGGATTGCAGC‑CC‑CCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCttt                                                                      <  1:1115530/71‑1 (MQ=255)
                                                                                                                                aGGATTGCAGC‑CC‑CCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCttt                                                                      >  2:1115530/1‑71 (MQ=255)
                                                                                                                                aGGATTGCAGC‑CC‑CCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACt               <  1:1245766/126‑1 (MQ=255)
                                                                                                                                aGGATTGCAGC‑CC‑CCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACt               >  2:1245766/1‑126 (MQ=255)
                                                                                                                                aGGATTGCAGC‑CC‑CCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTGGCCtt  <  2:839111/139‑1 (MQ=255)
                                                                                                                                     tGCAGC‑CC‑CCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGAt                    >  1:1583915/1‑116 (MQ=255)
                                                                                                                                     tGCAGC‑CC‑CCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGAt                    <  2:1583915/116‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                  
AATAGAGGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑CC‑CCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTGGCCTT  >  NC_000913/3100895‑3101160

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