Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,293,009 G→A 5.9% intergenic (‑87/‑518) hns ← / → tdk global DNA‑binding transcriptional dual regulator H‑NS/thymidine kinase/deoxyuridine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,293,0090GA 5.9% 211.9 / 2.6 69intergenic (‑87/‑518)hns/tdkglobal DNA‑binding transcriptional dual regulator H‑NS/thymidine kinase/deoxyuridine kinase
Reads supporting (aligned to +/- strand):  ref base G (31/32);  new base A (2/2);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.03e-01

ATTCTCTTGCCTGCGCACGAAGAGTACGGATGTTGTTCAGAATTTTAAGTGCTTCGCTCATTGTAGTAATCTCAAACTTATATTGGGGTGGTTTGTTGAGGTAATAATAGAGCCTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAGGAATTCAAAATTGTTCTTTATTTTGTGCCGCCAATAAATATCTTTTCATAAAATTAGCCAGAAAAGACGCGGCATATAGCCCTATTTACACCGATGATTTCGCAGCACGTGAGGTTAAAACTTCCTGATTCATGTCACATTTTATG  >  NC_000913/1292862‑1293155
                                                                                                                                                   |                                                                                                                                                  
aTTCTCTTGCCTGAGCACGAAGAGTACGGATTTTGTTCAGAATTTTAAGTGCTTCGATCATTGTAGTAATCTCAAACTTATATTGGGGTGGTTTGTTGAGGTAATAATAGAGCCTTAAATTCAGTTGTGCAATAGCCAGTAATGTAAg                                                                                                                                                    >  8:104347/1‑148 (MQ=255)
 ttCTCTTGCCTGCGCACGAAGAGTACGGATGTTGTTCAGAATTTTAAGTGCTTCGCTCATTGTAGTAATCTCAAACTTATATTGGGGTGGTTTGTTGAGGTAATAATAGAGCCTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAGGa                                                                                                                                                  <  6:157367/149‑1 (MQ=255)
    tctTGCCTGCGCACGAAGAGTACGGATGTTGTTCAGAATTTTAAGTGCTTCGCTCATTGTAGTAATCTCAAACTTATATTGGGGTGGTTTGTTGAGGTAATAATAGAGCCTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAGGAAtt                                                                                                                                               >  8:254274/1‑149 (MQ=255)
    tctTGCCTGCGCACGAAGAGTACGGATGTTGTTCAGAATTTTAAGTGCTTCGCTCATTGTAGTAATCTCAAACTTATATTGGGGTGGTTTGTTGAGGTAATAATAGAGCCTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAGGAAtt                                                                                                                                               >  4:216676/1‑149 (MQ=255)
          cTGCGCACGAAGAGTACGGATGTTGTTCAGAATTTTAAGTGCTTCGCTCATTGTAGTAATCTCAAACTTATATTGGGGTGGTTTGTTGAGGTAATAATAGAGCCTTAAATTCATTTGTGCAATAGCCAGGAATGTAAGGAATTCAAAAt                                                                                                                                         >  4:261246/1‑149 (MQ=255)
              gcACGAAGAGTACGGATGTTGTTCAGAATTTTAAGTGCTTCGCTCATTGTAGTAATCTCAAACTTATATTGGGGTGGTTTGTTGAGGTAATAATAGAGCCTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAGGAATTCAAAATTGtt                                                                                                                                     <  3:145873/149‑1 (MQ=255)
                aCGAAGAGTACGGATGTTGTTCAGAATTTTAAGTGCTTCGCTCATTGTAGTAATCTCAAACTTATATTGGGGTGGTTAGTTGAGGTAATAATAGAGACTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAGTAATTCAAAATTGTTCt                                                                                                                                   >  4:253815/1‑149 (MQ=255)
                   aaGAGTACGGATTTTGTTCAGAATTTTACGTGCTTCGCTCATTGTAGTAATCTCAAACTTATATTGGGGTGGTTTGTTGAGGTAATAATAGAGCCTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAGGAATTCAAAATTGTTCTTTa                                                                                                                                <  8:137597/149‑1 (MQ=255)
                    agagTACGGATGTTGTTCAGAATTTTAAGTGCTTCGCTCATTGTAGTAATCTCAAACTTATATTGGGGTGGTTTGTTGAGGTAATAATAGAGCCTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAGGAATTCAAAATTGTTCTTTAt                                                                                                                               <  3:216676/149‑1 (MQ=255)
                                ttgttCAGAATTTTAAGTGCTTCGCTCATTGTAGTAATCTCAAACTTATATTGGGGTGGTTTGTTGAGGTAATAATAGAGCCTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAGGAATTCAAAATTGtt                                                                                                                                     <  4:290176/131‑1 (MQ=255)
                                  gttCAGAATTTTAAGTGCTTCGCTCATTGTAGTAATCTCAAACTTATATTGGGGTGGTTTGTTGAGGTAATAATAGAGCCTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAGGAATTCAAAATTGTTCTTTATTTTGTGCCGggta                                                                                                                  >  8:114934‑M2/1‑144 (MQ=255)
                                   ttCAGAATTTTAAGTGCTTCGCTCATTGTAGTAATCTCAAACTTATATTGGGGTGGTTTGTTGAGGTAATAATAGAGCCTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAGGAATTCAAAATTGTTCTTTATTTTGTGCCGggtaat                                                                                                                >  7:124272‑M2/1‑143 (MQ=255)
                                   ttCAGAATTTTAAGTGCTGCGCTCATTGTAGTAATCTCAAACTTATATTGGGGTGGTTTGTTGAGGTATTAATAGAGCCTTAAATTCAGTTGTGCAATAGCCAGGAATTTAAGGAATTCAAAATTGTTCTTTATTTTGTGCCGggtaat                                                                                                                >  6:83968‑M2/1‑143 (MQ=255)
                                     cAGAATTTTAAGTGCTTCGCTCATTGTAGTAATCTCAAACTTATATTGGGGTGGTTTGTTGAGGTAATAATAGAGCCTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAGGAATTCAAAATTGTTCtt                                                                                                                                  <  8:7634/129‑1 (MQ=255)
                                     cAGAATTTTAAGTGCTTCGCTCATTGTAGTAATCTCAAACTTATATTGGGGTGGTTTGTTGAGGTAATAATAGAGCCTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAGGAATTCAAAATTGTTCtt                                                                                                                                  >  7:7634/1‑129 (MQ=255)
                                      aGAATTTTAAGTGCTTCGCTCATTGTAGTAATCTCAAACTTATATTGGGGTGGTTTGTTGAGGTAATAATAGAGCCTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAGGAATTCAAAATTGTTc                                                                                                                                    <  6:202776/126‑1 (MQ=255)
                                      aGAATTTTAAGTGCTTCGCTCATTGTAGTAATCTCAAACTTATATTGGGGTGGTTTGTTGAGGTAATAATAGAGCCTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAGGAATTCAAAATTGTTc                                                                                                                                    >  5:202776/1‑126 (MQ=255)
                                        aaTTTTAAGTGCTTCGCTCATTGTAGTAATCTCAAACTTATATTGGGGTGGTTTGTTGAGGTAATAATAGAGCCTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAGGAATTCAAAATTGTTCTTTATTTTGTGCCGggtaatgactc                                                                                                           >  7:254274‑M2/1‑138 (MQ=255)
                                                gTGCTTCGCTCATTGTAGTAATCTCAAACTTATATTGGGGTGGTTTGTTGAGGTAATAATAGAGCCTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAGGAATTCAAAATTGTTCTTTATTTTGTGCCGggtaatgactccaacttat                                                                                                   <  6:303971‑M2/149‑20 (MQ=255)
                                                gTGCTTCGCTCATTGTAGTAATCTCAAACTTATATTGGGGTGGTTTGTTGAGGTAATAATAGAGCCTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAGGAATTCAAAATTGTTCTTTATTTTGTGCCGggtaatgactccaacttat                                                                                                   >  8:209852‑M2/1‑130 (MQ=255)
                                                 ttcttCGCTCATTGTAGTAATCTGAAACTTATATTGGGGGGGTTTGTTGATGTAATAATAGAGCCTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAGGAATTCAAAATTGTTCTTTATTTTGTGCCGggtaatgactccaacttatt                                                                                                  >  6:141526‑M2/3‑129 (MQ=255)
                                                  gCTTCGCTCATTGTAGTAATCTCAAACTTATATTGGGGTGGTTTGTTGAGGTAATAATAGAGCCTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAAGAATTCAAAATTGTTCTTTATTTTGTgc                                                                                                                        >  3:206871/1‑126 (MQ=255)
                                                  gCTTCGCTCATTGTAGTAATCTCAAACTTATATTGGGGTGGTTTGTTGAGGTAATAATAGAGCCTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAAGAATTCAAAATTGTTCTTTATTTTGTgc                                                                                                                        <  4:206871/126‑1 (MQ=255)
                                                  gCTTCGCTCATTGTAGTAATCTCAAACTTATATTGGGGTGGTTTGTTGAGGTAATAATAGAGCCTTAAATACAGTTGTGCAATAGCCAGGAATGTAAAGAATTCAAAATTGTTCTTTATTTTGTgc                                                                                                                        <  6:148123/126‑1 (MQ=255)
                                                  gCTTCGCTCATTGTAGTAATCTCAAACTTATATTGGGGTGGTTTGTTGAGGTAATAATAGAGCCTTAAATACAGTTGTGCAATAGCCAGGAATGTAAAGAATTCAAAATTGTTCTTTATTTTGTgc                                                                                                                        >  5:148123/1‑126 (MQ=255)
                                                   cTTCGCTCATTGTAGTAATCTCAAACTTATATTGGGGTGGTTTGTTGAGGTAATAATAGAGCCTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAGGAATTCAAAATTGTTCTTTATTTTGTGCCGggtaatgactccaacttattg                                                                                                 >  5:303971‑M2/1‑127 (MQ=255)
                                                   cTTCGCTAATTGTAGTAATCTCAAACTTATATTGGGGTTGTGTGTTGAGGTAATAATAGAGCCTTAACTTCAGTTTTGCATTAGACAGGAATGTAAGGAATTCAAAATTGTTCTTTATTTTGTGCCGtgtaatgactccaacttattg                                                                                                 >  8:112076‑M2/1‑127 (MQ=255)
                                                      cGCTCATTGTAGTAATCTCAAACTTATATTGGGGTGGTTTGTTGAGGTAATAATAGAGCCTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAGGAATTCAAAATTGTTCTTTATTTTGTGCCGggtaatgactccaacttattgata                                                                                              >  3:261246‑M2/1‑124 (MQ=255)
                                                       gCTCATTGTAGTAATCTCAAACTTATATTGGGGTGGTTTGTTGAGGTAATAATAGAGCCTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAGGAATTCAAAATTGTTCTTTATTTTGTGCCGggtaatgactccaacttattgatagt                                                                                            >  2:290895‑M2/1‑123 (MQ=255)
                                                             tGTAGTAATCTCAAACTTATATTGGGGTGGTTTGTTGAGGTAATAATAGAGCCTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAGGAATTCAAAATTGTTCTTTATTTTGTGCCGggtaatgactccaacttatt                                                                                                  <  1:170189‑M2/137‑21 (MQ=255)
                                                             tGTAGTAATCTCAAACTTATATTGGGGTGGTTTGTTGAGGTAATAATAGAGCCTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAGGAATTCAAAATTGTTCTTTATTTTGTGCCGggtaatgactccaacttatt                                                                                                  >  2:170189‑M2/1‑117 (MQ=255)
                                                               tattaATCTGAAACTTATATTGGGGTGGTTTGTTGTTGTAATATTAGAGACTTAAATTCATTTGTGCAATAGCCTGGAATGTAATTAATTCAAAATTGTTCTTTATTTTGTGCCGggtaatgactccaacttattgatagtgttttat                                                                                     >  3:253815‑M2/1‑115 (MQ=255)
                                                                 gtaATCTCAAACTTATATTGGGGTGGTTTGTTGAGGTAATAATAGAGCCTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAGGAATTCAAAATTGTTCTTTATTTTGTGCCGggtaatgactccaacttattgatagtgttttatgt                                                                                   <  7:327094‑M2/148‑36 (MQ=255)
                                                                     tctcAAACTTATATTGGGGTGGTTTGTTGAGGTAATAATAGAGCCTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAGGAATTCAAAATTGTTCTTTATTTTTTGCCGggtaatgactccaacttattgatagtgttttatgtt                                                                                  <  5:24350‑M2/145‑37 (MQ=255)
                                                                      ctcAAACTTATATTGGGGTGGTTTGTTGAGGTAATAATAGAGCCTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAGGAATTCAAAATTGTTCTTTATTTTGTGCCGggtaatgactccaacttattgatagtgttttatgttcaga                                                                              <  7:321954‑M2/148‑41 (MQ=255)
                                                                       tcAAACTTATATTGGGGTGGTTTGTTGAGGTAATAATAGAGCCTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAGGAATTCaaaa                                                                                                                                          <  4:313517/87‑1 (MQ=255)
                                                                       tcAAACTTATATTGGGGTGGTTTGTTGAGGTAATAATAGAGCCTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAGGAATTCaaaa                                                                                                                                          >  3:313517/1‑87 (MQ=255)
                                                                         aaaCTTATATTGGGGTGGTTTGTTGAGGTAATAATAGAGCCTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAGGAATTCAAAATTGTTCTTTATTTTGTGCCGggtaatgactccaacttattgatagtgttttatgttcagataat                                                                          <  8:353280‑M2/149‑45 (MQ=255)
                                                                             ttATATTGGGGTGGTTTGTTGAGGTAATAATAGAGCCTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAGGAATTCAAAATTGTTCTTTATTTTGTGCCGggtaatgactccaacttattgatagtgttttatgttcagataatgccc                                                                      <  6:122469‑M2/149‑49 (MQ=255)
                                                                             ttATATTGGGGTGGTTTGTTGAGGTAATAATAGAGCCTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAGGAATTCAAAATTGTTCTTTATTTTGTGCCGggtaatgactccaacttattgatagtgttttatgttcagataat                                                                          >  2:192984‑M2/1‑101 (MQ=255)
                                                                             ttATATTGGGGTGGTTTGTTGAGGTAATAATAGAGCCTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAGGAATTCAAAATTGTTCTTTATTTTGTGCCGggtaatgactccaacttattgatagtgttttatgttcagataat                                                                          <  1:192984‑M2/145‑45 (MQ=255)
                                                                             ttATATTGGGGTGGTTTGTTGAGGTAATAATAGAGCCTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAGGAATTCAAAATTGTTCTTTATTTTGTGCCGggtaatgactccaacttattgatagtgtt                                                                                         <  5:89824‑M2/130‑30 (MQ=255)
                                                                             ttATATTGGGGTGGTTTGTTGAGGTAATAATAGAGCCTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAGGAATTCAAAATTGTTCTTTATTTTGTGCCGggtaatgactccaacttattgatagtgtt                                                                                         >  6:89824‑M2/1‑101 (MQ=255)
                                                                                tatTGGGGTGGTTTGTTGAGGTAATAATAGAGCCTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAGGAATTCAAAATTGTTCTTTATTTTGTGCCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgat                                                                   >  8:327094‑M2/1‑98 (MQ=255)
                                                                                tatTGGGGTGGTTTGTTGAGGTAATAATAGAGCCTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAGGAATTCAAAATTGTTCTTTATTTTGTGCCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccga                                                                    >  3:2367‑M2/1‑98 (MQ=255)
                                                                                    ggggTGGTTTGTTGAGGTAATAATAGAGCCTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAGGAATTCAAAATTGTTCTTTATTTTGTGCCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgact                                                               <  2:65919‑M2/149‑56 (MQ=255)
                                                                                          gtTTGTTGAGGTAATAATAGAGCCTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAGGAATTCAAAATTGTTCTTTATTTTGTGCCGggtaatgactccaa                                                                                                        >  7:176616‑M2/1‑88 (MQ=255)
                                                                                          gtTTGTTGAGGTAATAATAGAGCCTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAGGAATTCAAAATTGTTCTTTATTTTGTGCCGggtaatgactccaa                                                                                                        <  8:176616‑M2/102‑15 (MQ=255)
                                                                                            ttgttgAGGTAATAATAGAGCCTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAGGAATTCAAAATTGTTCTTTATTTTGTGCCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatg                                                       <  1:55202‑M2/149‑64 (MQ=255)
                                                                                               ttgAGGTAATAATAGAGCCTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAGGAATTCAAAATTGTTCTTTATTTTGTGCCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcag                                                    >  5:122469‑M2/1‑83 (MQ=255)
                                                                                                          ctcgagcctTACATTCAGTTGTGCAATAGCCAGGAATGTAAGGAATTCAAAATTGTTCTTTATTTTGTGCCGtgaaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgat                                          >  6:10759‑M2/4‑72 (MQ=255)
                                                                                                           taGAGCCTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAGGAATTCAAAATTGTTCTTTATTTTGTGCCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattt                                        <  8:102660‑M2/149‑79 (MQ=255)
                                                                                                            agagCCTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAGGAATTCAAAATTGTTCTTTATTTTGTGCCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgatttt                                       <  1:341105‑M2/149‑80 (MQ=255)
                                                                                                              agCCTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAGGAATTCAAAATTGTTCTTTATTTTGTGCCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttg                                      >  2:341105‑M2/1‑68 (MQ=255)
                                                                                                                 cTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAGGAATTCAAAATTGTTCTTTATTTTGTGCCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttg                                                            >  1:155316‑M2/1‑65 (MQ=255)
                                                                                                                 cTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAGGAATTCAAAATTGTTCTTTATTTTGTGCCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttg                                                            <  2:155316‑M2/123‑59 (MQ=255)
                                                                                                                  ttAAATTCAGTTGTGCAATAGCCAGGAATGTAAGGAATTCAAAATTGTTCTTTATTTTGTGCCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaac                                 <  4:122092‑M2/149‑86 (MQ=255)
                                                                                                                      aTTCAGTTGTGCAATAGCCAGGAATGTAAGGAATTCAAAATTGTTCTTTATTTTGTGCCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttg                                                            >  7:335168‑M2/1‑60 (MQ=255)
                                                                                                                      aTTCAGTTGTGCAATAGCCAGGAATGTAAGGAATTCAAAATTGTTCTTTATTTTGTGCCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttg                                                            <  8:335168‑M2/118‑59 (MQ=255)
                                                                                                                      aTTCAGTTGTGAAAAAGCCAGGATTGTAACGAATTCAAAATTGTTCTTTTTTTTGTGCCGggtaatgactccaacttattgatagggttttatgttcagataatgccctctgactttgtcatgaagctccaccgattttgagaacgaca                             >  6:24350‑M2/1‑60 (MQ=255)
                                                                                                                                 cAATAGCCAGGAATGTAAGGAATTCAAAATTGTTCTTTATTTTGTGCCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgt                  <  2:91435‑M2/149‑101 (MQ=255)
                                                                                                                                         aGGAATGTAAGGAATTCAAAATTGTTCTTTATTTTGTGCCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccg          <  7:37999‑M2/149‑109 (MQ=255)
                                                                                                                                           gAATGTAAGGAATTCAAAATTGTTCTTTATTTTGTGCCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtg        <  6:65093‑M2/149‑111 (MQ=255)
                                                                                                                                                tAAGGAATTCAAAATTGTTCTTTATTTTGTGCCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccagg   <  7:233136‑M2/149‑116 (MQ=255)
                                                                                                                                                tAAGGAATTCAAAATTGTTCTTTATTTTGTGCCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccagg   <  5:74712‑M2/149‑116 (MQ=255)
                                                                                                                                                tAAGGAATTCAAAATTGTTCTTTATTTTGTGCCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccagg   <  5:22783‑M2/149‑116 (MQ=255)
                                                                                                                                                tAAGGAATTCAAAATTGTTCTTTATTTTGTGCCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccag    >  5:65093‑M2/1‑34 (MQ=255)
                                                                                                                                                tAAGGAATTCAAAATTGTTCTTTATTTTGTGCCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttg                                      >  3:223181‑M2/1‑34 (MQ=255)
                                                                                                                                                tAAGGAATTCAAAATTGTTCTTTATTTTGTGCCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttg                                      <  4:223181‑M2/114‑81 (MQ=255)
                                                                                                                                                  aGGAATTCAAAATTGTTCTTTATTTTGTGCCGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggt  >  1:65919‑M2/1‑32 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ATTCTCTTGCCTGCGCACGAAGAGTACGGATGTTGTTCAGAATTTTAAGTGCTTCGCTCATTGTAGTAATCTCAAACTTATATTGGGGTGGTTTGTTGAGGTAATAATAGAGCCTTAAATTCAGTTGTGCAATAGCCAGGAATGTAAGGAATTCAAAATTGTTCTTTATTTTGTGCCGCCAATAAATATCTTTTCATAAAATTAGCCAGAAAAGACGCGGCATATAGCCCTATTTACACCGATGATTTCGCAGCACGTGAGGTTAAAACTTCCTGATTCATGTCACATTTTATG  >  NC_000913/1292862‑1293155

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: