Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,457,1050AG47.4% 51.2 / 31.7 57I441V (ATC→GTC) mplUDP‑N‑acetylmuramate:L‑alanyl‑gamma‑D‑glutamyl‑ meso‑diaminopimelate ligase
Reads supporting (aligned to +/- strand):  ref base A (4/26);  new base G (27/0);  total (31/26)
Fisher's exact test for biased strand distribution p-value = 3.11e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGA‑TCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCG‑CAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTC  >  NC_000913/4456967‑4457255
                                                                                                                                           |                                                                                                                                                        
caggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGA‑TCCATCAGAAAc                                                                                                                                             <  2:207180/151‑1 (MQ=255)
 aggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGGGGGG‑TCCATCAGaaa                                                                                                                                              >  2:200661/1‑149 (MQ=255)
   gtggcagaagtggcagaaGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTTGGGGGG‑TCCATCCGAAAatgc                                                                                                                                          >  2:1014390/1‑151 (MQ=255)
      gcagaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGA‑TCCATCAGAAACTGCTgg                                                                                                                                       <  2:888898/151‑1 (MQ=255)
        agaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGGGGGG‑TCCATCAGAAAATGCTgg                                                                                                                                       >  1:1028343/1‑149 (MQ=255)
        agaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑GTTTGGGGGGG‑TTCATCCGaaaacgct                                                                                                                                         >  2:973235/1‑142 (MQ=255)
           agtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGA‑TCCATCAGAAACTGCTGGATggg                                                                                                                                  <  1:636357/151‑1 (MQ=255)
           agtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑GTTTTGGGGGG‑GCCCTCAAAAAAAGCTggggggg                                                                                                                                  >  1:577758/1‑146 (MQ=255)
           agtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑GTTTGGGGGGG‑GCCATTAGAAAAACCtgggtggg                                                                                                                                  >  1:913009/1‑151 (MQ=255)
              ggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGA‑TCCATCAGAAACTGCTGGATGGGCTg                                                                                                                               <  1:158462/151‑1 (MQ=255)
                cagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTTGGGGGG‑TCCATCAGAAAATGCTGGAAGGGCTgg                                                                                                                              >  2:708405/1‑150 (MQ=255)
                 agaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGA‑TCCATCAGAAACTGCTGGATGGGCTGGCg                                                                                                                            <  1:756343/151‑1 (MQ=255)
                    agCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑GTTTTGTGGGG‑GCCATCAGAAACTGCtgggtggggcgggg                                                                                                                            >  1:533849/1‑142 (MQ=255)
                      ccTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGA‑TCCATCAGAAACTGCTGGATGGGCTGGCgaagaa                                                                                                                       <  1:262639/151‑1 (MQ=255)
                              cAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGTTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGA‑TCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGaa                                                                                                                <  1:529889/150‑1 (MQ=255)
                                  cTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGG‑TCCATCAGAAAATGCTGGATGGGGTGGGGAAGAAAGCGGAAGcc                                                                                                             >  2:156312/1‑149 (MQ=255)
                                  cTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑ATTTGGGGGGG‑TCCCCCAGAAAATGCTGGGTGGGGTGGGGAAGAAGGGGGGAGAcgc                                                                                                           >  1:899485/1‑151 (MQ=255)
                                    gCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGA‑TCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCcgcg                                                                                                          <  2:678980/150‑1 (MQ=255)
                                    gCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGA‑TCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCcgc                                                                                                           <  2:390628/149‑1 (MQ=255)
                                       cacTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTTGTGGGG‑TCCATCAGAAAATGCTGGGTGGGGTGGCGGAGAAGGGGGGAGCcgcg                                                                                                          >  2:1041114/1‑147 (MQ=255)
                                       cacTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGG‑TCCATCAGAAACAGCTGGGTGGGGTGGGGGAGAAAGGGGGAGCcgccccg                                                                                                       >  1:611434/1‑146 (MQ=255)
                                         cTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTTGGGGGG‑GTCATCCGAAACTGCTGGGTGGGGTGGCGGAGAAAGCGGGAggc                                                                                                             >  1:779508/1‑140 (MQ=255)
                                           ggAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGG‑TTCATCAGAAACTGCTGGATGGGGTGGGGAAGAAGGCGGGAGCCGCG‑CAGGaaat                                                                                                 >  2:234665/1‑149 (MQ=255)
                                                ggCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGA‑TCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCG‑CAGTAATTCGGcc                                                                                            <  2:855187/151‑1 (MQ=255)
                                                ggCGATGTGGATACGCTGGCAGATATGGTGATGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGG‑TCCCTCAGAAAATGCTGGGTGGGCTGGCGGAGAAGGCGGGAGCCGCG‑CAGGAATTCGGcc                                                                                            >  1:805121/1‑151 (MQ=255)
                                                    aTGTGGATACGCTGGCAGATATGGTGGTGAATACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGTGA‑TCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCG‑CAGTAATTCGGCCTTAg                                                                                        <  1:417850/151‑1 (MQ=255)
                                                        ggATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGG‑TCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAAGCGGGAGCCGCG‑CAGGAAATCGGCCTTAGc                                                                                       >  1:835391/1‑148 (MQ=255)
                                                        ggATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGA‑TCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCG‑CAGTAATTCGGCCTTAGCCa                                                                                     <  1:363661/150‑1 (MQ=255)
                                                              gCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGA‑TCCATCAGAAACTGCTTGATGGGCTGGCGAAGAAGGCGGAAGCCGCG‑CAGTAATTCGGCCTTAGCCAGAtata                                                                               <  2:649244/150‑1 (MQ=255)
                                                                tGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGATTTTGGGGGGG‑CCCA‑CCGAAAACGCTGGGTGGGGTGGGGGAGAAAGGGGAAACCGCG‑GAGGAAATCTGccctttgccg                                                                                    >  1:355324/1‑135 (MQ=255)
                                                                 ggCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGA‑TCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCG‑CAGTAATTCGGCCTTAGCCAGATATAGt                                                                             <  1:994530/149‑1 (MQ=255)
                                                                  gCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGA‑TCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCG‑CAGTAATTCGGCCTTAGCCAGATATAGTATg                                                                          <  2:1028384/151‑1 (MQ=255)
                                                                  gCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGA‑TCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCG‑CAGTAATTCGGCCTTAGCCAGATATAGTATg                                                                          <  1:750943/151‑1 (MQ=255)
                                                                   cAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGG‑TCCATCAGAAACTGCTGGGTGGGGTGGCGGAGAAAGCGGGAGCCGCG‑CAGGAATTCAGCCTTAGCCAGGAATAATATg                                                                          >  1:106090/1‑150 (MQ=255)
                                                                    aGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGA‑TCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCG‑CAGTAATTCGGCCTTAGCCAGATATAGTATGCg                                                                        <  1:295414/151‑1 (MQ=255)
                                                                     gATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGG‑TCCCTCAGAAACTGCTGGATGGGGTGGCGGAGAAAGCGGAAGCCGCG‑CAGGAATTCGGCCTTAGCCCGat                                                                                  >  2:351154/1‑140 (MQ=255)
                                                                       tatGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑GTTTGGGGGGG‑TCCATCAGCAACTGCTGGGTGGGGTGGCGAAGAAGGCGGAAGGCGCC‑CAGTAATTCAGACTTTACCACGTATTGTAttggt                                                                       >  1:560043/1‑145 (MQ=255)
                                                                        atGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGA‑TCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCG‑CAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAgg                                                                    <  2:355329/151‑1 (MQ=255)
                                                                          ggtggtGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGA‑TCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCG‑CAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGcc                                                                  <  2:333662/151‑1 (MQ=255)
                                                                           gtggtgAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGG‑TCCACCAGAAAATGTGGGGTGGGCTGGGGGAGAAGGGGGGAGCCGCG‑CAGTAAATCGGGCCTAGCCCGAGATAGTATGCGGAg                                                                     >  1:578257/1‑147 (MQ=255)
                                                                                 aaaaCCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGA‑TCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCG‑CAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAg                                                           <  2:925485/151‑1 (MQ=255)
                                                                                       gCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑CTTTGGTGGGG‑GCCATCAAAAAATGCTGGGTGGGGTGGGGAAGAAGGCGGGAGCCGCG‑CAGGAAATCGGCCCTAGCCAGATATAGTATGCGGAGGCCGGGTAAGGCggt                                                      >  1:530108/1‑148 (MQ=255)
                                                                                          cAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGG‑ACCCTCAGAAACTGCTGGGTGGGCTGGCGAAGAAGGGGGAAGCCGCG‑CCGTAATTCGGCCCTAGCCAGATATAGTATGCGTAGGCCGGGTaaagaggt                                                      >  2:333068/1‑141 (MQ=255)
                                                                                          cAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGA‑TCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCG‑CAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCAc                                                   <  2:345454/150‑1 (MQ=255)
                                                                                          cAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGGGGGA‑TCCATCAGAAAATGGTGGATTGGCTGGCGAAGAAAGAGGAAGCCGCC‑CAGGAATTCGGCCCTAGCCAGATATAGTATTCGTAGGCCGGGTAAAGCGGTCAc                                                   >  2:299132/1‑150 (MQ=255)
                                                                                           aGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGG‑GCCATCAGAAACTGCTGGGTGGGGGGGGGAAGAAAGGGGGAGGCGCG‑CAGGAAATCGGCCCTAGCCCGATATAGTATGCGTAGGCCGGGTAAGGCGGTCAc                                                   >  1:884499/1‑149 (MQ=255)
                                                                                                     ccATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGG‑TACATCACAAACTGCTGGATGGGGTGGGGGAGAAGGCGGGAGCCCCGCCGGGAATTGGGCCTT‑CCCAGCATTAAGATGCGTAGGCCCGATCAGGGCTTTACGCCCCATACGGc                                       >  1:682586/1‑151 (MQ=255)
                                                                                                      cATATTCTGGTGATGAGCAACGGCGG‑TTTTTGTGGGACTTAATCTCAAACTG‑TCGGTGGGGCGGGGGAGAAAGGGGAAGCCGCC‑CAATAATTTTGCCTTTGTGCGACATAGAACGCGTATGCCCGGGAAAGTGTTCTCGCCGCACCCGGc                                       >  1:803533/1‑150 (MQ=255)
                                                                                                      cATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGA‑TCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCG‑CAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                      <  2:844939/151‑1 (MQ=255)
                                                                                                        tatTCTGGTGATGAGCAACGGCGG‑TTTTTGTGGGG‑GCCATCAGAAACTTCTGGGTGGGGTGGCGGAGAAGGCGGGAGCCGCG‑CAGGAAATTGGCCTTAGCCCGATATAGTATGCGTAGGCCGGGTAAGGCGGTCACGCCGCATCCGGCAtt                                    >  1:920749/1‑151 (MQ=255)
                                                                                                         atTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGA‑TCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCG‑CAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAtt                                    <  1:335984/150‑1 (MQ=255)
                                                                                                             tGGTGATGAGCAACGGCGG‑TTTTGGTGGGG‑TCCATCAGAAACTGCTGGGGGGGCTGGCGAAGAAGGGGGGAGCCGCG‑CAGGAATTTGGCCTTAGCCAGAGATAGTATGCGTAGGGCGGGTAAGGCGTTTACCCCCCATCCCGCAttt                                   >  1:465899/1‑147 (MQ=255)
                                                                                                                 gatgaGCAACGGCGG‑TTTTGGTGGGA‑TCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCG‑CAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAg                           <  2:770894/151‑1 (MQ=255)
                                                                                                                      gCAACGGCGG‑TTTTGGTGGGA‑TCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCG‑CAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCt                      <  1:957559/151‑1 (MQ=255)
                                                                                                                                   ttGGTGGGA‑TCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCG‑CAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCtt          >  2:1036255/1‑151 (MQ=255)
                                                                                                                                      gtggGA‑TCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCG‑CAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGAt       >  1:225/1‑151 (MQ=255)
                                                                                                                                       tggGA‑TCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCG‑CAGTAATTCGGCCTTAGCCAGATATAGTTTGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGAt       <  2:492244/150‑1 (MQ=255)
                                                                                                                                           a‑TCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCG‑CAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTc  >  1:331157/1‑151 (MQ=255)
                                                                                                                                           |                                                                                                                                                        
CAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGG‑TTTTGGTGGGA‑TCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCG‑CAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTC  >  NC_000913/4456967‑4457255

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