Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,710,703 G→A 17.6% intergenic (‑228/+96) eptB ← / ← yhjX KDO phosphoethanolamine transferase, Ca(2+)‑inducible/pyruvate‑inducible inner membrane protein, putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,710,7030GA17.6% 221.4 / 34.3 97intergenic (‑228/+96)eptB/yhjXKDO phosphoethanolamine transferase, Ca(2+)‑inducible/pyruvate‑inducible inner membrane protein, putative transporter
Reads supporting (aligned to +/- strand):  ref base G (35/44);  new base A (9/8);  total (44/53)
Fisher's exact test for biased strand distribution p-value = 5.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAA  >  NC_000913/3710605‑3710802
                                                                                                  |                                                                                                   
gAAAAAGTGTGGGGAAAGAGTGTGGCTGAGGGTAAAACTTGGAGAAAGGCGAGTTTGTAAAGAGAATTTAGCGTTACAAAAAGCAGTAAAAAGGGACAg                                                                                                     >  1:195679/1‑99 (MQ=255)
gAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGt                                                                                                    >  2:2001219/1‑100 (MQ=255)
     aGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAAtt                                                                                               <  2:2541296/100‑1 (MQ=255)
          ggggAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGccc                                                                                          <  2:2155842/100‑1 (MQ=255)
          ggggAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAATTAACTGGccc                                                                                          <  2:304061/100‑1 (MQ=255)
            ggAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCaa                                                                                        >  2:1768364/1‑100 (MQ=255)
            ggAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCaa                                                                                        >  2:352099/1‑100 (MQ=255)
               aaGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAATTAACTGGcccttttt                                                                                     <  2:1567534/100‑6 (MQ=255)
               aaGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAATTAACTGGcccttttt                                                                                     <  2:238788/100‑6 (MQ=255)
                agagTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAATTAACTGGcccttttta                                                                                    <  2:2538925/100‑7 (MQ=255)
                 gagTGTGCATGAAGCACAAATTTACAAAAATGCAAATTTTAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAATTAACTGGccctttttat                                                                                   <  1:1630472/100‑8 (MQ=255)
                   gtgtGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTAc                                                                                 >  2:2118160/1‑100 (MQ=255)
                   gtgtGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTAc                                                                                 <  2:219072/100‑1 (MQ=255)
                   gtgtGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAATTAACTGGccctttttactg                                                                                 <  2:589791/100‑10 (MQ=38)
                   gtgtGCATGAAGCACAAAAAAACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAATTAACTGGccctttttactg                                                                                 <  1:1636427/100‑10 (MQ=37)
                    tgtgCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAATTAACTGGccctttttactgc                                                                                <  2:82976/100‑11 (MQ=38)
                    tgtgCATGAAGCAAAAAATAACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACg                                                                                <  1:2491174/100‑1 (MQ=255)
                     gtgCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGa                                                                               >  1:2319014/1‑100 (MQ=255)
                      tgCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAg                                                                              >  1:385193/1‑100 (MQ=255)
                       gCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGt                                                                             <  2:573889/100‑1 (MQ=255)
                        cATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTc                                                                            <  1:532167/100‑1 (MQ=255)
                        cATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTc                                                                            >  2:2308303/1‑100 (MQ=255)
                        cATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTc                                                                            <  2:1281734/100‑1 (MQ=255)
                         aTGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCa                                                                           >  2:1781422/1‑100 (MQ=255)
                          tGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCAt                                                                          >  1:2059958/1‑100 (MQ=255)
                           gAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGAAGTAAAAAGGGCCAGTTAATTGGCCAAAATTTTACGAGTCAtt                                                                         >  1:461510/1‑100 (MQ=255)
                            aaGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTg                                                                        <  1:160591/100‑1 (MQ=255)
                            aaGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTg                                                                        <  2:1702548/100‑1 (MQ=255)
                              gCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGtt                                                                      <  2:879916/100‑1 (MQ=255)
                                acaAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTta                                                                    <  1:1941732/100‑1 (MQ=255)
                                 caAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTtat                                                                   >  1:1732153/1‑100 (MQ=255)
                                     tttACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACAt                                                               >  1:185867/1‑100 (MQ=255)
                                         cAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCg                                                           >  2:999939/1‑100 (MQ=255)
                                         cAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCg                                                           >  2:1783638/1‑100 (MQ=255)
                                         cAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCg                                                           >  1:1221401/1‑100 (MQ=255)
                                          aaaaaTGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGt                                                          <  2:2573463/100‑1 (MQ=255)
                                           aaaaTGTAAGTTTGAAAAGAAACTTAAGCGTTACAAAAATCAGTAAAAAGGTCCAGTTAGTTGGCCCAAATTGTACGAGTCATTGTTTATACATTACg                                                           <  1:2265883/98‑1 (MQ=255)
                                           aaaaTGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGtt                                                         >  2:253711/1‑100 (MQ=255)
                                           aaaaTGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGtt                                                         >  1:2067130/1‑100 (MQ=255)
                                            agaTGGAAGTTGAAAGAAAGACTTTAGCGTGACAAAAATCAGTAAAAATGGCCAGTTAATTGGCCCAATTTTTACGAGTCATTGTTTATACATTTCGTTg                                                        <  1:1183348/98‑1 (MQ=255)
                                            aaaTGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGttn                                                        <  2:1652785/100‑2 (MQ=255)
                                             aaTGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGa                                                       <  2:1608058/100‑1 (MQ=255)
                                               tGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCGAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATg                                                     >  2:392720/1‑100 (MQ=255)
                                                 cAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCg                                                   >  1:1629361/1‑100 (MQ=255)
                                                 cAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTAATTGTTTATACATTTCGTTGATGCg                                                   >  1:663350/1‑100 (MQ=255)
                                                 cAAGTTTGAAAAGAGACCTTAGCGTTACAAAAAGCAATAAAAAAGGCGAGGTAATTGGGCCAACTTATACGAGTCAATGTTTATACATTTAGCTGATGCg                                                   >  1:790956/1‑100 (MQ=255)
                                                  aaGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGt                                                  >  2:204967/1‑100 (MQ=255)
                                                    gTTTGAACAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGttt                                                <  1:1392473/100‑1 (MQ=255)
                                                    gTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGttt                                                <  2:1721973/100‑1 (MQ=255)
                                                     tttGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTg                                               <  2:325616/100‑1 (MQ=255)
                                                       tGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCGAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGAt                                             <  2:855959/100‑1 (MQ=255)
                                                        gAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATg                                            <  1:460146/100‑1 (MQ=255)
                                                        gAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATg                                            <  2:2141817/100‑1 (MQ=255)
                                                          aaaGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGca                                          >  2:1890814/1‑100 (MQ=255)
                                                            agagACTTTATCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATGGTCTATACATTTCGTTGATGCGTTTGATGcaca                                        <  1:1702713/100‑1 (MQ=255)
                                                            agagACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGcaca                                        >  2:1413631/1‑100 (MQ=255)
                                                            agagACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGcaca                                        <  2:1916886/100‑1 (MQ=255)
                                                               gaCTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACgg                                     >  2:482498/1‑100 (MQ=255)
                                                                aCTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGa                                    <  1:343874/100‑1 (MQ=255)
                                                                aCTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGa                                    <  1:482496/100‑1 (MQ=255)
                                                                  tttAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTGACGAGTCATTGTTTATACATGTCGTTGATGCGTTTGATGCACACGGAAg                                  <  2:1732857/100‑1 (MQ=255)
                                                                   ttAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGc                                 <  2:543662/100‑1 (MQ=255)
                                                                   ttAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGc                                 <  2:573533/100‑1 (MQ=255)
                                                                       cGTTACAAAAAGCAGTAAAAAGGGCCATTTAATTGGACCAAATTTTACGAGTCAGTGTTTATACATTTCGTTGATGCGTTTTATGCACACGGAAg                                  <  1:1447762/95‑1 (MQ=255)
                                                                        gTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAg                            <  2:1923688/100‑1 (MQ=255)
                                                                        gTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAg                            >  2:452074/1‑100 (MQ=255)
                                                                         ttACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGc                           >  1:730197/1‑100 (MQ=255)
                                                                         ttACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGc                           <  2:1011168/100‑1 (MQ=255)
                                                                          tACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGATGCTGAAGc                           >  1:646061/1‑99 (MQ=255)
                                                                           aCAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGccc                         >  1:548464/1‑100 (MQ=255)
                                                                           aCAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGccc                         >  1:2571985/1‑100 (MQ=255)
                                                                            cAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCa                        >  2:105513/1‑100 (MQ=255)
                                                                             aaaaaGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAg                       <  2:2367865/100‑1 (MQ=255)
                                                                             aaaaaGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAg                       <  1:766616/100‑1 (MQ=255)
                                                                              aaaaGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGt                      <  1:1138227/100‑1 (MQ=255)
                                                                                 aGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGc                   <  2:878526/100‑1 (MQ=255)
                                                                                 aGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATACACACGGAAGCTGAAGCCCAGTAGc                   <  1:60741/100‑1 (MQ=255)
                                                                                  gCAGTAAAAAGGGCCAATTAATTGGCCCAACTTTTACGAGTCATTGTTTATACATTTCGTTGATGCTTTTGCTGCACACGGAAGCTGAAGGCCAg                       >  1:2007124/1‑95 (MQ=255)
                                                                                   gtaaaatttgggCCAATTAACTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTc                 >  2:2011026/10‑100 (MQ=38)
                                                                                   gtaaaatttgggCCAATTAACTGGCCCAAATTTTACGAGTCATTGAAAATACATTTCGTTGATGCGTTTGATGCACACGGAAGATGAAGCCCAGTAGCTc                 >  1:465446/10‑100 (MQ=37)
                                                                                   gtaaaatttgggCCAATTAACTGGCACAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCCCACGGAAGCTGAAGCCCAGTAGCTc                 >  2:875337/10‑100 (MQ=255)
                                                                                    ngTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGCTGCGTTGGATGCACACGGAAGCTGAAGCCCAGTAGCTcg                >  1:1818925/2‑100 (MQ=255)
                                                                                     gTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTcgc               >  1:1371947/1‑100 (MQ=255)
                                                                                     aaaatttgggCCAATTAACTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTcgc               >  1:786042/8‑100 (MQ=255)
                                                                                      aaatttgggCCAATTAACTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTcgcg              >  2:1398824/7‑100 (MQ=255)
                                                                                       aatttgggCCAATTAACTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCgg             >  1:744088/6‑100 (MQ=255)
                                                                                        aaaaGGGCCAGTTTATTGGCCCAAATTTTACGAGCCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGc            <  1:608974/100‑1 (MQ=255)
                                                                                        aaaaGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGc            >  2:1936309/1‑100 (MQ=255)
                                                                                        aaaaGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGc            <  2:1132954/100‑1 (MQ=255)
                                                                                          ttgggCCAATTAACTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTg          >  1:422587/3‑100 (MQ=255)
                                                                                          ttgggCCAATTAACTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTg          >  2:133212/3‑100 (MQ=255)
                                                                                          aaGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTg          <  2:239920/100‑1 (MQ=255)
                                                                                          aaGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTg          <  1:738815/100‑1 (MQ=255)
                                                                                            gggCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAg        <  2:2067143/100‑1 (MQ=255)
                                                                                            gggCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAg        <  2:1134351/100‑1 (MQ=255)
                                                                                             ggCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGc       <  1:352097/100‑1 (MQ=255)
                                                                                             ggCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGc       >  2:1493792/1‑100 (MQ=255)
                                                                                             ggCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGc       <  2:1044048/100‑1 (MQ=255)
                                                                                                  gTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTaa  <  2:2440472/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
GAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAA  >  NC_000913/3710605‑3710802

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