Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,698,281 A→G V300V (GTT→GTC malI ← transcriptional repressor of Mal regulon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,698,2810AG98.8% 301.1 / ‑4.5 85V300V (GTT→GTCmalItranscriptional repressor of Mal regulon
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (36/48);  total (37/48)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

AAGATTGCGTGAATGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCT  >  NC_000913/1698185‑1698376
                                                                                                |                                                                                               
aaGATTGCGTGAATGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACa                                                                                              <  2:981426/100‑1 (MQ=255)
 aGATTGCGTGAATGCGTCTCTTCATGGGTGATTTTTTGCATTATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAg                                                                                             >  2:845164/1‑100 (MQ=255)
 aGATTGCGTGAATGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAg                                                                                             <  2:1991660/100‑1 (MQ=255)
  gATTGCGTGAATGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAgg                                                                                            <  2:2326729/100‑1 (MQ=255)
  gATTGCGTGAATGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAgg                                                                                            >  1:1406464/1‑100 (MQ=255)
    ttGCGTGAATGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGta                                                                                          <  1:1331807/100‑1 (MQ=255)
    ttGCGTGAATGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGta                                                                                          <  1:909025/100‑1 (MQ=255)
    ttGCGTGAATGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGta                                                                                          <  1:828773/100‑1 (MQ=255)
      gCGTGAATGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGtata                                                                                        >  1:141238/1‑100 (MQ=255)
           aaTGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATatcatc                                                                                   <  2:696868/100‑1 (MQ=255)
            aTGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATAtcatca                                                                                  <  1:647528/100‑1 (MQ=255)
             tGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCaa                                                                                 <  2:1961543/100‑1 (MQ=255)
              gCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAg                                                                                <  2:6961/100‑1 (MQ=255)
               cGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAgt                                                                               >  1:948465/1‑100 (MQ=255)
                gTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAgtg                                                                              <  2:32427/100‑1 (MQ=255)
                gTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAgtg                                                                              <  1:845165/100‑1 (MQ=255)
                 tctcTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAgtgt                                                                             >  1:2542989/1‑100 (MQ=255)
                    ctTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGtggt                                                                          <  2:2569106/100‑1 (MQ=255)
                     ttCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGtggtt                                                                         >  2:1635419/1‑100 (MQ=255)
                      tCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGtggttg                                                                        <  1:2847797/100‑1 (MQ=255)
                       cATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGtggttgg                                                                       >  2:985353/1‑100 (MQ=255)
                        aTGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGtggttggt                                                                      >  1:1434558/1‑100 (MQ=255)
                         tGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGgtg                                                                     <  2:802096/100‑1 (MQ=255)
                          gggTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCGTCAAGTGTGGTTGgtgt                                                                    >  2:51189/1‑100 (MQ=255)
                           ggTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTc                                                                   <  1:1035553/100‑1 (MQ=255)
                           ggTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTc                                                                   <  2:1662332/100‑1 (MQ=255)
                           ggTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTc                                                                   >  2:193939/1‑100 (MQ=255)
                           ggTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTc                                                                   >  1:2781214/1‑100 (MQ=255)
                           ggTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTc                                                                   <  2:1434552/100‑1 (MQ=255)
                             tGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGc                                                                 <  2:1800977/100‑1 (MQ=255)
                             tGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGc                                                                 <  2:1750934/100‑1 (MQ=255)
                              gATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCa                                                                <  2:2638166/100‑1 (MQ=255)
                              gATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCa                                                                <  1:2853172/100‑1 (MQ=255)
                              gATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCa                                                                >  1:252419/1‑100 (MQ=255)
                                 tttttGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCg                                                             >  1:176741/1‑100 (MQ=255)
                                   tttGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGt                                                           <  1:193940/100‑1 (MQ=255)
                                   tttGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGt                                                           <  1:804454/100‑1 (MQ=255)
                                      gCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGCTATATCATCAAGTGTGGTTGGTGTCGCATCGGTaaa                                                        >  1:2029615/1‑100 (MQ=255)
                                        atcatGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATg                                                      <  2:1133523/100‑1 (MQ=255)
                                        atcatGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATg                                                      >  2:1232391/1‑100 (MQ=255)
                                         tcatGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATgc                                                     >  1:51707/1‑100 (MQ=255)
                                         tcatGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATgc                                                     >  2:2178980/1‑100 (MQ=255)
                                         tcatGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATgc                                                     >  2:1420501/1‑100 (MQ=255)
                                           atGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATgccn                                                   <  2:1225109/100‑2 (MQ=255)
                                            tGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATgccgc                                                  >  1:1341773/1‑100 (MQ=255)
                                            tGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATgccgc                                                  <  2:948464/100‑1 (MQ=255)
                                             gCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATgccgcc                                                 >  2:1048256/1‑100 (MQ=255)
                                                aTCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGc                                              >  1:2625419/1‑100 (MQ=255)
                                                      aaGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAAcc                                        <  1:172403/100‑1 (MQ=255)
                                                      aaGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAAcc                                        <  2:1137874/100‑1 (MQ=255)
                                                          gtgATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGtt                                    <  2:141238/100‑1 (MQ=255)
                                                          gtgATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGtt                                    <  1:2543812/100‑1 (MQ=255)
                                                          gtgATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGtt                                    <  1:2543943/100‑1 (MQ=255)
                                                          gtgATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGtt                                    >  1:2768651/1‑100 (MQ=255)
                                                           tcttccgatcttCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTg                                   <  1:1173290/90‑1 (MQ=255)
                                                             aTACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCt                                 <  2:2356107/100‑1 (MQ=255)
                                                              tACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTc                                <  1:2433555/100‑1 (MQ=255)
                                                              tACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTc                                >  1:1699766/1‑100 (MQ=255)
                                                               aCCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCa                               <  1:2602322/100‑1 (MQ=255)
                                                               aCCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCa                               >  2:1679094/1‑100 (MQ=255)
                                                                ccAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCaa                              >  1:2662957/1‑100 (MQ=255)
                                                                 cAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCaaa                             >  2:558413/1‑100 (MQ=255)
                                                                   aGTTCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGt                           <  1:1150817/100‑1 (MQ=255)
                                                                      tCCCGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAAc                        <  1:608596/100‑1 (MQ=255)
                                                                       cccGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACg                       <  1:1635431/100‑1 (MQ=255)
                                                                       cccGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACg                       >  1:2355829/1‑100 (MQ=255)
                                                                        ccGCGCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACg                       <  2:2235570/99‑1 (MQ=255)
                                                                         cgcgCTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGAt                     <  2:2781193/100‑1 (MQ=255)
                                                                             cTGGCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGAc                 >  2:1929357/1‑100 (MQ=255)
                                                                               ggCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTc               >  1:199846/1‑100 (MQ=255)
                                                                                gCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTcc              >  1:1168722/1‑100 (MQ=255)
                                                                                gCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTcc              <  2:2748956/100‑1 (MQ=255)
                                                                                gCGTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTcc              <  1:260922/100‑1 (MQ=255)
                                                                                  gTGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGc            <  1:1392962/100‑1 (MQ=255)
                                                                                   tGCTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCt           >  2:1806269/1‑100 (MQ=255)
                                                                                     cTGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCtt          <  2:1110778/99‑1 (MQ=255)
                                                                                      tGGCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTc        >  2:566801/1‑100 (MQ=255)
                                                                                       ggCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTcc       >  2:1669417/1‑100 (MQ=255)
                                                                                       ggCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTcc       >  2:1906639/1‑100 (MQ=255)
                                                                                       ggCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTcc       >  2:1466943/1‑100 (MQ=255)
                                                                                        gCCCAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTccc      <  1:1232395/100‑1 (MQ=255)
                                                                                         cccAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTcccc     >  2:816549/1‑100 (MQ=255)
                                                                                          ccAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCg    <  1:2630981/100‑1 (MQ=255)
                                                                                          ccAGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCg    >  1:1985186/1‑100 (MQ=255)
                                                                                            aGGTGACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCCGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCt  <  1:51189/100‑1 (MQ=255)
                                                                                                |                                                                                               
AAGATTGCGTGAATGCGTCTCTTCATGGGTGATTTTTTGCATCATGCGATCCGCAAGTGTGATACCAAGTTCCCGCGCTGGCGTGCTGGCCCAGGTAACAGGTATATCATCAAGTGTGGTTGGTGTCGCATCGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCT  >  NC_000913/1698185‑1698376

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