Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,238,562 (T)7→6 intergenic (+265/‑18) ygjR → / → alx putative NAD(P)‑dependent dehydrogenase/putative membrane‑bound redox modulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,238,5560T.98.5% 275.6 / ‑2.6 68intergenic (+259/‑24)ygjR/alxputative NAD(P)‑dependent dehydrogenase/putative membrane‑bound redox modulator
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (32/35);  total (32/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCA  >  NC_000913/3238418‑3238685
                                                                                                                                          |                                                                                                                                 
tgtgAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGt                                                                                                                                   >  1:1614473/1‑139 (MQ=255)
 gtgAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGtt                                                                                                                                  >  1:555037/1‑139 (MQ=255)
 gtgAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGtt                                                                                                                                  >  1:597701/1‑139 (MQ=255)
  tgAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGttt                                                                                                                                 >  1:1538370/1‑139 (MQ=255)
     aaaaaCACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGtttttt                                                                                                                              <  2:1538370/139‑1 (MQ=255)
               ccGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAgg                                                                                                                   <  2:555037/139‑1 (MQ=255)
                 ggTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGaa                                                                                                                 <  1:1234579/139‑1 (MQ=255)
                 ggTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGaa                                                                                                                 <  1:754907/139‑1 (MQ=255)
                               ccGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGAACTTCTATGAATACt                                                                                                   <  1:1122510/139‑1 (MQ=255)
                                   aaGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGAACt                                                                                                               <  2:1181311/123‑1 (MQ=255)
                                   aaGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTCGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGAACt                                                                                                               >  1:1181311/1‑123 (MQ=255)
                                        aTACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGaa                                                                                                                 <  1:258758/116‑1 (MQ=255)
                                        aTACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGaa                                                                                                                 >  2:258758/1‑116 (MQ=255)
                                         tACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGAACt                                                                                                               <  2:406368/117‑1 (MQ=255)
                                         tACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGAACt                                                                                                               >  1:406368/1‑117 (MQ=255)
                                              aGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGAACt                                                                                                               <  2:977889/112‑1 (MQ=255)
                                              aGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGAACt                                                                                                               >  1:977889/1‑112 (MQ=255)
                                                aCGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGaa                                                                                                                 <  2:1557704/108‑1 (MQ=255)
                                                aCGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGaa                                                                                                                 >  1:1557704/1‑108 (MQ=255)
                                                aCGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCt                                                                                  >  2:1424575/1‑139 (MQ=255)
                                                     ttCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGAACTTCt                                                                                                            >  2:195321/1‑108 (MQ=255)
                                                     ttCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGAACTTCt                                                                                                            <  1:195321/108‑1 (MQ=255)
                                                          tttttGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGttttttatg                                                                                                                           >  1:1359833/1‑86 (MQ=255)
                                                          tttttGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGttttttatg                                                                                                                           <  2:1359833/89‑4 (MQ=255)
                                                          tttttGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGttttt                                                                                                                               <  2:959010/85‑1 (MQ=255)
                                                          tttttGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGttttt                                                                                                                               >  1:959010/1‑85 (MQ=255)
                                                          tttttGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGAACTTCTATGa                                                                                                        <  1:409966/107‑1 (MQ=255)
                                                          tttttGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGAACTTCTATGa                                                                                                        >  2:409966/1‑107 (MQ=255)
                                                               gttgtAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGc                                                                   <  1:922720/139‑1 (MQ=255)
                                                                ttgtAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCt                                                                  >  2:1093190/1‑139 (MQ=255)
                                                                    aaGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGAACt                                                                                                               >  1:517574/1‑90 (MQ=255)
                                                                    aaGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGAACt                                                                                                               <  2:517574/90‑1 (MQ=255)
                                                                        gagaCCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCtgttgttgt                                                          <  2:672484/139‑1 (MQ=255)
                                                                          gaCCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCtgttgtt                                                            >  1:985763/1‑135 (MQ=255)
                                                                          gaCCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCtgttgtt                                                            <  2:985763/135‑1 (MQ=255)
                                                                           aCCTTGCCGGAAGGCGAGGACTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCTGTTGTTGTCGc                                                       >  2:1550135/1‑139 (MQ=255)
                                                                               tGCCGGACGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCG‑TTCTTTATGTGTAAGGAACt                                                                                                               <  1:957878/79‑1 (MQ=255)
                                                                               tGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGAACt                                                                                                               >  2:957878/1‑79 (MQ=255)
                                                                               tGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCtgttgt                                                             >  1:1378547/1‑129 (MQ=255)
                                                                               tGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCtgttgt                                                             <  2:1378547/129‑1 (MQ=255)
                                                                                gCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAgg                                                                                                                   <  1:1377792/74‑1 (MQ=255)
                                                                                gCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAgg                                                                                                                   >  2:1377792/1‑74 (MQ=255)
                                                                                gCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCa                                                                                            <  1:896982/97‑1 (MQ=255)
                                                                                gCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCa                                                                                            >  2:896982/1‑97 (MQ=255)
                                                                                     aaGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCTGTTGTTGTCGCCATTATg                                                <  1:665889/136‑1 (MQ=255)
                                                                                     aaGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCTGTTGTTGTCGCCATTATg                                                >  2:665889/1‑136 (MQ=255)
                                                                                         cGAGGCCTATGTATATAAAGCAGCGTCTGGCGTCTTCCGTCGTTGGTCG‑TTTTTTATGTGTAAGGATCTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGCAGATGCGCTGTTGTTTTCGCCATTATGCTGGCTa                                         <  1:1550135/139‑1 (MQ=255)
                                                                                           aGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCtgttgttgt                                                          >  2:46869/1‑120 (MQ=255)
                                                                                           aGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCtgttgttgt                                                          <  1:46869/120‑1 (MQ=255)
                                                                                           aGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCTGTTGTTGTCGCCATTATGCTGGCTATc                                       >  1:512699/1‑139 (MQ=255)
                                                                                              tCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGa                                     <  1:788197/139‑1 (MQ=255)
                                                                                              tCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGAACt                                                                                                               >  1:1001317/1‑64 (MQ=255)
                                                                                              tCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGAACt                                                                                                               <  2:1001317/64‑1 (MQ=255)
                                                                                                  tGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGAACTTCt                                                                                                            <  1:1172270/63‑1 (MQ=255)
                                                                                                  tGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGAACTTCt                                                                                                            >  2:1172270/1‑63 (MQ=255)
                                                                                                  tGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGACCtgt                                <  1:1093190/139‑1 (MQ=255)
                                                                                                       aaaaaGCAGCGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAgg                                                                                                                   >  2:1349711/1‑51 (MQ=255)
                                                                                                       aaaaaGCAGCGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAgg                                                                                                                   <  1:1349711/51‑1 (MQ=255)
                                                                                                        aaaaGCAGCGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGACCTGTTGTTGc                          >  1:1047941/1‑139 (MQ=255)
                                                                                                        aaaaGCAGCGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGACCTGTTGTTGc                          >  1:1606833/1‑139 (MQ=255)
                                                                                                             cagcGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGACCTGTTGTTGCAgggg                     <  2:512699/139‑1 (MQ=255)
                                                                                                              agcGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCTGTTGTTGTCGCCATTATGCt                                              <  1:1359809/113‑1 (MQ=255)
                                                                                                              agcGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCTGTTGTTGTCGCCATTATGCt                                              >  2:1359809/1‑113 (MQ=255)
                                                                                                              agcGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGAcc                                   <  2:571450/124‑1 (MQ=255)
                                                                                                              agcGGCTGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGAcc                                   >  1:571450/1‑124 (MQ=255)
                                                                                                                    tGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATTGGGCGCATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGAc                                    <  2:1118499/117‑1 (MQ=255)
                                                                                                                    tGACGTCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGAc                                    >  1:1118499/1‑117 (MQ=255)
                                                                                                                         tCTTCCGACGTTGGCCG‑TTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGca         <  2:1606833/139‑1 (MQ=255)
                                                                                                                                aCGTTGGCCG‑TTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCtgttgtt                                                            >  1:1420633/1‑81 (MQ=255)
                                                                                                                                aCGTTGGCCG‑TTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCtgttgtt                                                            <  2:1420633/81‑1 (MQ=255)
                                                                                                                                aCGTTGGCCG‑TTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGcca  >  1:580641/1‑139 (MQ=255)
                                                                                                                                aCGTTGGCCG‑TTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGcc   <  1:1404844/138‑1 (MQ=255)
                                                                                                                                aCGTTGGCCG‑TTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGcc   >  2:1404844/1‑138 (MQ=255)
                                                                                                                                   ttGGCCG‑TTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCtg                                                                 <  1:1218399/73‑1 (MQ=255)
                                                                                                                                   ttGGCCG‑TTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCtg                                                                 >  2:1218399/1‑73 (MQ=255)
                                                                                                                                     ggCCG‑TTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGACCTGTTGTTGCAggg                      <  2:864637/114‑1 (MQ=255)
                                                                                                                                     ggCCG‑TTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGACCTGTTGTTGCAggg                      >  1:864637/1‑114 (MQ=255)
                                                                                                                                          |                                                                                                                                 
TGTGAAAAAACACATCCGGTCACCGGGCAACCCGAAAGGAATACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCA  >  NC_000913/3238418‑3238685

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: ‑