Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 665,554 C→T E254K (GAA→AAA)  mrdB ← cell wall shape‑determining protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913665,5540CT98.2% 207.8 / ‑3.4 56E254K (GAA→AAA) mrdBcell wall shape‑determining protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (24/31);  total (24/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

CCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTCAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACGCGCTGGCGCTGGT  >  NC_000913/665406‑665697
                                                                                                                                                    |                                                                                                                                               
ccAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTaa                                                                                                                                               <  1:846080/151‑1 (MQ=255)
      gtggtTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCtgtg                                                                                                                                         >  2:1132774/1‑151 (MQ=255)
        ggtTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGAc                                                                                                                                       <  1:513260/151‑1 (MQ=255)
              cgcTCTGGCGCGTATCCACAGCCCGTGCATGCTCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGt                                                                                                                                  >  1:916471/1‑150 (MQ=255)
               gcTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGc                                                                                                                                >  2:702365/1‑151 (MQ=255)
                ctctGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTg                                                                                                                                 <  2:793296/149‑1 (MQ=255)
                    ggcggcTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTg                                                                                                                            <  1:222264/150‑1 (MQ=255)
                        gcTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGcagcc                                                                                                                       >  1:294855/1‑151 (MQ=255)
                          tATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGcagcc                                                                                                                       >  1:848939/1‑149 (MQ=255)
                               acaGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCttt                                                                                                                <  2:559502/151‑1 (MQ=255)
                                caGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCttt                                                                                                                <  1:1109574/150‑1 (MQ=255)
                                    ccGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCcgc                                                                                                           >  1:171412/1‑151 (MQ=255)
                                    ccGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCcgc                                                                                                           <  1:655873/151‑1 (MQ=255)
                                        gcATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGGAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTa                                                                                                        >  1:174457/1‑150 (MQ=255)
                                        gcATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTa                                                                                                        >  2:440310/1‑150 (MQ=255)
                                          aTGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCGCTAACCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATc                                                                                                     <  2:190790/151‑1 (MQ=255)
                                          aTGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAAt                                                                                                      >  2:132341/1‑150 (MQ=255)
                                          aTGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATc                                                                                                     >  1:1142461/1‑151 (MQ=255)
                                           tGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATcc                                                                                                    >  1:614896/1‑151 (MQ=255)
                                             aTCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATccgc                                                                                                  <  2:580848/151‑1 (MQ=255)
                                                agcagAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCgg                                                                                               <  1:25199/151‑1 (MQ=255)
                                                     aaTGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCa                                                                                          <  1:634265/151‑1 (MQ=255)
                                                        gTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATa                                                                                       <  1:166477/151‑1 (MQ=255)
                                                         tAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAAt                                                                                        >  1:416706/1‑149 (MQ=255)
                                                         tAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATa                                                                                       <  2:467781/150‑1 (MQ=255)
                                                          agagagCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGc                                                                                     >  1:973421/1‑151 (MQ=255)
                                                          agagagCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGc                                                                                     <  2:969639/151‑1 (MQ=255)
                                                             gagcgagcAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAAt                                                                                  <  2:238541/151‑1 (MQ=255)
                                                                   gcAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAg                                                                             <  2:358963/150‑1 (MQ=255)
                                                                     agaatcagaatGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACt                                                                          <  2:139273/151‑1 (MQ=255)
                                                                        atcagaatGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAAACTGaa                                                                       <  2:299128/151‑1 (MQ=255)
                                                                            gaatGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTCAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGaataat                                                                   <  2:531136/151‑1 (MQ=255)
                                                                                 cccACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGAt                                                              >  1:257428/1‑151 (MQ=255)
                                                                                  ccACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATa                                                             <  2:1010095/151‑1 (MQ=255)
                                                                                   cACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGAt                                                              >  2:496575/1‑149 (MQ=255)
                                                                                     cTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAg                                                            <  1:784857/149‑1 (MQ=255)
                                                                                       aaTCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGccc                                                         <  2:1141158/150‑1 (MQ=255)
                                                                                        aTCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCcgc                                                       >  2:776581/1‑151 (MQ=255)
                                                                                           ccAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGcc                                                    <  2:130725/151‑1 (MQ=255)
                                                                                           ccAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGcc                                                    <  2:861937/151‑1 (MQ=255)
                                                                                               ctctTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGt                                                >  2:766799/1‑151 (MQ=255)
                                                                                                tctTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTg                                               >  1:103318/1‑151 (MQ=255)
                                                                                                     cGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCt                                          <  1:1050842/151‑1 (MQ=255)
                                                                                                     cGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCt                                          <  1:845007/151‑1 (MQ=255)
                                                                                                         aGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGAtt                                      >  1:927074/1‑151 (MQ=255)
                                                                                                           tACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTcc                                    >  2:565962/1‑151 (MQ=255)
                                                                                                            aCCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCg                                   <  2:803067/151‑1 (MQ=255)
                                                                                                                  aaGATAAAGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTcc                              >  1:1062223/1‑150 (MQ=255)
                                                                                                                       aaaGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGc                        <  1:991017/151‑1 (MQ=255)
                                                                                                                       aaaGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGc                        <  1:123020/151‑1 (MQ=255)
                                                                                                                       aaaGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGc                        <  1:1086310/151‑1 (MQ=255)
                                                                                                                        aaGTCAGTATGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGc                        >  2:767238/1‑150 (MQ=255)
                                                                                                                                aTGGCGTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTAcg               <  2:1090880/151‑1 (MQ=255)
                                                                                                                                     gTTCGGGGAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACgcgctg          <  2:257966/151‑1 (MQ=255)
                                                                                                                                          gggAGAAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACgcgctggcgct     >  1:126400/1‑151 (MQ=255)
                                                                                                                                              gaAATTTAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACGCGCTGGCGCTGGt  >  1:515037/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
CCAAAGGTGGTTTGCGCTCTGGCGGCTATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTCAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACGCGCTGGCGCTGGT  >  NC_000913/665406‑665697

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