Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,162,830 G→A 6.8% E65K (GAA→AAA)  lpoB → OM lipoprotein stimulator of MrcB transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,162,8300GA 6.8% 180.0 / 2.7 59E65K (GAA→AAA) lpoBOM lipoprotein stimulator of MrcB transpeptidase
Reads supporting (aligned to +/- strand):  ref base G (31/23);  new base A (2/2);  total (34/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.37e-01

TCGCCGGGTGTGTGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAGCACGAAGATCAAACTGCACCGCCTGCGCCGCATATTCGCCATTATGACTGGAATGGCGCAATGCAGCCGATGGTCAGTAAGATGCTTGGGGCTGACGGGGTGACTGCGGGTAGCGTCCTGCTGGTTGATAGCGTTAACAACCGTACTAACGG  >  NC_000913/1162687‑1162978
                                                                                                                                               |                                                                                                                                                    
tcGCCGGGTGTGTGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAGCACg                                                                                                                                                      >  7:49755/1‑144 (MQ=255)
tcGCCGGGTGTGTGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCACCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAGCACGAAGAt                                                                                                                                                 <  3:67884/149‑1 (MQ=255)
   ccGGGTGTGTGGGTCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACCACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCACGCCCAATTGAGCACGAAGATCa                                                                                                                                               >  3:191344/1‑148 (MQ=255)
   ccGGGTGTGTGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAGCACGAAGATCa                                                                                                                                               >  5:93953/1‑148 (MQ=255)
    cGGGGGTGTGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCACCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAGCACGAAGATCaaa                                                                                                                                             <  5:194484/149‑1 (MQ=255)
      ggTGTGTGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAGCACGAAGATCAAACt                                                                                                                                           <  4:205045/149‑1 (MQ=255)
         gtgtgGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAGCACGAAGATCAAACTGCa                                                                                                                                        <  7:179119/149‑1 (MQ=255)
                  aaCGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAGCACGAAGATCAAACTGCACCGCCTgcg                                                                                                                               >  7:334906/1‑149 (MQ=255)
                    cgggAACCTGCCCCGGTAGAAGAAGTGAACCCAGCGCCGGAACACCCAGCCGAGCCACAACAGCCTGGCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAGCACGAAGATCAAACTGCACCGCCTGcgcc                                                                                                                             <  7:350346/146‑1 (MQ=255)
                     gTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAGCACGAAGATCAAACTGCACCGCCTGcgccg                                                                                                                            >  8:124405/1‑149 (MQ=255)
                      tGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACACCAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAGCACGAAGATCAAACTGCACCGCCTGcgccg                                                                                                                            >  2:180230/1‑148 (MQ=255)
                       gAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAGCACGAAGATCAAACTGCACCGCCTGCGCCGCa                                                                                                                          >  7:75354/1‑149 (MQ=255)
                       gAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAGCACGAAGATCAAACTGCACCGCCTGCGCCGCa                                                                                                                          >  8:116434/1‑149 (MQ=255)
                           cTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGACCAATTGAGCACGAAGATCAAACTGCACCGCCTGCGCCGCatat                                                                                                                       >  6:215666/1‑148 (MQ=255)
                              cACCGGTAGAAGAAGTGAACCCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAGCACGAAGATCAAACTGCAACGCCTGCGCCGCATAATCGc                                                                                                                   >  6:209033/1‑149 (MQ=255)
                               cccGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCAAGGCCCAATAGAGCACTAAGATCAAACTGCACCGCCTGCGCCGCatataccc                                                                                                                   >  4:175471/2‑144 (MQ=255)
                                   gTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGCGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAGCACGAAGATCAAACTGCACCGCCTGCGCCGCATATTCGCCATTa                                                                                                              <  6:93953/149‑1 (MQ=255)
                                    tAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAGCACGAAGATCAAACTGCACCGCCTGCGCCGCATATTCGCCATTAt                                                                                                             >  3:119354/1‑149 (MQ=255)
                                       aagaagTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAGCACGAAGATCAAACTGCACCGCCTGCGCCGCATATTCGCCATTATGAc                                                                                                          >  8:167015/1‑149 (MQ=255)
                                           agTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGCCGATCCCGCAGCAGCCAGGCCCAATTGAGCACAAAGATCAAACTGCACCGCCTGCGCCGCATATTCGCCATTATGACTgg                                                                                                       >  6:7722/1‑148 (MQ=255)
                                           agTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAGCACAAAGATCAAACTGCACCGCCTGCGCCGCATATTCGCCATTATGACTgg                                                                                                       >  6:2879/1‑148 (MQ=255)
                                            gTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAGCACGAAGATCAAACTGCACCGCCTGCGCCGCATATTCGCCATTATGACTGGaa                                                                                                     <  8:334906/149‑1 (MQ=255)
                                                      cgcCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAGCACGAAGATCAAACTGCACCGCCTGCGCCGCATATTCGCCATTATGACTGGAATGGAGCAATg                                                                                           >  5:260208/1‑149 (MQ=255)
                                                            aacaacCCGCCGAGCCACAACAGCCAGACCCCACAGACCCCTCGGAGCCGACGATCCCGCAGCACCCAGGCCCACGAGAGCAGGAAGATCAAACTGCCCCGCCTGCGCCGCATATTCGCCATTATGACTGGACTGGCGCCAGGCAGCCg                                                                                     >  5:170620/1‑149 (MQ=255)
                                                            aacaacCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAGCACGAAGATCAAACTGCACCGCCTGCGCCGCATATTCGCCATTATGACTGGAATGGCGCAATGCAGcc                                                                                      >  5:325508/1‑148 (MQ=255)
                                                                            acaacaGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAGCACGAAGATCAAACTGCACCGCCTGCGCCGCATATTCGCCATTATGACTGGAATGGCGCAATGCAGCCGATGGTCAGTAAGATGc                                                                     >  1:193056/1‑149 (MQ=255)
                                                                             caacaGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAGCACGAAGATCAAACTGCACCGCCTGCGCCGCATATTCGCCATTATGACTGGAATGGCGCAATGCAGCCGATGGTCAGTAAGATGCt                                                                    >  5:76785/1‑149 (MQ=255)
                                                                              aacaGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAGCACGAAGATCAAACTGCACCGCCTGCGCCGCATATTCGCCATTATGACTGGAATGGCGCAATGCAGCCGATGGTCAGTAAGATGCtt                                                                   >  5:154487/1‑149 (MQ=255)
                                                                              aacaGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAGCACGAAGATCAAACTGCACCGCCTGCGCCGCATATTCGCCATTATGACTGGAATGGCGCAATGCAGCCGATGGTCAGTAAGATGCtt                                                                   >  3:219796/1‑149 (MQ=255)
                                                                                 aGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAGCACGAAGATCAAACTGCACCGCCTGCGCCGCATATTCGCCATTATGACTGGAATGGCGCAATGCAGCCGATGGTCAGTAAGATGCTTggg                                                                <  1:171946/149‑1 (MQ=255)
                                                                                  gCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAGCACGAAGATCAAACTGCACCGCCTGCGCCGCATATTCGCCATTATGACTGGAATGGCGCAATGCAGCCGATGGTCAGTAAGATGCTTgggg                                                               <  7:124405/149‑1 (MQ=255)
                                                                                    cTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAGCACGAAGATCAAACTGCACCGCCTGCGCCGCATATTCGCCATTATGACTGGAATGGCGCAATGCAGCCGATGGTCAGTAAGATGCTTGGGGCt                                                             >  3:33655/1‑149 (MQ=255)
                                                                                         cccACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAGCACGAAGATCAAACTGCACCGCCTGCGCCGCATATTCGCCATTATGACTGGAATGGCGCAATGCAGCCGATGGTCAGTAAGATGCTTGGGGGTGACgg                                                        >  2:67931/1‑149 (MQ=255)
                                                                                         cccACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAGCACGAAGATCAAACTGCACCGCCTGCGCCGCATATTCGCCATTATGACTGGAATGGCGCAATGCAGCCGATGGTCAGTAAGATGCTTGGGGCTGACgg                                                        >  1:5788/1‑149 (MQ=255)
                                                                                          ccACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGTCCCAATTGAGCACGAAGATCAAACTGCACCGCCTGCGCCGCATATTCGCCATTATGACTGGAATGGCGCAATGCAGCCGATGGTCAGTAAGATGCTTGGGGCTGACggg                                                       <  5:209033/149‑1 (MQ=255)
                                                                                            acaGTGCCCGCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGCGCACGAAGATCAAACTGCACCGCCTGCGCCGCATATGCGCCATTATGACTGGAATGGCGCAATGCAGCCGATGGTCAGTAAGATGCTTGGGGCTGACggggg                                                     <  5:336227/149‑2 (MQ=255)
                                                                                              aGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAGCACGAAGATCAAACTGCACCGCCTGCGCCGCATATTCGCCATTATGACTGGAATGGCGCAATGCAGCCGATGGTCAGTAAGATGCTTGGGGCTGAc                                                          >  2:92565/1‑142 (MQ=255)
                                                                                              aGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAGCACGAAGATCAAACTGCACCGCCTGCGCCGCATATTCGCCATTATGACTGGAATGGCGCAATGCAGCCGATGGTCAGTAAGATGCTTGGGGCTGAc                                                          <  1:92565/142‑1 (MQ=255)
                                                                                                   ccTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAGCACGAAGATCAAACTGCACCGCCTGCGCCGCATATTCGCCATTATGACTGGAATGGCGCAATGCAGCCGATGGTCAGTAAGATGCTTGGGGCTGACGGGGTGACTGCg                                              >  4:156470/1‑149 (MQ=255)
                                                                                                   ccTCGGTGACGACGATCCCGCAGCAGCCTGGCCCAATTGAGCACGAAGATCAAACTGCACCGCCTGCGCCGCATATTCGCCATTATGACTGGAATGGCGCAATGCAGCCGATGGTCAGTAAGATGCTTGGGGCTGACGGGGTGACTGCg                                              >  2:23598/1‑149 (MQ=255)
                                                                                                       ggTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAGCACGAAGATCAAACTGCACCGCCTGCGCCGCATATTCGCCATTATGACTGGAATGGCGCAATGCAGCCGATGGTCAGTAAGATGCTTGGGGCTGACGGGGTGACTGCggg                                            <  7:351/147‑1 (MQ=255)
                                                                                                          gCCGAAGATCCCGCAGCCGCCAGGCCCAATTGAGCACGAAGATCAAACTGCCCCGCCTGCGCCGCATATTCGCCATTATGAGTGGAATGGCGCAATGCAGCCGATGGTCAGGAAGATGATTGGGGCTGACGGGGTGAGTGAGGGTAGCg                                       >  6:75190/1‑149 (MQ=255)
                                                                                                           ccGACGATCCCGCTGCAGGCAGGCCCAATTGAGCACGAAGATCAAAATGCACCGCCTGCGCCGCATATTCGCCATTATGACTGGAATGGCGCAATGCAGCCGATGGTCAGTAAGATGCTTTGGGCTGACGGGGTGACTGCGGGTAGCg                                       <  4:191344/148‑1 (MQ=255)
                                                                                                             gacgaTCCCGCAGCAGCCAGGCCCAATTGAGCACGAAGATCAAACTGCACCGCCTGCGCCGCATATTCGCCATTATGACTGGAATGGCGCAATGCAGCCGATGGTCAGTAAGATGCTTGGGGCTGACGGGGTGACTGCGGGTAGCGTcc                                    >  5:197863/1‑149 (MQ=255)
                                                                                                                gaTCCCGCCGCCGCCAGGCCCAATTGAGCCCGAAGATCAAACTGCACCGCCTGCGCCGCATATTCGCCATGATGACTGGAATGGCGCAATGCAGCCGAGGGTCAGTAAGATGCTTGGGGCTGCCGGGGTGACTGCGGGTCGCGTCgtgc                                 >  2:138513/1‑149 (MQ=255)
                                                                                                                     cGCAGCAGCCAGGCCCAATTGAGCACGAAGATCAAACTGCACCGCCTGCGCCGCATATTCGCCATTATGACTGGAATGGCGCAATGCAGCCGATGGTCAGTAAGATGCTTGGGGCTGACGGGGTGACTGCGGGTAGCGTCCTGCTGGtt                            <  2:193056/149‑1 (MQ=255)
                                                                                                                      gcagcaGCCAGGCCCAATTGAGCACAAAGATCAAACTGCACCGCCTGCGCCGCATATTCGCCATTATGACTGGAATGGCGCAATGCAGCCGATGGTCAGTAAGATGCTTGGGGCTGACGGGGTTACTGCGGGTAGCGTCCTGCTGGTTg                           <  5:7722/149‑1 (MQ=255)
                                                                                                                      gcagcaGCCAGGCCCAATTGAGCACAAAGATCAAACTGCACCGCCTGCGCCGCATATTCGCCATTATGACTGGAATGGCGCAATGCAGCCGATGGTCAGTAAGATGCTTGGGGCTGACGGGGTTACTGCGGGTAGCGTCCTGCTGGTTg                           <  5:2879/149‑1 (MQ=255)
                                                                                                                         gcagcCAGGCCCAATTGAGCACGAAGATCAAACTGCACCGCCTGCGCCGCATATTCGCCATTATGACTGGAATGGCGCAATGCAGCCGATGGTCAGTAAGATGCTTGGGGCTGACGGGGTGACTGCGGGTAGCGTCCTGCTGGTTGATa                        <  4:119354/149‑1 (MQ=255)
                                                                                                                             ccAGGCCCAATTGAGCACGAAGATCAAACTGCACCGCCTGCGCCGCATATTCGCCATTATGACTGGAATGGCGCAATGCAGCCGATGGTCAGTAAGATGCTTGGGGCTGACGGGGTGACTGCGGGTAGCGTCCTGCTGGTTGATAGCGt                    <  1:180230/149‑1 (MQ=255)
                                                                                                                             ccAGGCCCAATTGAGCACGAAGATCAAACTGCACCGCCTGCGCCGCATATTCGCCATTATGACTGGAATGGCGCAATGCAGCCGATGGTCAGTAAGATGCTTGGGGCTGACGGGGTGACTGCGGGTAGCGTCCTGCTGGTTGATAGCGt                    <  7:167015/149‑1 (MQ=255)
                                                                                                                             ccAGGCCCAATTGAGCACGAAGATCAAACTGCACCGCCGGCGCCGCATTTTCGCCATTATGACTGGAATGGCGCAATGCAGCCGATGGGCAGTAAGATGCTTGGGGCTGACGGGGGGACTGCGGGTAGCGTCCTGCTGGTTGATAGCGt                    <  3:175471/149‑1 (MQ=255)
                                                                                                                              cAGGCCCAATTGAGCACGAAGATCAAACTGCACCGCCTGCGCCGCATATTCGCCATTATGACTGGAATGGCGCAATGCAGCCGACGGTCAGTAAGATGCTTGGGGCTGACGGGGTGACTGCGGGTAGCGTCCTGCTGGTTGATAGCGt                    >  7:20512/1‑148 (MQ=255)
                                                                                                                               aGGCCCAATTGAGCACGAAGATCAAACTGCACCGCCTGCGCCGCATATTCGCCATTATGACTGGAATGGCGCAATGCAGCCGATGGTCAGTAAGATGCTTGGGGCTGACggg                                                       <  7:207952/112‑1 (MQ=255)
                                                                                                                               aGGCCCAATTGAGCACGAAGATCAAACTGCACCGCCTGCGCCGCATATTCGCCATTATGACTGGAATGGCGCAATGCAGCCGATGGTCAGTAAGATGCTTGGGGCTGACggg                                                       >  8:207952/1‑112 (MQ=255)
                                                                                                                                  cccAATTGAGCACGAAGATCAAACTGCACCGCCTGCGCCGCATATTCGCCATTATGACTGGAATGGCGCAATGCAGCCGATGGTCAGTAAGATGCTTGGGGCTGACGGGGTGACTGCGGGTAGCGTCCTGCTGGTTGATAGCGTTaac                >  6:130181/1‑148 (MQ=255)
                                                                                                                                   ccAATTGAGCACGAAGATCAAACTGCACCGCCTGCGCCGCATATTCGCCATTATGACTGGAATGGCGCAATGCAGCCGATGGTCAGTAAGATGCTTGGGGCTGACGGGGTGACTGCGGGTAGCGTCCTGCTGGTTGATAGCGTTaacaa              <  3:297636/149‑1 (MQ=255)
                                                                                                                                    cAATTGAGCACGAAGATCAAACTGCACCGCCTGCGCCGCATATTCGCCATTATGACTGGAATGGCGCAATGCAGCCGATGGTCAGTAAGATGCTTGGGGCTGACGGGGTGACTGCGGGTAGCGTCCTGCTGGTTGATAGCGTTaacaac             <  8:49755/149‑1 (MQ=255)
                                                                                                                                         gAGCACGAAGATCAAACTGCACCGCCTGCGCCGCATATTCGCCATTATGACTGGAATGGCGCAATGCAGCCGATGGTCAGTAAGATGCTTGGGGCTGACGGGGTGACTGc                                               >  2:354513/1‑110 (MQ=255)
                                                                                                                                         gAGCACGAAGATCAAACTGCACCGCCTGCGCCGCATATTCGCCATTATGACTGGAATGGCGCAATGCAGCCGATGGTCAGTAAGATGCTTGGGGCTGACGGGGTGACTGc                                               <  1:354513/110‑1 (MQ=255)
                                                                                                                                               gAAGATCAAACTGCACCGCCTGCGCCGCATATTCGCCATTATGACTGGAATGGCGCAATGCAGCCGATGGTCAGTAAGATGCTTGGGGCTGACGGGGTGACTGCGGGTAGCGTCCTGCTGGTGGATAGCGTTAACACCCGTACTAACgg  >  5:177530/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
TCGCCGGGTGTGTGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAGCACGAAGATCAAACTGCACCGCCTGCGCCGCATATTCGCCATTATGACTGGAATGGCGCAATGCAGCCGATGGTCAGTAAGATGCTTGGGGCTGACGGGGTGACTGCGGGTAGCGTCCTGCTGGTTGATAGCGTTAACAACCGTACTAACGG  >  NC_000913/1162687‑1162978

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