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

CAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTC  >  NC_000913/4456967‑4457255
                                                                                                                                           |                                                                                                                                                        
caggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAAc                                                                                                                                             <  2:207180/151‑1 (MQ=255)
 aggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCAGaaa                                                                                                                                              >  2:200661/1‑149 (MQ=255)
   gtggcagaagtggcagaaGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCATCCGAAAatgc                                                                                                                                          >  2:1014390/1‑151 (MQ=255)
      gcagaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTgg                                                                                                                                       <  2:888898/151‑1 (MQ=255)
        agaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCAGAAAATGCTgg                                                                                                                                       >  1:1028343/1‑149 (MQ=255)
        agaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTGGGGGGGTTCATCCGaaaacgct                                                                                                                                         >  2:973235/1‑142 (MQ=255)
           agtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATggg                                                                                                                                  <  1:636357/151‑1 (MQ=255)
           agtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTTGGGGGGGCCCTCAAAAAAAGCTggggggg                                                                                                                                  >  1:577758/1‑146 (MQ=255)
           agtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTGGGGGGGGCCATTAGAAAAACCtgggtggg                                                                                                                                  >  1:913009/1‑151 (MQ=255)
              ggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTg                                                                                                                               <  1:158462/151‑1 (MQ=255)
                cagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCATCAGAAAATGCTGGAAGGGCTgg                                                                                                                              >  2:708405/1‑150 (MQ=255)
                 agaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCg                                                                                                                            <  1:756343/151‑1 (MQ=255)
                    agCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTTGTGGGGGCCATCAGAAACTGCtgggtggggcgggg                                                                                                                            >  1:533849/1‑142 (MQ=255)
                      ccTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCgaagaa                                                                                                                       <  1:262639/151‑1 (MQ=255)
                              cAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGTTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGaa                                                                                                                <  1:529889/150‑1 (MQ=255)
                                  cTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAATGCTGGATGGGGTGGGGAAGAAAGCGGAAGcc                                                                                                             >  2:156312/1‑149 (MQ=255)
                                  cTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGATTTGGGGGGGTCCCCCAGAAAATGCTGGGTGGGGTGGGGAAGAAGGGGGGAGAcgc                                                                                                           >  1:899485/1‑151 (MQ=255)
                                    gCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCcgcg                                                                                                          <  2:678980/150‑1 (MQ=255)
                                    gCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCcgc                                                                                                           <  2:390628/149‑1 (MQ=255)
                                       cacTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAAATGCTGGGTGGGGTGGCGGAGAAGGGGGGAGCcgcg                                                                                                          >  2:1041114/1‑147 (MQ=255)
                                       cacTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACAGCTGGGTGGGGTGGGGGAGAAAGGGGGAGCcgccccg                                                                                                       >  1:611434/1‑146 (MQ=255)
                                         cTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGGTCATCCGAAACTGCTGGGTGGGGTGGCGGAGAAAGCGGGAggc                                                                                                             >  1:779508/1‑140 (MQ=255)
                                           ggAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTTCATCAGAAACTGCTGGATGGGGTGGGGAAGAAGGCGGGAGCCGCGCAGGaaat                                                                                                 >  2:234665/1‑149 (MQ=255)
                                                ggCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGcc                                                                                            <  2:855187/151‑1 (MQ=255)
                                                ggCGATGTGGATACGCTGGCAGATATGGTGATGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCCTCAGAAAATGCTGGGTGGGCTGGCGGAGAAGGCGGGAGCCGCGCAGGAATTCGGcc                                                                                            >  1:805121/1‑151 (MQ=255)
                                                    aTGTGGATACGCTGGCAGATATGGTGGTGAATACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGTGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAg                                                                                        <  1:417850/151‑1 (MQ=255)
                                                        ggATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAAGCGGGAGCCGCGCAGGAAATCGGCCTTAGc                                                                                       >  1:835391/1‑148 (MQ=255)
                                                        ggATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCa                                                                                     <  1:363661/150‑1 (MQ=255)
                                                              gCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTTGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGAtata                                                                               <  2:649244/150‑1 (MQ=255)
                                                                tGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGATTTTGGGGGGGCCCACCGAAAACGCTGGGTGGGGTGGGGGAGAAAGGGGAAACCGCGGAGGAAATCTGccctttgccg                                                                                    >  1:355324/1‑135 (MQ=255)
                                                                 ggCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGt                                                                             <  1:994530/149‑1 (MQ=255)
                                                                  gCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATg                                                                          <  2:1028384/151‑1 (MQ=255)
                                                                  gCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATg                                                                          <  1:750943/151‑1 (MQ=255)
                                                                   cAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCTGGGTGGGGTGGCGGAGAAAGCGGGAGCCGCGCAGGAATTCAGCCTTAGCCAGGAATAATATg                                                                          >  1:106090/1‑150 (MQ=255)
                                                                    aGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCg                                                                        <  1:295414/151‑1 (MQ=255)
                                                                     gATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCCTCAGAAACTGCTGGATGGGGTGGCGGAGAAAGCGGAAGCCGCGCAGGAATTCGGCCTTAGCCCGat                                                                                  >  2:351154/1‑140 (MQ=255)
                                                                       tatGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTGGGGGGGTCCATCAGCAACTGCTGGGTGGGGTGGCGAAGAAGGCGGAAGGCGCCCAGTAATTCAGACTTTACCACGTATTGTAttggt                                                                       >  1:560043/1‑145 (MQ=255)
                                                                        atGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAgg                                                                    <  2:355329/151‑1 (MQ=255)
                                                                          ggtggtGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGcc                                                                  <  2:333662/151‑1 (MQ=255)
                                                                           gtggtgAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCACCAGAAAATGTGGGGTGGGCTGGGGGAGAAGGGGGGAGCCGCGCAGTAAATCGGGCCTAGCCCGAGATAGTATGCGGAg                                                                     >  1:578257/1‑147 (MQ=255)
                                                                                 aaaaCCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAg                                                           <  2:925485/151‑1 (MQ=255)
                                                                                       gCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGCTTTGGTGGGGGCCATCAAAAAATGCTGGGTGGGGTGGGGAAGAAGGCGGGAGCCGCGCAGGAAATCGGCCCTAGCCAGATATAGTATGCGGAGGCCGGGTAAGGCggt                                                      >  1:530108/1‑148 (MQ=255)
                                                                                          cAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGACCCTCAGAAACTGCTGGGTGGGCTGGCGAAGAAGGGGGAAGCCGCGCCGTAATTCGGCCCTAGCCAGATATAGTATGCGTAGGCCGGGTaaagaggt                                                      >  2:333068/1‑141 (MQ=255)
                                                                                          cAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCAc                                                   <  2:345454/150‑1 (MQ=255)
                                                                                          cAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGATCCATCAGAAAATGGTGGATTGGCTGGCGAAGAAAGAGGAAGCCGCCCAGGAATTCGGCCCTAGCCAGATATAGTATTCGTAGGCCGGGTAAAGCGGTCAc                                                   >  2:299132/1‑150 (MQ=255)
                                                                                           aGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGGCCATCAGAAACTGCTGGGTGGGGGGGGGAAGAAAGGGGGAGGCGCGCAGGAAATCGGCCCTAGCCCGATATAGTATGCGTAGGCCGGGTAAGGCGGTCAc                                                   >  1:884499/1‑149 (MQ=255)
                                                                                                     ccATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTACATCACAAACTGCTGGATGGGGTGGGGGAGAAGGCGGGAGCCCCGCCGGGAATTGGGCCTTCCCAGCATTAAGATGCGTAGGCCCGATCAGGGCTTTACGCCCCATACGGc                                       >  1:682586/1‑151 (MQ=255)
                                                                                                      cATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGACTTAATCTCAAACTGTCGGTGGGGCGGGGGAGAAAGGGGAAGCCGCCCAATAATTTTGCCTTTGTGCGACATAGAACGCGTATGCCCGGGAAAGTGTTCTCGCCGCACCCGGc                                       >  1:803533/1‑150 (MQ=255)
                                                                                                      cATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                      <  2:844939/151‑1 (MQ=255)
                                                                                                        tatTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGGCCATCAGAAACTTCTGGGTGGGGTGGCGGAGAAGGCGGGAGCCGCGCAGGAAATTGGCCTTAGCCCGATATAGTATGCGTAGGCCGGGTAAGGCGGTCACGCCGCATCCGGCAtt                                    >  1:920749/1‑151 (MQ=255)
                                                                                                         atTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAtt                                    <  1:335984/150‑1 (MQ=255)
                                                                                                             tGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCTGGGGGGGCTGGCGAAGAAGGGGGGAGCCGCGCAGGAATTTGGCCTTAGCCAGAGATAGTATGCGTAGGGCGGGTAAGGCGTTTACCCCCCATCCCGCAttt                                   >  1:465899/1‑147 (MQ=255)
                                                                                                                 gatgaGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAg                           <  2:770894/151‑1 (MQ=255)
                                                                                                                      gCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCt                      <  1:957559/151‑1 (MQ=255)
                                                                                                                                   ttGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCtt          >  2:1036255/1‑151 (MQ=255)
                                                                                                                                      gtggGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGAt       >  1:225/1‑151 (MQ=255)
                                                                                                                                       tggGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTTTGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGAt       <  2:492244/150‑1 (MQ=255)
                                                                                                                                           aTCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTc  >  1:331157/1‑151 (MQ=255)
                                                                                                                                           |                                                                                                                                                        
CAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTC  >  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: