Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,329,894 (G)6→7 coding (535/879 nt) yjcZ → YjcZ family protein, yhjH motility defect suppressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,329,8881.G98.4% 229.7 / ‑3.6 64coding (529/879 nt)yjcZYjcZ family protein, yhjH motility defect suppressor
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (33/30);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

GCCTACTGGTGCAGGCAGGTGATCAGAGCGAGCAAATTTTTACCCGCTGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCAGAGAGTATCAATGCGGCACCGACCACGCGTCACTA  >  NC_000913/4329772‑4330018
                                                                                                                     |                                                                                                                                    
gCCTACTGGTGCAGGCAGGTGATCAGAGCGAGCAAATTTTTACCCGCTGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGAtgcct                                                                                                               <  2:1435153/139‑1 (MQ=255)
          gcaggcagGTGATCAGAGCGAGCAAATTTTTACCCGCTGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTggg                                                                                                     >  1:294576/1‑139 (MQ=255)
             ggcaggTGATCAGAGCGAGCAAATTTTTACCCGCTGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑‑gggg                                                                                                                                >  2:1383867/1‑108 (MQ=255)
             ggcaggTGATCAGAGCGAGCAAATTTTTACCCGCTGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑‑gggg                                                                                                                                <  1:1383867/108‑1 (MQ=255)
              gcaggTGATCAGAGCGAGCAAATTTTTACCCGCTGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGc                                                                                                 >  2:1265182/1‑139 (MQ=255)
              gcaggTGATCAGAGCGAGCAAATTTTTACCCGCTGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑‑ggggg                                                                                                                               <  2:27497/108‑1 (MQ=255)
              gcaggTGATCAGAGCGAGCAAATTTTTACCCGCTGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑‑ggggg                                                                                                                               >  1:27497/1‑108 (MQ=255)
                     aTCAGAGCGAGCAAATTTTTACCCGCTGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAg                                                                                          <  2:294576/139‑1 (MQ=255)
                          agcgagcAAATTTTTACCCGCTGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTg                                                                                     <  2:91018/139‑1 (MQ=255)
                          agcgagcAAATTTTTACCCGCTGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTg                                                                                     >  2:210852/1‑139 (MQ=255)
                            cgagcAAATTTTTACCCGCTGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGc                                                                                   >  2:1059030/1‑139 (MQ=255)
                            cgagcAAATTTTTACCCGCTGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGc                                                                                   >  2:361232/1‑139 (MQ=255)
                             gagcAAATTTTTACCCGCTGGCGGGCTGGGGCGTATAACGACTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCg                                                                                  <  1:361232/139‑1 (MQ=255)
                              agcAAATTTTTACCCGCTGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGAtgcctg                                                                                                              >  1:336811/1‑110 (MQ=255)
                              agcAAATTTTTACCCGCTGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGAtgcctg                                                                                                              <  2:336811/110‑1 (MQ=255)
                               gcAAATTTTTACCCGCTGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGAtgcct                                                                                                               >  2:749405/1‑108 (MQ=255)
                               gcAAATTTTTACCCGCTGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGAtgcct                                                                                                               <  1:749405/108‑1 (MQ=255)
                               gcAAATTTTTACCCGCTGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCgtt                                                                                >  2:600334/1‑139 (MQ=255)
                               gcAAATTTTTACCCGCTGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCgtt                                                                                >  1:1115057/1‑139 (MQ=255)
                                          cccGCTGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTg                                                                     <  2:1115057/139‑1 (MQ=255)
                                             gCTGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGAtt                                                                  >  1:29262/1‑139 (MQ=255)
                                              cTGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGAttt                                                                 >  2:258616/1‑139 (MQ=255)
                                               tGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGttgttg                                                                            <  2:1114613/127‑1 (MQ=255)
                                               tGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGttgttg                                                                            >  1:1114613/1‑127 (MQ=255)
                                                         gggCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAg                                                      <  1:600334/139‑1 (MQ=255)
                                                          ggCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGc                                                     <  2:29262/139‑1 (MQ=255)
                                                           gCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCg                                                    >  1:423166/1‑139 (MQ=255)
                                                             gTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGac                                                  >  1:218075/1‑139 (MQ=255)
                                                              tataACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGAtgcct                                                                                                               >  1:954797/1‑77 (MQ=255)
                                                              tataACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGAtgcct                                                                                                               <  2:954797/77‑1 (MQ=255)
                                                              tataACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGaca                                                 <  2:475549/139‑1 (MQ=255)
                                                                taACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGAtgcctg                                                                                                              >  2:1137064/1‑76 (MQ=255)
                                                                taACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGAtgcctg                                                                                                              <  1:1137064/76‑1 (MQ=255)
                                                                taACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGAcaca                                               <  2:218075/139‑1 (MQ=255)
                                                                taACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGAcaca                                               <  2:1364874/139‑1 (MQ=255)
                                                                taACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑‑GGGCGGGCATTTGGCGAtgcctg                                                                                                              <  1:1005451/76‑1 (MQ=255)
                                                                taACGCCTGGTCGTTACCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGAtgcctg                                                                                                              >  2:1005451/1‑76 (MQ=255)
                                                                  aCGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACacaac                                             <  1:210852/139‑1 (MQ=255)
                                                                     ccTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTg                                                                         <  2:568248/108‑1 (MQ=255)
                                                                     ccTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTg                                                                         >  1:568248/1‑108 (MQ=255)
                                                                          tCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCa                                     <  1:1059030/139‑1 (MQ=255)
                                                                             ttGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCagag                                  >  2:642952/1‑139 (MQ=255)
                                                                               gCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGAtgcctg                                                                                                              >  1:841418/1‑61 (MQ=255)
                                                                               gCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGAtgcctg                                                                                                              <  2:841418/61‑1 (MQ=255)
                                                                               gCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCagagag                                >  2:1251173/1‑139 (MQ=255)
                                                                                  tGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCAGAGAGTAt                             <  1:642952/139‑1 (MQ=255)
                                                                                   gggCGCTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGAtgcct                                                                                                               >  1:535732/1‑56 (MQ=255)
                                                                                   gggCGCTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGAtgcct                                                                                                               <  2:535732/56‑1 (MQ=255)
                                                                                   gggCGCTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCAGAGAGTATc                            >  2:598975/1‑139 (MQ=255)
                                                                                    ggCGCTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCAGAGAGTATCa                           >  2:164087/1‑139 (MQ=255)
                                                                                      cgcTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGTAGAGAGTAt                             <  1:1012055/135‑1 (MQ=255)
                                                                                      cgcTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCAGAGAGTAt                             >  2:1012055/1‑135 (MQ=255)
                                                                                        cTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCAGAGAGTATCAATGc                       >  2:204325/1‑139 (MQ=255)
                                                                                         tGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCAGAGAGTATCAATGCg                      >  1:3700/1‑139 (MQ=255)
                                                                                          gTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGAcaca                                               <  2:605088/113‑1 (MQ=255)
                                                                                          gTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGAcaca                                               >  1:605088/1‑113 (MQ=255)
                                                                                           ttttATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCAGAGAGTATCAATGCGGc                    <  2:3700/139‑1 (MQ=255)
                                                                                               aTCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCAGAGAGTATCAATGCGGCa                   <  2:34074/136‑1 (MQ=255)
                                                                                               aTCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCAGAGAGTATCAATGCGGCa                   >  1:34074/1‑136 (MQ=255)
                                                                                                      tGGAGGAGTTGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCAGAGAGTATCAATGCGGCACCGACCAcgc         <  1:204325/139‑1 (MQ=255)
                                                                                                           gagTTGCGCTGGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCAGAGAGTATCAATGCGGCACCGACCACGCGTCa     <  2:492628/139‑1 (MQ=255)
                                                                                                             gTTGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCAGAGAGTATCAATGCGGCACCGACCACGCGTCACTa  >  2:1462447/1‑139 (MQ=255)
                                                                                                              ttGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAg                                                      <  1:1348347/86‑1 (MQ=255)
                                                                                                              ttGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAg                                                      >  2:1348347/1‑86 (MQ=255)
                                                                                                              ttGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCa                                     >  1:1409234/1‑103 (MQ=255)
                                                                                                              ttGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCa                                     <  2:1409234/103‑1 (MQ=255)
                                                                                                              ttGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCAGAGAGTATCAATGCGGCACCGACCACGCGTCa     >  1:165059/1‑135 (MQ=255)
                                                                                                              ttGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCAGAGAGTATCAATGCGGCACCGACCACGCGTCa     <  2:165059/135‑1 (MQ=255)
                                                                                                                     |                                                                                                                                    
GCCTACTGGTGCAGGCAGGTGATCAGAGCGAGCAAATTTTTACCCGCTGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCAGAGAGTATCAATGCGGCACCGACCACGCGTCACTA  >  NC_000913/4329772‑4330018

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: